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Summary 

From  the  equations  of  niotion  of  a  rotating  rocket,  closed  solutions 
are  obtained  for  the  ballistic  dispersion  dizring  boost  under  the  assumption 
that  the  angular  acceleration  is  constant  -  Part  II  and  Appendix  A. 

These  solutions  are  evaliuated  numerically  for  the  two  cases: 

(a)  launching  spin  zero 

(b)  angular  acceleration  zero, 

and  the  results  are  presented  and  discussed  in  Part  I.  By  comparing  the 
dispersion  v^'ith  that  of  a  non-rotating  round,  the  effectiveness  of  spin 
as  a  means  of  reducing  dispersion  can  be  assessed  in  each  case. 

In  case  (a)  it  is  found  that  reductions  in  the  dispersion  by  a 
factor  of  3  arc  theoretically  possible  by  of i -setting  the  nozzles  of  a 
multiple  boost  system  tangentially  by  less  than  5^-  To  achieve  larger 
reductions  by  this  method,  the  nozzle  off-set  must  be  increased  in 
proportion  to  the  square  of  the  reduction  sought. 

In  case  (b) ,  the  constant  spin  case,  it  is  found  that  a  spin  of 
Just  less  than  1  rev/ sec  is  in  general  sufficient  to  reduce  the  dispersion 
by  a  factor  of  3-  For  higher  spins  the  reduction  varies  linearly  with 
the  spin. 


In  Part  III  closed  solutions  are  obtained  for  the  dispersion  of  the 
dart  duo  to  unclean  separation  and  aerodynamic  malalignments.  Prom  thoso 
solutions  the  maximum  dispersions  tnat  can  arise  are  then  deduced.  The 

results  are  presented  in  Part  I,  and  it  is  concluded  that  dispersion  of 
the  dart  from  these  causes  will  be  small  compared  v^ith  the  dispersion  at 
the  end  of  boost. 
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PART  I 

THE  EFFECT  OF  SPIN  ON  DISPERSION  DURING 

EOOST  AND  AFTER  SEPARATION 


1  Introduction 

■ - ...  ■■  ■  ■  ■  ■ 

1*1  General  remarks 

The  trajectory  followed  by  a  guided  missile  during  the  initial 
uncontrolled  pai^t  of  its  flight  cannot  be  predicted  with  certainty  on 
account  of  a  number  of  effects  of  unknown  magnitude  that  influence  its 
motion,  e.g.  thrust  variation  of  boosts.  The  trajectories  of  a  large 
number  of  similar  rounds,  i.e,  rounds  manufactured  to  the  same  design, 

Y/ill  therefore  give  a  certain  distribution  of  trajectories  about  the 
mean  (or  standard)  trajectory.  A  measure  of  the  *  spread*  of  this 
distribution,  e^g.  the  standard  deviation  or  the  radius  of  the  95^. zone, 
is  called  the  dispersion  of  the  round. 

IVhcn  the  dispersion  is  known  the  percentage  of  rounds  whose 
trajectories  are  aeeeptable  to  the  particular  guidance  system  considered 
can  be  determined;  or  alternatively,  the  requirements  of  a  particular 
guidance  system*  can  be  stated  for  which  only  a  certain  small  percentage 
(say  5/^)  of  the  missiles  is  lost.  In  this  respect  some  guidance  systems, 
such  as  a  non-directional  command  link,  are  ideal  because  they  could  be 
made  to  operate  for  all  positions  of  the  round  in  space.  At  the  other 
extreme  there  is  beam- riding  guidance,  which  operates  only  as  long  as 
the  missile  lies  within  the  narrow  cone  of  the  radar-tracking  beam.  In 
this  case  an  ancillary  guidance  systeni  is  necessary  to  shepherd  the 
rounds  into  the  tracking  beam.  The  larger  the  angle  of  this  gathering 
beam  hov/ever  the  longer  the  transfer  time  of  the  missile  to  the  narrow 
beam  and  hence  the  greater  the  minimum,  engagement  range.  On  this 
account  it  is  important  to  utilise  a  gathering  beam  whose  angle  is  as 
small  as  possible.  This  demands  a  knov;iedge  of  the  dispersion  of  the 
round  and  the  relative  importance  of  the  various  factors  upon  v/hieh  the 
dispersion  depends. 

In  Ref.  3  the  contributions  of  these  various  factors  to  the  total 
dispersion  are  discussed  for  Seaslug:  the  ballistic  dispersion  is  found 
to  be  the  most  important  single  factor  upon  which  the  width  of  the 
gathering  beam  will  depend.  This  fact  has  received  recent  support 
from  firings  of  302/STV*s:  a  random  dispersion  of  about  8*^  for  the  95!^ 
zone  was  observed.  (This  value  is  higher  than  was  expected  and  suggests 
that  previous  estimates  of  the  ballistic  dispersion  for  a  wrap-round 
boost  arc  in  error  by  a  factor  of  2  or  3.  There  is  evidence  however  that 
the  thrust  variation  for  7^”  boost  motors  is  greater  than  previously 
thought,  and  this  may  v/ell  explain  the  discrepancy  between  the  estimated 
and  observed  ballistic  dispersions.) 

Ref.  2  deals  with  the  ballistic  dispersion  of  various  non-spinning 
boost  configurations  and  shows  that  for  a  given  layout  the  most  profitable 
method  of  reducing  the  ballistic  dispersion  during  boost  is  to  increase 
either  the  aerodynamic  stability  of  the  round  or  the  length  of  launcher 
from  which  it  is  fired.  In  practice  there  is  clearly  a  limit  to  which 
such  steps  can  be  taken  without  introducing  formidable  constructional 
problems.  In  addition  the  theory  shows  that  beyond  a  certain  point 
these  measures  become  less  and  less  effective.  Also,  the  gain  from 
increased  stability  will  be  off-set  in  the  limit  by  certain  dispersions, 
e,g,  wind  error  dispersion,  which  increase  with  the  stability  of  the 
configuration. 
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A  further  means  of  reducing  dispersion  that  v/arrants  consideration 
is  that  of  spinning  the  round.  It  would  he  expected  that  even  a  small 
spin  amounting  to  one  revolution  v/ithin  the  first  half -wavelength  of  yaw 
v/ould  prevent  the  accumulation  of  malalignment  dispersions  in  any  one 
direction.  This  statement  presupposes  that  the  malalignments  themselves 
remain  independent  of  the  spin.  It  must  he  assumed^  for  instance,  that 
with  a  tandem  boost  the  ^unknown  factor*,  i*e.  the  factor  between  the 
observed  dispersion  and  that  predicted  from  the  measured  angular  displace¬ 
ment  of  the  nozzle  axis  from  the  rocket  axis,  is  the  same  for  a  rotating 
round  as  for  a  non-rotating  round.  This  point  has  received  practical  ' 
confirmation  from  a  series  of  firings  of  3**  rockets^. 

The  introduction  of  spin  into  the  problem  makes  the  theoreticaJL 
treatment  much  more  complicated  and  leads  to  integrals  which  are  not  at 
all  suitable  for  numerical  evaluation  on  account  of  their  highly 
oscillatory  integrands.  Approximate  evaluations  of  these  integrals 
for  the  ballistic  dispersion  have  been  obtained  in  Ref,  1  in  certain 
special  cases,  and  depend  on  the  launching  spin  being  sufficiently  large. 
This  form  of  solution  -  an  asymptotic  solution  -  is  modified  in  the 
present  report  (App.  A)  to  cover  other  cases  of  possible  guided-weapon 
interest  with  somewhat  improved  accuracy.  In  order  to  deal  with  zero 
launching  spin  -  the  case  most  relevant  to  present  guided-weapon  work  - 
a  nov/  solution  has  been  derived  that  is  exact 'for  a  neutrally  stable  round 
and  approximately  correct  in  general  for  sufficiently  large  rates  of 
change  of  spin  at  launch.  The  discovery  of  this  solution,  called  the 
approximate  solution  in  App,  A,  v/as  one  of  the  most  difficult  obstacles 
to  be  overcome  before  a  complete  study  of  the  effect  of  spin  could  be 
undertaken.  IWien'  the  spin  and  rate  of  change  of  spin  at  launch  are 
both  large  the  asymptotic  and  approximate  solutions  arc  in  agreement. 

The  evaluation  of  the  dispersion'  of  the  dart  after  separation  has 
also  required  a  great  deal  of  new  work,  involving  the  formulation  of  the 
equations  of  motion  \7ith  the  inclusion  of  the  aerodynamic  malalignments  in 
Part  II  and  their  evaluation  and  simplification  in  Part  III.  Part  1 
presents  the  .results  obtained  from  the  theoretical  investigations  of 
Part  II  and  Part  HI,  quoting  the  appropriate  formulae. 

It  has  not  been  necessary  to  lay  down  any  definite  missile  configur¬ 
ation  but  only  to  consider  the  values  over  v/hich  certain  fundamental 
parameters  are  likely  to  range.  Results  are  given  for  a  wide  range  of 
these  parameters  and  so  do  not  relate  to  any  specific  design,  although 
their  application  to  a  Seaslug-type  missile  is  considered. 

It  has  been  assumed  throughout  that  the  round  has  symmetry  of  order 
three  or  more  (as  defined  in  Ref.  1.5)* 

1.2  Methods  of  imparting  spin  during  boost 

No  detailed  consideration  is  given  in  this  note  to  the  various 
ways  and  means  by  which  spin  could  be  imparted  to  a  round.  Instead 
results  are  obtained  for  two  types  of  spin-form,  to  vdiich  most  spin-forms 
occurring"  in  practice  arc  likely  to  approximate. These  are 

(a)  constant  angular  acceleration  withi  zero  rate  of  rotation  at 

lauhch,  and 

(b)  constant  rato  of  rotation  after  launch. 

Type  (a)  can  be  identified  closely  with  the  spin  arising  from  off-set 
boost  nozzles  of  a  multiple  boost  system,  and  typo  (b)  with  the  spin  of  a 
round  projected  from  a  spiral  or  rotating  launcher,  particularly  during 
the  first  part  of  the  boost  period.  The  method  of  spinning  by  off- sotting 
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boost  nozzles  has  been  applied  to  test  vehicles  with  wrap-round  boosts 
and  is  the  only  method  in  use  at  the  present  time.  This  method  might 
also  be  applicable  to  a  tandem  boost  with  multiple  nozzles,  or  to  a  boost 
of  the  bircumferential  type  (for  description’  see  Ref.  2),  particularly 
if  the  stabilising  fins  are  offset  to  delay  their  damping  out  of  the 
spin  until  later  in  the  boost  period. 

The  method  of  imparting  spin  by  off-set  fins  alone  is  unsatisfactory, 
because  the  rate  at  which  the  spin  builds  up  is  then  smallest  when  the 
velocity  is  low,  and  hence  '//hen  the  largest  dispersion  occurs:  besides 
this  disadvantage  the  structural  requirements  of  the  fins  would  become 
severe  with  increasing  size  and  off-set  angle.  Rounds  with  large 
stabilising  fins  often  generate  a  small  spin  of  1  or  2  revs/sec  due  to 
a  slight  deformation  in  the  fins,  but  this  spin  does  not  build  up  until 
near  the  end  of  the.  boost  and  so  does  not  produce  any  effective  reduction 
in  the  dispersion. 

The  spin-forms  arising  from  off-set  boost  nozzles,  and  spin  at  launch 
and  combinations  of  both  are  the  subject  of  para,  2.  The  maximuiri  value 
that  the  spin  attains  during  boost  can  be  related  to  the  damping  properties 
of  the  round  in  spin  and  to  the  magnitude  of  the  boost  couple. 

1.3  Causes  of  dispersion  during  boost 

The  various  causes  of  dispersion  can  be  divided  into  two  classes: 

(a)  malalignments  inherent  in  the  round,  and 

(3)  external  causes,  e.g.  initial  launching  conditions  and  wind. 

The  dispersions  due  to  (3)  depend  on  the.  spin  only  in  so  far  as  it 
produces  a  processional  motion.  For  sufficiently  small  spin  this  effect 
is  negligible:  the  exact  criteria  are  contained  in  assumptions  B,1  and 
B.2  page  1^%  and  require  .that 


p2r2  <  <  nV 


and 


|3r 


<  < 


nV/s 


where  23  =  M  of  I  in  roll/ivi  of  I  in  pitch 

r  =  spin 

V  =  velocity 

a  =  acceleration 

^2y2  _  aerodynamic  restoring  moiaent/incidence 

M  of  I  in  pitch 


These  conditions  almost  certainly  hold  for  spins  of  the  magnitude 
arising  in  practice  with  guided  v/eapons,  say  of  less  than  5  revs/sec. 

When  the  precessional  effects  are  neglected,  the  motion  becomes  identical, 
with  that  of  the  non-rotating  round  and  so  is  not  considered  in  this  note. 
See  Ref ,  1,1, 
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The  dispersions  due  to  (pc)  on  the  other  hand  are  all  reduced  by 
spin,  except  for  the  dispersion  arising  from  malalignment  of  the  princi¬ 
pal  axis  of  inertia,  v/hich  increases  v;ith  increasing  spin.  It  is 
therefore  very  necessary  to  form  some  idea  of  the  value  of  the  spin  for 
which  the  inertia  axis  dispersion  becomes  comparable  with  the  remaining 
dispersions  due  to  (a),  as  no  advantage  is  gained  with  higher  spins. 

1.4  Dispersion  due  to  malalignment  of  the  principal  axis  of  inertia 

The  effect  of  a  malaligned  inertia  axis  on  the  motion  of  a  rocket 
is  equivalent  to  that  of  a  destabilising  couple.  If  ccq  is  the  small 
angle  made  by  the  principal  longitudinal  axis  of  inertia  with  the  axis 
of  symmetry,  the  component  of  angular  momentum  in  the  transverse  plane 
has  magnitude  AraQ,  v/here  A  is  the  M  of  I  in  pitch,  and  lies  in  the 
plane  containing  the  longitudinal  axes  of  symmetry  and  inertia.  The 
rate  of  change  of  this  vector,  v/hich'is  equal  to  the  couple  produced, 
has  magnitude  At^oq  and  lies  in  the  transverse  plane  in  the  destabilising 
sense.  For  a  round  v/ith  constant  spin,  i.e.  type  (b),  the  inertia 
axis  dispersion  is  therefore  identical  v/ith  that  of  a  constant  destabilising 
couple  I^aQ^  and  can  be  readily  compared  vriLth  that  due  to  a  constant 
boost  destabilising  couple  of  magnitude  jGp  j  •  The  dispersions  from  these 
two  causes  are  then  equal  when  the  spin  is  r*^^  where 


1^1  /  “C  • 


It  is  shown  later  in  para.  5*13  that  the  dispersion  due  to  a  constant 
destabilising.,  moment  varies  inversely  v/ith.  r  (for  values  of  r  not 
near  zero).  It  can  then  be  proved  that  'when  r  =  r^  the  resultant 
of  these  twe  dispersions  is  actually  a  minimum.  This  result. was  found 
to  be  in  reasonable  agreement  with  the  dispersions  observed  in  the  series 
of  firings  of  3”  solid  fuel  rockets  from  spiral  projectors^.  A  minimum 
dispersion  was  obtained  with  a  launcher  pitch  of  8  feet,  giving  a 
launching  spin  of  I6  revs/sec. 

If  h  is  the  moment  arm,  lGp|/thrust,  r^  can  be  expressed  by 


ro^  -  ha  /  kp2  , 


where  is  the  radius  of  gyration  in  pitch.  For  a  round  with  a  tandem 
boost  ni  is  roughly  equal  to  the  thrust  malalignment  angle  times  the 
distance  of  the  boost  nozzle  exit  plane  from  the  C;G.  Hence  h  is 
approximately  V'3  kp  and 


« 

r^^  =  V3  a  ttj  /  kp  aQ  . 


The  values  of  ajj.ttQ  for  the  35yo  zone  obtained  from  the  3"  rocket  firings 
were  O.OOi+B  and  0.00075  radians.  If  these  values  are  taken  to  be 
appropriate  to  a  guided  weapon,  the  value  of  ro  for  a  round  of  30  ft 
length  and  600  ft/sec^  acceleration  is  then  l+.h-  revs/sec. 
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For  a  ndssile  with  a  wrap-round  boost  h  and  depend  on  the 
positioning  of  the  boosts  along  the  body.  From  values  given  in  Ref.  2 
for  a  t3^ical  V/rap-ro.und  configuration  v^ith  a  10*^  venturi  radial  off-set 
(see  Figs.  3  and  6  in  Ref.  2),  h/kp^  is  found  to  be  about  0,00077  ft”"^ 
for  all  boost  positions  considered.  Taking  =  0,00075  radians  and 
a  =  600  ft/seo2,  the  above  equation  gives  2^-,0  revs/sec  for  r^. 

The  dispersion  arising  from  any  malaligrjneiit  is  clearly  greatest 
when  the  velocity  is  least,  i.e,  at  the  beginning  of  the  boost  period. 

It  is  therefore  necessary  to  limit  the  spin  only  during  this  initial 
period.  In  practice  this  amounts  to  a  limit  on  the  launching  spin  and, 
to  a  lesser  extent,  on  the  initial  angular  acceleration.  The  value 
of  4  revs/sec  estimated  above  for  this  limit  is  large  enough  to  show  that 
the  use  of  spin  merits  further  consideration:  the  object  of  this  note 
is  then  to  find  out  how  effectively  the  remaining  dispersions  due  to  (a) 
are  reduced  by  smaller  Initial  spins.  In  view  of  the  uncertainty  of  this 
value  we  proceed  under  the  assumption  that  for  the  values  of  spin  under 
consideration  the  dispersion  Irom  the  inertia  axis  malalignment  is 
negligibly  small  compared  with  the  other  dispersions  due  to  (a).  The 
reductions  in  dispersion  recorded  may  therefore  be  slightly  optimistic. 

1.5  Dispersions  due  to  boost  and  aerodynamic  malalignments 

The  principal  remaining  dispersions  arise  from 

A.  destabilising  moments  due  to 

(i)  boost  malalignment 

(ii)  aerodynamic  malalignment^  and 

B.  transverse  forces  due  to 

(i)  boost  malalignment 

(ii)  aerodynamic  malalignment. 

During  the  boost  period  the  magnitudes  of  A, (l)  and  B.(i)  will 
vary  in  some  unknown  manner  with  irregularities  in  the  boost  thrust. 

If  however  a  large  number  of  firings  is  taken  into  consideration  A.(i) 
and  B.(i)  can  be  replaced  by  constant  values,  depending  on  their 
statistical  distribution  for  'che  batch  of  rounds,  and  the  resulting 
dispersions  can  be  evaluated. 

The  variation  of  A.  (ii)  and  B.(ii)  during  boost  is  more  definite, 
provided  acro-olastic  effecuS  arc  neglected,  being  proportional  to  the 
square  of  the  velocity.  The  resulting  dispersions  can  then  be 
represented* in  terms  of  integrals,  which  are  however  not  soluble  in 
any  convenient  ‘closed  form*;  and  numerical  evaluation  would  be  tedious 
on  account  of  the  highly  oscillatory  integrands.  In  practice  it  should 
be  possible  to  keep  the  aerodynamic  malalignments  to  within  a  limit  for 
which  the  resulting  dispersions  arc  small  compared  vd.th  the  ballistic 
dispersions  A.  (i)  and  D.(i).  In  Ref.  2  it  has  been  shown  that  in  the 
case  of  an  unspun  round  the  dispersion  from  A.(ii)  will  be  small  compared 
with  that  from  all  other  causes  provided  the  aerodynamic  malalignment 
angle  is  kept  well  belovv^  1°.  -The  introduction  of  spin  will  reduce 
aerodynamic  dispersions  more  than  ballistic  dispersions,  because  A, (ii; 
and  B.(ii)  are  very  much  smaller  at  launch,  when  the  spin  might  stoll 
be  building  up,  than  at  the  end  of  boost  (say  100  times  greater)  when  the 
spin  is  established.  For  this  reason  a  more  detailed  examination  of 
the  aerodynamic  dispersions  has  not  been  undertaken. 
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The  boost  malalignment  dispersions  A.(i)  and  B. (i)  are  considered 
in  para.  4  for  the  case  of  constant  angular  acceleration  with  zero 
launching  spin  (type  (a)),  emd  in  para.  5  for  the  case  of  constant  spin 
(type  (b)).  It  is  assumed  that  the  launch  is  perfect^  i.e.  that  the 
malalignments  do  not  alter  the  initial  conditions  at  launch.  By 
launching  a  round  already  spinning  much  greater  reductions  in  dispersion 
are  obtainable  than  by  spinning  only  after  launch  within  practical 
limits.  Since  haviever  certain  dispersions  are  not  reduced  by  spin, 
e.g.  dispersion  due  to  v/ind  error,  there  is  a  limit  to  the  reduction 
it  is  worthwhile  trying  to  achieve.  This  point  is  illustrated  in  detail 
in  para.  6  for  the  case  of  a  Seaslug  type  missile.  If  this  limit  can 
be  attained  with  say  off-set  nozzles,  then  no  advantage  is  gained  by 
launching  with  spin  from  say  a  spiral  or  spinning  launcher.  Besides 
this  point,  it  is  likely  that  in  most  cases  the  structural  problems 
associated  with  a  spiral  launcher  would  militate  against  its  introduction; 
such  considerations  as  these  are  not  taken  to  be  v/ithin  the  scope  of  this 
note. 

1.6  Dispersions  arising  after  separation 

In’  paras,  7  and  8  consideration  is  given  to  the  dispersions  arising 
after  boost  separation  and  before  the  round  comes  under  control.  The 
object  of  these  paragraphs  is  to  deduce  limits  Y/ithin  v/hich  the  aero¬ 
dynamic  malalignments  of  the  dart  must  lie  if  large  dispersions  are  to 
be  avoided:  dispersions  arising  from  ^unclean*  separation  are  also 
discussed.  The  method  of  approach  naturally  requires  a  number  of 
approximating  assumptions  to  be  made:  it  is  assumed,  for  example,  that 
the  spin  is  constant.  This  assuraption  is  hov^-ever  justifiable,  since 
in  theory  the  spin  v/ill  quickly  tend  to  a  constant  value,  depending 
on  the  malaligniiient  incidences  of  the  aerodynamic  surfaces.  The  other 
important  parameters  which  determine  thj  motion  are  the  lift  and  stability 
properties  of  the  round. 

It  is  found  that  should  the  spin  be  zero,  the  dispersions  can 
increase  indefinitely,  and  so  a  limit  must  be  imposed  on  the  time  of 
uncontrolled  flight  if  large  dispersions  are  to  be  avoided.  The 
presence  of  even  a  small  spin  reduces  the  dispersions  appreciably. 


2  Spin- form 

2.1  Two  methods  of  imparting  spin 

Spin  may  be  imparted  to  a  round  by  either  or  both  of  the  following 
methods. 

(a)  By  off-setting  the  noz^es  of  a  multiple  boost  assembly, 

(b)  By  projecting  the  round  from  a  spiral  launcher. 

By  method  (a)  the  round  is  angularly  accelerated  from  the  instant  it  ceases 
to  be  constrained  by  the  launcher,  and  continues  to  accelerate  iintil  the 
aerodynamic  spin  damping  moment  builds  up.  During  the  latter  half  of 
the  boost  stage  the  spin  will  in  general  decrease  slightly  according  to 
the  magnitude  of  this  damping  moment,  'i/Wien  the  decrease  in  boost  thrust 
sets  in,  the  spin  decreases  rapidly.  With  method  (b)  the  spin  is 
appreciably  constant  just  after  launch,  but  decreases  Yri.th  increasing 
velocity  and  damping  moment. 

Since  the  dispersion  is  most  affected  by  the  spin-form  immediately 
after  launch,  it  is  desirable  that  the  spin  should  build  up  as  quickly 
as  possible.  In. practice  however  the  angular  acceleration  and  rate  of 
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rotation  are  limited  for  structural  reasons.  With  a  spiral  launcher  the 
spin  is  of  course  greatest  at  launch,  while  the  angular  acceleration  when 
on  the  launcher  is 


(1) 


where  and  are  the  spin  and  velocity  at  launch  and  a  is  the 
acceleration.  T.his  equation  shows  that  h  is  proportional  to  a  and  so 
will  be  roughly  constant  in  practice^  reaching  a  maximuin  by  launch. 

The  spin-forms  produced  by  methods  (a)  or  (b)  or  combinations  of 
both  are  derived  in  the  following  sub-paragraphs, 

2.2  Solution  for  spin 

2.21  Prom  the  equation  of  motion  in  spin,  namel^y, 


Gr  ^  Tr  =  G,  (2) 


where  G-  is  the  magnitude  of  the  boost  couple 
•  C  is  the  M  of  I  in  roll 

and  is  the  magnitude  of  the  aerodynamic  damping  moment  in  spin 

defined  in  11.3.22*^, 

we  see  that  at  launch  when  Pr  is  small  compared  to  G-,  the  initial 
acceleration  is 


=  G/G, 


(3) 


and  that  until  Pj^  becomes  comparable  with  G,  the  spin  equation  is 


r  =  r^  +  (g/g) (t  -  tg) , 


(4) 


where  t^  is  the  time  at  the  instant  of  launch. 

2.22  To  solve  for  the  spin  over  the  whole  boost  range,  it  is 
assumed  that 


%  Y  pVr  , 


R'' 


(5) 


read  as  Part  II  paragraph  3.22. 

11. 


CONFIDENTIAL 


CONFIDENTIAL 


Technical  Note  No.  G.V/'.  177 


where  is  constant  over  the  range  v/here  is  comparable  v/ith  G. 

This  assumption  has  been  justified  by  an  aji^ysis  of  S.T.V.l  spin-forms. 
Prom  three  similar  spin-forms  the  maximum  value  attained  by  the  spin 
was  read  off,  and  Yj^  found  by  the  relation 


■Cr  =  „ 

r=0 


Using  this  value  of  Yr  (0,30  lb  ft)  and  taking  into  aeeount  the  fall-off 
in  G-  just  before  separation,  the  spin  obtained  by  integration  of  equation 
(2)  showed  complete  agreement  v/ith  the  firings. 

2.23  The  solution  of  equation  (2),  when  the  acceleration  is  assumed 
constant,  ean  be  written  non-dimensionally  as  (sec  II. 9.3) 


r/rG  = 


-T'' 


+  e 


[e(t)  -  E  (To)]  , 


(6) 


where  T^  '=  Y ps/C 
r(.  =  k  G/C 

X 

E(x)  =  j"  du 

o 

k  =  y2C/aYp 

and  the  time  t  is  given  by 


(7) 


(8) 


t  =  kT. 


(9) 


Por  launeher  lengths  of  the  order  of  those  occurring  in  practice,  the  terms 
in  Tq  in  this  eq’oation  are  negligibly  small,  on  aeeount  of  the  smallness 
of  the  aerodynamie  damping  at  launeh. 

2.21f  In  Pig.  1  all  curves  are  plotted  for  s^  =  0.  The  curve  A  is 


-t2 

r/r(.  =  e  E(t) 


and  so  denotes  the  spin-form  v/hen  r^  =  0.  This  equation  shov/s  that  the 
spin  at  any  instant  is  proportional  to  G  and  hence  to  the  nozzle 
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offset  angle  in  the  tangential  plane  Ajj.  The  maximum  value  reached  by  the 
spin  is  seen  from  the  graph  to  be 


=0.54  k  Fq  by  equations  (3)  and  (7),  (lO) 
and  occurs  at  T  =  0.93, 


i.e.  at  time 

t  =  0.93  k 

by  equation  (9) . 

Combining  this  curve  v/ith  any  of  the  B  curves, 


shown  for  r^/rQ.  equal  to  0,50^  0«75  Q-nd  1.00,  we  find  the  spin  obtained 
when  both  methods  (a)  and  (b)  are  used  together. 

Curves  B  determine  the  spin  obtained  by  method  (b)  alone,  namely. 


2.3  Estimation  of  spin  damping  moment 

2.31  In  order  to  determine  Yp  for  a  given  design  of  missile  it 
is  assTomed  that 

(i)  the  contribution  to  the  damping  moment  from  the  wings  (or 
fins)  is  very  much  greater  than  fran  any  other  part  of  the  body, 

and 

(ii)  the  lift  distribution  across  the  v/ings  (or  fins)  +  body  is 
elliptic. 

It  can  then  be  shov/n  that 


Yr  =  c  kR  (2di)^, 


(11) 


where  for  cruciform  wings 


1  - 


26, 

% 


sinifO^ 

2% 


/8 
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J.-  ^  * 

sin  0^  =  d^di 

(13) 

2d^  =  mximuiTi  ^body  +  boosts*  diameter 

2d^  =  wing-span 

and 

=  lift ,  constant  =  lif t/incidence  x 

(velocity) 

This  formula  should  give  good  results  at  subsonic  velocities  for  wings  of 
high  aspect  ratio.  For  very  low  aspect  ratios  however,  assumptions 
(ij  and  (ii),  particularly  (ii),  are  scarcely  justified  and  this  equation 
for  c  becomes  unreliable,  although  would  still  be  expressible  in  the 
form  of  equation  (ll) . 

For  the  S.T.V.l  equation  (l2)  was  found  to  lead  to  a  value  of  Yp 
in  reasonable  agreement  with  tha.t  determined  from  the  maximum  spin, 
see  para.  2.22.  The  wing-span  is  3  feet,  and  the  maximim  radial 
distance  swept  through  by  the  boosts  is  10.5  in;  taking  k-^,  obtained 
by  lov^-speed  v/ind-tunnel  tests^  to  be  0.2+8  Ib.ft"'^,  equation  (ll) 
gives 


=  0.38  lb. ft. 

2.32  An  approximate  expression  for  k,  defined  by  equation  (8), 
can  be  obtained  by  taking 


C  =  m  do^/2 


app  r  o  xinja  t  ely . 

Then  by  equations  (S)  and  (ll) 


k  =  f(0J//^ 


(14) 


v/here 


f(0o)  =V2  sin  0Q  /ir  -  20q/x  +  sin  24.0q/2J'; 

e  =  k^/m  (15) 

■  m  =  mass  of  projectile. 

f(eo)  is  plotted  in  Fig.  22. 

3  Dispersions  arising  from  boost  malalignments 

3.1  Definitions  and  notation 

The  angular  deviation  is  the  angle  betv/een  the  axis  of  the  launcher 
and  the  direction  of  motion  of  the  C.G.  of  the  round.  The  dispersion 
is  a  statistical  measure  for  the  angular  deviations  of  a  number  of  rounds, 
e.g,  the  root  mean  square  or  the  sise  of  the  95^  zone. 

The  angular  displacement  is  the  angle  at  the  launcher  between  its 
axis  and  the  direction  of  the  C.G.  of  the  round.  If  the  angiilar 
deviation  increases  monotonioally  with  distance  it  follows  that  the 
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angular  displacement  is  less  than  the  angular  deviation;  the  two  in 
fact  become  equal  at  -infinity.  (in  practice  the  angular  deviation  v/ill 
not  necessarily  be  a  mono  tonic  function,  but  it  can  usually  be  considered 
as  such  plus  oscillations  of  small  airplitude  which  do  not  invalidate 
this  result)  .  The  value  of  the  angular  deviation  at  any  point  therefore 
sets  an  upper  limit  to  the  width  of  beam  required  to  gather  a  round 
coming  under  control  at  that  point,  and  so  is  as  useful  a  quantity  as 
the  angular  displacement.  Besides  this,  the  integration  of  the  angular 
deviation  to  give  the  angular  displacement  can  bO  carried  out  only  in  a 
few  special  cases  and  so  in  general  is  not  readily  obtained.  No  further 
mention  will  be  made  of  the  angular  displacement  in  this  note. 

For  a  non- rota  ting  round  the  motion  lies  in  a  fixed  plane,  through 
the  axis  of  the  launcher,  and  so  the  angular  deviation  can  be  denoted 
by  a  single  quantity.  liVhen  however  the  round  is  rotating  the  motion 
can  be  resolved  on  to  tv/o  fixed  perpendicular  planes  through  the  launcher 
axiSp  and  tv/o  quantities  are  needed  to  define  the  angular  deviations  in 
these  tv/o  pianos.  If  these  quantities  fom  the  real  and  imaginary 
parts  of  a  complex  number  Z,  the  total  angular  deviation  is  approximately 
|Z|  fer'  small  angles  and  is  independent  of  the  orientation  of  the  two 
reforanee  planes. 

If  Z/p,  Zq.  are  the  complex  angular  deviations  duo  to  raalalignod 
boost  force  and  malaligned  boost  couple  respectively,  wo  can  write 


Zrtj)  =  PipZQi(s)  (l6) 

(17) 


Z^(s)  and  Zq.(s)  are  functions  of  the  range  s,  and  of  certain  parameters 
defining  the  aerodynamic  and  ballistic  properties  of  the  missile:  they 
are  dimensionless  quantities  obtained  by  solving  the  equations  of 
motion.  p  p  and  Pq.  are  constants  proportional  to  tlie  malalignment 
angles,  and  for  convenience  include  all  those  parameters  of  the  round 
not  occurring  inseparably  in  Z^(s)  and  Zq.(s).  and  are  of  course 

complex  quantities  whose  arguments  determine  the  orientation  about  the 
missile  axis  of  jTpJ  and  |Gp  |,  the  components  of  the  boost  thmst  and 
boost  couple  perpendicular  to  the  missile  axis.  The  expressions  for 
and  Imq.1  used  here,  and  to  v/hich  the  values  of  Zp(s)  and  ZQ.(s)/n 
given  later  correspond,  are 


1 

=  I'lplA 

(18) 

=  lGj,l/Tkp2, 

(19) 

whore  T  is  the  boost  thrust  and  kp  is  the  radius  of  gyration  of  the 
round  in  :*itchc 

It  is  Sometimes  found  convenient  to  express  |Tp|  and  jGpj  as  > 

jTpj  =:  T^pa^  (20) 

1Gp1  =  T6b«g,  (21) 
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v/here  Ag  is  the  offset  angle  of  each  nozzle  in  the  radial  plane 
is  the  distance  of  each  boost  thrust  axis  from  the  C.G. 
and  are  the  malalignment  angles  defined  by  these  equations. 

3.2  Parameters  in  problem 

In  order  to  obtain  any  reasonably  simple  closed  solution  for  the 
dispersion  functions  Zj(s)  and  Zq.(s),  quite  a  number  of  sitiq)lifying 
assumptions  have  been  introduced,  v/hich  leave  these  f'unctions  dependent 
on  only  the  few  most  influential  parameters.  The  spin-form,  for 
cxanq)le,  has  been  defined  to  the  extent  of  the  first  tv/o  terms  of 
its  series  expa.nsion,  i.e.  by  the  spin  at  launch  fo  and  the  angular 
acceleration  at  launch  r^.  This  is  justified  because  the  dispersion 
functions  tend  rapidly  to  a  limit.  It  is  found'  tliat  the  number  of 
parameters  can  then  be  reduced  to  four,  namely  n,  Sq,  nq  and  n2.  The 
linear  acceleration  occurs  in  nq  and  n2  and  is  not  an  additional  parameter 
in  the  dispersion  functions.  The  definitions  of  these  parameters  arc 
as  follows 

n  is  a  measure  of  the  stability  of  the  configuration  and  is 
defined  by 


n 


2  yS 


aerodynamic  restoring  raoment/incidence 
M  of  I  in  pitch 


(22) 


Sq  is  the  'effective  launcher  length'  defined  by 


^o 


(23) 


depend  on  the  spin  and  arc  defined  by 


nq  =  -  2  ^2^0 

(2A-) 

^2  =  ^c/a. 

(25) 

The  equation  for  the  spin  is 

^  =  ^o  +  ^o  (^  “  ■^o)  (26) 


and  corresponds  to  the  solution  obtained  in  para.  2.2  under  the  assumption 
that  the  aerodynamic  damping  moment  in  spin  is  negligibly  small.  Inte¬ 
gration  of  equation  (26)  leads  to  the  expression  for  the  total  angle  cr 
turned  through  by  the  round  about  its  axis, 

^(s)  -  <^o  =  ^2  (s  -  Sq)  +  nq  (/ s  -  V'sq).  (27) 
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From  the  straight-line  spin-forms  represented  by  equation  (26), 
two  sets  of  fan-ilies  v^ill  be  considered.  These  are  defined  by 

(a)  I'o  =  0 

(b)  =  0, 

and  correspond  respectively  to  methods  (a)  and  (b)  for  imparting  spin 
discussed  in  para,  2.1. 

Case  (a)  -  Launching;  spin  aero,  constant  angular  acceleration 

If.  1  Expressions  for  the  dispersion  functions 
4.11  2rp(s)  and  Zq(s) 

For  case  (a)  r^  is  zero  and  Zit*(s)  and  Zq(s)  depend  only  on  the 
three  parameters  n,  and  n2.  The  evaluation  of  these  functions  is 
presented  in  Part  II  para.  12o2  and  is  valid  only  for  small  values  of 
n/n2  (say  loss  than  “1);  these  results  are  now  quoted. 

V/e  can  write 


Zt(s)  =  nSjCs)  -  5t(s) 


Zq.(s)  =  nSQ.(s)  -  5q.(s) 


>  .  II. 12.1  (4), 


where* 


Sm(s)  =  -i 


f^J^(3o^^)  Jn(sQ,s) 

1  I - -  — - 


L  «i  ^l  J 


-  “1  3l 


O 

e 


11.11.3(11) 


^(s)  =  - 


I—  ttn 


+  1 


■p(0C,.  )  — 

-hiol  s^(s)  - 


—  CLi 


(s)j 


ICTf- 


11.11.3(12) 


5rn(  s)  — 


—  U  ^ 

•.  e  .  r 


2/3  L 


P(«s)  D(3s)' 

-  «1  ^  3l  J 


i[o-(s)-cro] 


N(aGQ)  -in(s-s-)  1^(3  Sq)  in(£5-SQ) 

L  3i  - 

.  11.10.2(8) 


ic-, 


'.Q 


^p(s)= 


©(ttg)  r>(3s) 


2i/3^  ^  ^1  ^1 


i[^(s)-o-J 


^(“o,)  -in(s-q)  D(3^  in(s-3-) 


u  a  - 


11.10.2(9), 


A  bar  over  a  quantity  denotes  its  complu>;  conjugate. 
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v/here 


le 


xYi  _ 


i[(r(u)-c'o]  D(YJ  -  D(Y 


u 


X 


for  X  =  a  or  3  , 


A.3. 32(39) 


where 


?^(s)  =  (2i/6^)  d(  )  -  D(6g)  e 


in(s-3Q) 


J...  11,11.2(10) 

(28) 


D(u)  =  e 


-i  ^  , 

2  r  1  +  i 


u 


i  ^ 

2  ' 
e  dx 


A.l.lfl 
and  (2 


(29), 


=  '/2(n2  +  n)/?^ 


3]_  =  \/2(n2  -  n)/'K 


Y^  =  2n2A 


=\/  2n/7C 


=  [(n„  +  n)  Vs  -  n2Vi7]  +  n)/2 


yii.io, 

2(6) 


Ys  =  \/2n2A  (Vs  -  V^) 


11.12.21(6) 


6^  =  \/  2nsA 


'^a=  - 


+  n 


7t 


11.11.2(7)- 

(9) 


_  11.12.21(7) 


Kft  = 


n 


2^2!o 


n2  -  n 


7!: 


11.12.21(8). 


As  s  increases  ttg,  3  g,  Yg,  6g  increase,  and  so  lD(ag)l,  lD(8g)l, 

|D(Yg)-|  and  lD(6g)  j  decrease  and  the  amplitudes  of  all  oscillatory  teras 
become  zero  at  infinity.  Prom  tables  of  D(u)  it  is  found  that  1d(u)1  has 

decreased  to  one-tenth  of  1  D(0)  |,  i.e.  of  l/V2,  when  u  =  4.5.  Ihe 
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amplitudes  of  the  rotational  terms,  i.e.  terms  in  are  there¬ 

fore  small  when  a^,  Pg,  and  Yg  are  greater  than  about  4.5.  When  n  is 
much  smaller  than  n2,  we  see  that  Pg  and  Yg  are  approximately  equal, 
and  are  greater  than  4.5  when 

Vs  ^  '^0  ^*5  \/'x/2n2 

^  +  5.6/Vn^. 

When  Sq  =  0,  n2  =  0.05  ft”^,  this  condition  requires  that  s  >  6l+0  ft. 

+in(s-So) 

The  amplitudes  of  the  terms  in  e  occurring  in  Zi'(sJ  and 

Zq^(s)  are  found  from  the  above  equations  to  decrease  like  |l^(5g)  j, 
where 


I)l(u) 


=  D(u). 

%u 


Prom  tables  it  is  seen  that  |D2_(u)  |  tends  to  zero  very  much  more  rapidly 
than  jD(u)|  as  u  tends  to  infinity,  and  equals  1d(0)1/10  at  about 
u  =  1,0.  This  means  that  the  effect  of  the  yawing  on  the  dispersion 
is  small  when  63  is  greater  than  1.0  i.e.  when 


s  >  n/2n. 


i.e,  after  the  first  half  wavelength  of  yaw.  Taking  n  =  O.OO5  ft"*^  this 
condition  requires  that  s  >  320  ft. 

It  has  been  shovm  in  this  paragraph  that,  for  most  practical  oases, 
Zj(s)  and  Zq.(s)  have  converged  sufficiently  near  to  Zj(eo)  and  '^n{oo)  after 
the  first  fev/-  hundred  feet  of  flight.  This  fact  is  illustratea  by  Pig.  9 
v^hich  shov/s  Zp(s)/n  for  the  typical  case  Sq  =  10  ft,  n  =  O.OO75  ft“^  and 
n2  =  0.05  ft”^.  The  corresponding  yav/  function  5Q.(s)/n  is  given  in  Pig, 10 
and  docs  not  converge  with  the  same  rapidity;  as  s  -♦  «>  ,  Sq.(s)  0. 

4.12  Zt(”)  and  Zq.(oo  ) 


The  equations  of  para.  4.11  simplify  somewhat  on  putting  s  =  <»  , 
Is  is  seen  that 


=  Cq.(“)  =  0  , 
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and  that 


5t(~)  = 


-  1 


Ja(so>~)  '^3(^0/°) 

rD(a3  )d(63  )  I)((3  )d(63  )-j 

0  0  0  0 

L  «l  Pi  J 

+  i 

L  “i^i  Pi^i  J 

(30) 


11.12.22(16) 


^(”)  =  - 


’Ja(so^”)  Jp(So,”)' 

''~h 


-J){%  )e(63  )  D(P  )d(63  ) 

_ -£L-  + _ -O— ^ _ CL. 


“1^ 


^  h 


icr^ 


(31) 


11.12.22(17), 


v/here  now 


Jx(so;»“)  =  -  i  ° 


[d(0)  -  H-\)]  / 


for  7\.  =  a  or  (3. 

y/hen  S3  =  0,  Vv'e  have  =0  ( ^  =  a  or  P)  and  hence 


Jx(o,-)  = 


1 0%. 


1  G 


D'(0)  / 


=  1 


1^0  / 
e  / 


t;  Yi- 


Since  Is  positive,  D(-Kp)  can  be  evaluated  from  tables  by  using  tho 
relation 


D('Kp) 


(1  +  i)  -  1(^3) . 


I..I3  Approximate  expressions  for  jZ^Cw)']  and  \Zq  (co )  1  /n 

i\l though  equations  (30)  and  (31)  can  be  evaluated  numerically  for 
given  values  of  the  parameters  n,  s^,  and  np  v^ithout  too  much  difficulty,, 
particiilarly  simple  expressions  for  Srji(“)  and  ^(“)  are  obtainable  from 
these  equations  for  sufficiently  large  values  of  n2.  It  is  shown  in 
11.12.23  that 


"^^■o  y  /  ^  (1  +  N 

Q  5^(co)  =  ^ -  M^s^) 

2\^n2n 


11.12.23(23) 
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^  2Vn^rr 


B(Sso) 


11.12.23(24), 


where  B(u)  and.  A(u)  are  the  real  and  imaginary  parts  of  D(u) 
It  is  then  easily  seen  that 


lZj(co)l  =  r)j^(co)l  =  7u/n/2n2  A(6g^) 

lZG.(“)l/n  =  =  7i;B(6g^)//2n2n 

for  large  n2. 


.  11.12.24(28) 

(32) 

.  11.12.24(29) 

(33) 


In  Table  I  the  values  of  1Zit(«>)1  and  |ZQ.(e3)l/n  given  by  equations 
(32)  and  (33)  are  compared  with  the  more  accurate  values  obtained  from 
equations  (30)  and  (3l). 


TABLE  .1 


n2 

n 

^0 

' accurate ' 
lZ,p(co)  1 

! 

A^So)/^^2^ 

•'  accurate' 

l^(“)i  /n 

ff^ 

ft“^ 

ft 

ft 

ft 

0.0025 

0 

0.25'  • 

’•  0.25 

99 

86 

30 

0.23 

0.23 

63 

54 

0.05 

0 

0.50 

0.49 

50 

37 

0.01 

30 

0.42 

0.37 

20 

13 

0.0025 

0 

0.18 

0.17 

70 

64 

0.10 

30 

0.17 

0.16 

45 

40 

0 

-  0.35 

0.35 

35 

30 

u,  u± 

30 

0.29 

0.27 

14 

11 

n  nno  c: 

0 

'0.14 

0.14 

57 

55 

30 

0.14 

0.13 

■  36 

33 

0.15 

0 

0.29 

0.29  . 

29 

26 

0.01 

30 

0.24 

0.23 

12 

10 

4.2  Dispersion  due  to  boost  destabilising  couple  -  Results 

In  Figs.  2  to  4  the  values  of  lZQ.(«)|/n,  obtained  from  the 
accurate  equations  of  para.  4.12,  are  plotted  against  n  in  the  range 
0.0025  <  n  ^  0.01  ft-1  for  Sq  =  D,  10,  20  and  30  ft  and  ■n2  =  O.O5,  0.10 
and  0,15  ft”^. 
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The  vailues  of  the  launching  velocities  corresponding  to  these  values 
of  Sq  are  0,  100,  1^)  and  170  ft/sec  when  the  acceleration  is  500  ft/soc^. 

The  values  of  boost  offset  angle  Ajj  corresponding  to  the  above 
values  of  n2  can  be  arrived  at  as  follov/s.  By  equations  (3)  and  (25) 


n^  —  Or/ Ca 


—  niAjq^j|/C, 


(34) 


where  is  the.  distance  of  the  centre  of  each  nozzle  exit  plane  from 
the  missile  axis,  and  m  is  the  mass  of  the  round.  Putting 


C  =  (roughly), 


equation  (30  'becomes 


n2  =  2Aj^ej^  .  •  (35) 

Vdien  equals  1  foot,  the  values  of  Aj^  corresponding  to  n2  =  0,05,  0,10 
and  0,15  ft“^  are  1.5°,  3.0”  and  4.5°  respectively. 

In  Pigs,  5  to  8  the  same  results  are  plotted  to  show  the  reduction 
of  lZQ.(«)l/n  with  increasing  no.  The  range  of  n2  has  been  continued 
down  to  zero,  by  taking  the  values  of  lZo(<»)|/n  for  the  non-rotating  round 
given  in  Ref.  2,  Pig.  9.  Per  values  of  n2  greater  than  O.I5,  1  Zq.(  CO  )l/n 
is  determined  by  equation  (33)  with  good  accuracy.  This  equation  shows 
that  the  dispersion  decreases  like  l/v5i^,  i.e,  inversely  as  the  square 
root  of  the  nozzle  offset  angle. 

The  absolute  value  of  the  dispersion  can  be  obtained  from  these 
graphs  if  IMq.1  is  known.  This  factor  depends  on  the  type  of  configuration 
considered,  and  in  particular  on  Ag  and  the  position  of  the  boosts  which 
affects  and  kp.  /m  analysis  _is  given  in  Ref.  2  of  the  various 
malalignments  contributing  to  ||aQ.j,and  typical  niamerical  values  of 
jpG-lkp^  (equal  to  lGp|/T  by  equation  (19))  are  given  in  Ref,  2,  para.  4..4 
for  tandem,  circumferential  and  wrap- round  boosts. 

4. 3  Dispersion  due  to  transverse  boost  force  '-  Results 

In  Pig.  11  ^m(“)|  is  plotted  against  h  -  for  the  particular  case 
n2  =  0.05  ft“^;  lZri(^”)  I  is  seen  to  increase  with  n.  The  same  figure 

shows  lZj(sq)|  where  Sj^  =  2,800  feet  for  a  non-rotating  round.  The 
dispersion  of  a  rotating  round  is  seen  to  be  less  than  that  of  the 
same  round  unrotated  by  a  factor  of  from  5  to  10  in  this  particular  case. 
In  fact,  for  a  non-rotating  round  (s)  ->■»  as  s  «> ,  and  for  values 
of  s  greater  than  2,800  feet  the  dispersion  is  increased  by  ^  log  s/2800 
approximately. 

A  rough  estimation  v/ill  nov/  be  made  of  the  largest  dispersion  likely 
to  arise  with  a  v/rap- round  boost.  The  tvro  main  causes  of  a  transverse 
boost  force  are: 

(i)  inequalities  in  the  inclinations  of  boost  thrusts  to  the  mean 
thrust  direction,  and 
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(ii)  inequalities  in  individual  thrusts. 

(i)  arises  from 

(a)  failure  of  the  individual  thrust  axes  to  coincide  vd.th  the 

nozzle  axis,  and 

(b)  inequalities  in 

The  order  of  magnitude  of  (i)  for  a  four  boost  assembly  corres¬ 
ponding  to  a  95?o  zone  is  about  0.001  radian. 

(ii)  is  proportional  to  the  radial  inclination  of  the  boosts  Ag. 
Taking  the  fractional  variation  in  thrust  of  a  four  boost  assembly  to 
have  a  3J/o  zone  of  I/40,  and  taking  Ag  =  20°,  the  contribution  from  (ii) 
turns  out  to  be -0.01  radian.  Thus  for  large  values  of  Ag,  say  10°  to 
20°,  (ii)  v^ill  predominate  and  if  we  write  as 


IittI  =  AgttQ, 

by  equations  (l8)  and  (20),  then  can  be  interpreted  as  the  fractional 
variation  in  thrust. 

V/ith  Ag  equal  to  20°,  is  about  0.01,  and  Pig.  11  gives 

|Zt(si)1  at  Sj^  =  2300  feet  to  be  about  2.0  for  the  non-rotating  round. 

Hence  the  dispersion  at  sp  =  2800  feet  is  0.02  radian  i.e.  just  over  1°. 

A  nozzle  offset  corresponding  to  n2  =  0,05  ft”^  v/ould  reduce  this 
dispersion  to  less  than  ■^°.  Larger  offsets  would  reduce  the  dispersion 
still  further. 

It  seems  that  dispersion  from  a'  transverse  boost  force  will  usually 
be  small  compared  with  dispersions  from  other  causes. 

5  •  Case  (b)  Constant  Spin 

5.1  Expressions  for  the  dispersion  functions 

5.11  Zqi(s)  and  Z(j(s) 

The  angular  acceleration  is  zero  in  this  case,  and  so  by  equations 
(24)  and  (25)  we  have  that  n2  =  0  and  nq  =  /^a  r^.  The  dispersion 

functions  Zq(s)  and  Zq.(s),  which  now  depend  only  on  the  parameters  n,  Sq 
and  nq,  are  evaluated  in  Part  II  para. -12.3  f or • sufficiently  large  values 
of  Tq.  The  solutions  obtained  take  on  different  forms  for  values  of  s 
less  than  and  greater  than  sp  =  (nq/2n)^.  (V/hen  nq  is  greater  than 
2v^  /x  N,  the  error  is  less  than  2/xN2  except  near  s  =3  so  where  it  is  less 
than  1/N) . 

Vdien  .  ■  V 


So  V®  ^  S3, 
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2  Zm(s)  = 


-D(aJ 


/  1-  _  _  D(^^s')  f  1  _  _1_ 

“1  Pi 


t 


•/s^  /  s 


-D(a^  )  D(-P3  ') 


1_ 


Pi'  -- 


1  i  sin  n(s-Sp,)  .  y  /  \ 
i _ i: _ 21.  +  1  n  ^(s) 


^^sc 


/s 


D(a,  ) 


L-  a-i 


Pi’  -J 


cos  n(s-s^) 
-  /s 


-  n  Si(s) 


11.12.31(30) 


2  e  °  Zr 


(s)  =  (I-  -  -i-) 

L  \Vs  '^^a'  Pi'  n''^s 


i[a-(s)-ob] 


'+  X 


■D(«s.)  D(H33')-i 


L  “l 


Pi’  J 


1  isinn(s-s  ) 

+  - ^ ^  -  in^(s) 

P 


l_  -fsp,  Vs 


+  X 


-D(a  )  d(-P3  On 

^o 


l_  «1 


Pi'  J 


cos  n(s-s  ) 


11.12.31(31) 

(36) 


and  when  so  s 


2  e”^°'o  z  (s)  = 


_  «-l  V^a  Vs/  \/s  Vs^ /_ 


^(«So)  D(-|3s^0- 


L  «1  ft.'  J 


“  1  ■ i  sin  n(s-s  )  .  ^  ,  V 

- -  +  xn^Js) 

y^a  ^fs  ^ 


P(“3.)  ^(-Pa  '  h 


L  Pi'  J 


"cos  n(s-SQ) 

_  Vs 


-  nSi(s) 


V  i[o'(sJ-c&] 
2D(o)e  P 


pi' 


1  2in  _ /j; 


;in(u-3p) 


11.12.34(45) 
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2  Zg{s)  .d0s')^1  .  1  ^ 


-“i  Vs  Pi'  Vs 


i{cr(s)-o-J 


+  V'D(«3o)  D(-Ps,'n 


Pi’  J 


—  VsT  /s 


+  1 


P(“Pso')' 

L  ai  ^  3i‘  J 


COS  n(s-s^) 
Vs 


+n 


^l(s) 


+  2i 


D(0) 


3  ' 
*^1 


1  2in  ,  , 

-  + -  D(6J 

6 


^in(i>-sp) 


Jsp 

11.12.34(46) 


where  Si(s)  and  ^2^^)  real  and  imaginaiy  parts  of  ^^(s)  defined 

by  equation  (28),  E>(u)  is  defined  by  equation  (29),  and 


6s  =j2r.s/'K 


«!  =  3i'  =  6i  =\j2r\/^ 


Vi^  =  =  -  ni/2n. 


All  arguments  of  D(u)  in  these  equations  are  positive,  and  the  bar  d«2notes 
the  complex  conjugate. 
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5*12  Ziji(oo)  and.  ZQ(eo) 


A  particularly  simple  solution  for  Z,j,(«)  and  Zq.(~)  is  obtainable 

when  =  ni/2n.  This  condition  will  almost  certainly  hold  for 

values  of  n^^  to  which  the  equations  of  para.  5.11  are  applicable.  In 
11.12.35  it  is  then  shovm  that,  on  taking  cTq  =  0, 


11.12.35(51) 

and  (52) 


11.12.35(53) 

and  (5Z*.), 

where  the  real  part  of  is  small  compared  with  the  imaginary  part. 

5.13  Approximate  formulae  for  Zj(s)  and  Zq.(s) 

If  in  addition  to  the  assumption  that  /s  «  7s^,  we  restrict  /s  to 
the  values  very  much  less  than  which  satisfy  ^ 

2n  (/s^  -  ^fs)^  »  1  , 


Z;j,(”)  =  iV2xn 


^Sq) 


n, 


Zg(”)  _  .1^  3(53^) 

n  n]_2  \i  n  n^ 


the  eqiiations  for  Zt(s)  and  Zq_(s)  given  in  para.  5.11  can  be  simplified. 
It  is  shown  in  11.12.32  that,  on  taking  cTq  =  0, 


R1  Zm(s)  = 


Si: 


incr(s) 


nq/s 


Imi  Zm(  s)  “ 


^1 


cos  n(s-So)- coserfs)  ^ 
- i ^ ^  -  n^.,(s} 

/s  - 


El 


2g(s) 


n 


^1 


Im 


Zg(s) 


.n 


ni  L 


^2(3)  - 


sinn  (s-Sq) 
n  v& 


(57) 


(38) 
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’ivhere  ^^(s)  and  ^2^^^ 


^l(s)  =  ^2%/n  [-A(6s^)  +  A( 6s)  cosn(s-SQ)  +3(6^)  sinn(s-SQ)] 
< 


11.12.32(2+0) 


^(s)  =  /2Vn  [b(6sq) 


+  A(6s)  sinn(s-so)  -  0030(3-3^)] 


11.12.32(2+1). 


In  11.12.33  it  is  shown  that 


Zr;(s) 


n 


4:  Im 


Zg(s) 


n 


The  above  equation  for  IinZQ^(s)/n  then  agrees  with  that  given  in  Ref. 1.6 

i.e. 


Zg(s) 

n 


G(6  Sq>6  s) 


(39) 


where 

S^>^s)  =B(6sq)  -  A3_(6sq) 


sin  -  (5„^-6  „ 


2(63)  cos  ^  (63^-63  2) 


11.12.32(43). 


If  we  let  s  tend  to  infinity  in  the  above  equations  for  Zj(s)  and 
Z(j(s),  it  is  found  that  their  limits,  Zr^{<”)  and  Zq,(«>)  ,  are  the  same  as 
those  of  para.  5.12.  The  reason  for  this  is  that  when  s  approaches 
^t(s)  and  Zq.(s)  arc  no  longer'  dependent  on  s  to  any  appreciable 
extent;  hence  the  expressions  for  Zj(s)  and  Zq,(s)  given  in  this 
paragraph  then  hold  for  all  values  of  s, 

5.2  Dispersion  due  to  boost  destabilising  couple  -  Results 

The  formula  used  for  evaluating  the  dispersion  due  to  a  boost 
destabilising  couple  is  equation  (391  above.  This  equation  should  be 
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ireasonably  accurate  for  valuea  of  n^  down  to  0,3  or  0,3  ft  depending 
on  n.  It  has  the  advantage  that  since  n^  does  not  occur  in  G( ^  > 

the  function  n^^l  Zq.(s)  1/n  can  be  plotted  in  terms  of  n  and  Sq  only  for 
various  values  of  the  range  s.  This  has  been  done  in  Pig,  12  for  n 
in  the  range  0.0025  to  0.01  ff^,  for  s^  =  0,  10,  20  and  30  ft  and  for 
1^0  "V^ues  of  s,  s  =  500  ft  and  s  =;<»  ,  It  appears  that  at  s  =  5OO  ft 
lZG.(s)j  has  already  converged  close  to  the  limit  j  Zq.(=o  )  |,  showing  that 
the  greater  part  of  the  dispersion  arises  in  the  first  few  hundred  feet. 

In  Pigs.  13  to  15  1Zg.(«)  l/n  is  plotted  against  nj  for  s  =  0,  10 
and  30  ft  and  the  three  values  of  n,  O.OO5,  O.OO75  and  0.01  ft“l.  The 
range  of  np  has  been  extended  dow’i  to  zero  by  using  the  values  of 

1  /n  for  a  non- rotating  round  given  in  Ref.  2,  Pig.  9.  It  is  seen 
■fchat  as  n.  ranges  from  0  to  1.5^  the  dispersion  of  a  non-spinning  round 
is  reduced  in  most  cases  by  a  factor  of  about  10.  Por  values  of  np 
above  1.5  the  reduction  will  be  proportional  to  ^/np  by  equation  (39), 
i.e.  the  dispersion  varies  inversely  with  the  rate  of  rotation  of  a  round. 

The  function  Z(p(G)/n  is  shovm  accurately  evaluated  in  Pig,  I6  fpr 
the  particular  case  v/hen  n  =  0.0075  ft"l,  s^  =  10  ft  and  np  =  1.5  ft“2. 

The  values  of  ZQ.(s)/n  determined  by  the  approximate  formula,  equation  (56), 
for  this  case  oscillato  from  side  to  side  of  the  correct  values  and  do 
not  differ  from  them  by  more  than  0,2  ft.  The  values  determined  by 
the  more  approximate  formula,  equations  (37)  ^-^nd  (38),  lie  on  the 
straight  line  R1  ZQ.(s)/n  =  2/(l.5)^  =  0,9  ft,  with  Im  ZQ.(s)/n  varying 
from  0  at  launch  to  6«1  ft  at  infinity. 

Pig.  17  shov/s  the  yav/  corresponding  to  Pig,  I6.  The  convergence 
of  the  yaw  to  zero  as  s  is  slow  compared  Y/ith  the  rate  of  convergence 
of  the  dispersion. 

6  Choice  of  Spin  during  Boost 

6.1  It  has  been  seen  in  paragraphs  4.  and  5  that  malalignment  dispersions 
can  be  considerably  reduced  by  sufficiently  large  spins.  The  question 
that  will  now  be  considered  is  -  ^By  v/hat  factor  is  it  worthwhile 
decreasing  the  dispersions  of  the  non- spinning  round?’  To  answer  this 
it  is  necessary  to  compare  the  malalignment  dispersion  for  a  non-rotating 
round  v/ith  its  other  dispersions  which  are  not  affected  by  the  initial 
spin- form.  This  v/ill  now  be  done  for  a  Scaslug-typo  missile  with 
a  vrrap-round  boost.  The  results  given  in  Refs.  2  and  3  state  that  for 
an  unspun  round  the  various  dispersions  likely  to  occiir  are: 

Malalignment  dispersion 


dispersien 


10°  venturi  offset  angle 

3^0  -  6° 

depending  on  position  of  boosts 

20°  venturi  offset  angle 

hP 

for  most  boost-positions 

Other  dispersions 

Wind  (lO’/scc  error) 

2° 

Separation  dispersion 

1°  -  2° 

(possible  maximum) 

Tracking  Beam  displacement 

.2° 

due  to  avoiding  action  of 
enemy 

2° 

due  to  variations  in  time 
of  burning . 

The  resultant  dispersion  is  found  to  be  from  5^  to  7^7^* 
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For  a  rotating  round  the  following  table  gives  the  resultant 
dispersion  when  the  iDalalignment  dispersion  is  reduced  by  a  factor  of 
2^3)  or  If* 

Factor  by  which  n^alalignment  Resultant  Dispersion 

dispersion  is  reduced 


1 

5°  - 

7i° 

2 

4°  - 

50 

3 

I 

q 

4i° 

4 

7I0  _ 

^2 

4i° 

It  is  seen  that  there  is  no  great  advantage  in  decreasing  the  mlaJLignment 
dispersion  by  more  than  a  factor  of  2  to  3*  This  conclusion  naturally 
depends  on  the  values  for  the  various  dispersions  taken  above:  if  for 
instance  the  deviations  due  to  wind  and  tracking  beam  displacement  have 
been  overestimated,  then  a  greater  deduction  in  the  malalignment  dispersion 
would  give  a  greater  proportional  drop  in  the  resultant  dispersion*  Hov/- 
ever  this  drop  does  not  appear  appreciable  unless  the  errors  have  been 
grossly  overestimated,  and  this  is  not  considered  likely*  For  example,  if 
the  wind  and  tracking  errors  were  halved  and  the  separation  dispersion 
taken  to  be  zero,  the  resultant  dispersion  would  be  - 


Factor  by  which  malalignment  Resultant  Dispersion 

dispersion  is  reduced 


1 

4°  - 

2 

olo  ^ 

3i° 

3 

2°  - 

O^o 

4 

2”  - 

2^0 

From  these  figures  it  is  seen  to  be  hardly  worthv/hile  aiming  at  a  factor 
of  reduction  greater  than  3* 

6*2  By  way  of  example  it  is  now  shoTO  how  the  ballistic  dispersion  of  an 
unspun  round  can  be  reduced  by  a  factor  of  3  by  employing  methods  (a)  and 
(b) ,  para*  2.1,  of  imparting  spin, 

6,21  Method  (a)  Offset  nozzles 

Fig*  18  gives  the  reduction  R  in  dispersion,  achieved  by  offsetting 
the  boost  nozzles  of  a  round,  in  terms  of  the  paraiBetcr  n2  which  is  related 
to  the  offset  angle  hy  equation  (35)  •  For  values  of  s^  other  than  10  ft 
the  corresponding  graphs  can  be  obtained  from  Figs*  6  to  8  and  will  not  be 
greatly  different  from  Fig.  18* 

For  R  =  3  and  n  corresponding  to  the  above  Scaslug  figu;“es,^say 
0,005  to  0*007  ft“^,  we  find  that  n2  lies  botv/oen  0*06  and  0,075  ft 
The  corresponding  nozzle  offset  angle  is. 2*0^  to  2*5^^  and  the  corres¬ 
ponding  angular  acceleration  at  launch  is  24-  to  30  rad/ sec  f  taking  *tho 
acceleration  to  be  4-00  ft/sec2* 

The  maximum  value  attained  by  the  spin  during  boost  can  be  found 
by  equations  (lO)  and  (l4*) ,  which  give 


29. 


CONFIDENTIAL 


CONFIDENTIAL 


Technical  Note  No.  G-.V/’,  177 


where  a  is  the  acceleration,  Oq  and  t  arc  defined  by  equations  (13) 
and  (15),  and  f(0o)  function  plotted  in  Pig.  22.  Taking 

Qq  ='’':/6,  i.c.  djy^dq  =  0.50,  I  -  0.0001-0  ft“l  (a  value  appropriate  to  the 
302/STV) ,  we  find  that 


=  25  to  31  rads/sec, 

i.e. 


=  h  to  5  revs/scc, 

for  an  offset  angle  of  2.0°  to  2,3°- 

6'. 22  .  Method  (b)  -  launching;  Yilth  spin  from  a  spiral  launcher 

Fig.  19  chov/s  R  against  v/hich  is  proportional  to  the  constant 
spin  Tq,  in  the  case  when  Sq  =  10  ft,  and  n  =  0.005,  0,0075  and  0,01  ff*^. 

Wien  R  =  5  and  n  ranges  f^oin  0.005  to  0,007  ft""^,  we  find  that 
n-L  lies  betv/ecn  0,50  and  0,58  ft'“2.  The  value  of  the  spin  corresponding 
to  an  acceleration  of  4*00  ft/sde^  then  turns  out  to  be 


Tq  k*2  to  5 A  rads/ see, 


i.e  is  soincv/hat  less  than  1  rcv/scc. 


The  angular  acceleration  of  the  round  when  moving  up  the  laiaicher 
is  greatest  when  the  acceleration  is  greatest ,  i.e.  at  laiinch,  and  is 
equal  to  an2^/27s^  by  equations  (l)  and  (25).  For  Gq  =  10  ft  and  the 
above  values  of  a  and  nq  ,  the  maximum  angular  acceleration  of  the  round 
is  from  19  to  224-aradG/sec'^. 


6,25  To  summarise  we  can  say  that  v/ith  both  methods  (a)  and  (b)  the 
round  must  be  subjected  to  angular  accelerations  of  the  same  order,  about 
5  to  revs/sec^,  but  that  the  maximuia  rate  of  rotation  for  method  (b)  is 
only  1  rev/scc  compared  v/ith  k  revs/scc  for  method  (a). 


7  Dispersion  follov/ing  separation  -  Values  of  parexmeters 

•The  chief  parameters  affecting  the  dispersion  of  the  dart  after 
separation  are  its  lift  and  stability  properties  and  its  spin-form.  In 
this  paragraph  we  consider  in  turn  the  values  that  these  parameters  arc 
likely  to  assume  in  pi'^acticc, 

7.1  Lift 


If  NmQ-g  be  the  maximum  sca-lcvel  lift  that  the  dart  can  develop  at 
velocity  V  and  v^ing  incidence  5  ,  then  vro  can  v/rito 


Lift  =  kLV^S  =  Nn^g, 


(40) 


where  kjj  is  supposed- constant.  '• 
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If  ni^  is  the  mass  of  the  dart  at  separation  and  if  we  define  t  by 


^  =  ^iJ^o  > 


we  cam  use  equation  (1+j0)  to  express  ■&  in  terns  of  N,  thus 

I  =  Ng/V%  . 


TaJcing  N  to  range  from  25  to  100  vfhen  V  =  2500  ft/sec  and  6  =  20°, 
we  find  that  Z  lies  between  0,0001-  and  O.OOI6  ft“^. 

7.2  Stability 

The  stability  constant  of  the  dart  is  defined  as 


n  =  iTCd/kp, 


where  d  =  distance  of  C.P.  aft  of  C.G. 

.kp  =  radius  of  gyration  in  pitch  of  dart  at  separation. 

This  definition  coincides  with  that  given  for  the  complete  weapon  dtiring 
boost  in  para.  3.2. 

If  we  put 


d  =  eL 


kp  =  L/2V5 


appro  xdma  t  oly , 


vfhere  L  is  the  length  of  the  projectile,  we  have 


r?  =  12ee/L, 


(U) 


For  values  of  L  between  I6  and  25  feet  and  of  Z  between  0.0004  and 
0.0016  ft'l,  we  find  that  n  lies  in  the  range  O.OO4  to  0,01  ft  when 
e  =  1/12. 

7,3  Spin~f orm 

During  the  fall-off  of  thrust  at  the  end  of  the  boost  period  and 
during  the  ■  consequent  separation,  the  spin  will  rapidly  decrease  towards 
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zero  on  account  of  the  large  damping  affect  of  the  v/ings.  In  practice 
a  small  spin  usually  remains  due  to  aileron  or  wing  malalignment  incidence 
this  spin  may  be  in  the  opposite  direction  to  the  original  spin. 

It  was  foxind  in  a  number  of  S.T.V.l  firings  that  the  spin  of  the 
dart  remained  fairly  constant  and  rarely  exceeded  1  rev/scc, 

A  solution  for  the  ^in/velocity  ratio  is  obtainable  when  the 
retardation  is  assumed  proportional  to  (velocity)2.  Integration  of 
equation  111.4.12(4)  yields 


r 

V 


r,  -YR(s-so)/cn 

1  -  c 

'4 

Ldt  \V/4 

_C _ 

YrVo  ' 


(h2) 


v/here  the  siiffix  o  nov/  denotes  the  value  of  a  quantity  at  separation, 
and  where 


C  =  M  of  I  of  dart  in  roll 


=  dai]5)ing  moment  in  spin/Vr  (supposed  eonstant)* 

Equation  (^2)  shows  that  as  s-s^  inereases  above  about  r/V  tends 

rapidly  to  the  constant  value 


which  is  zero  when  no  malalignment  incidences  arc  present. 

In  view  of  the  above  remarks  the  ratio  spin/velocity  will  be 
supposed  constant  and  denoted  by  Y. 

8  Estimation  of  dispersions  arising  after  separation 
8*1  Definitions,  etc. 

The  dispersion  of  a  roulid  at  any  point  after  separation  is  defined 
as  the  angle  between  the  direction  of  motion  of  the  C.G.  of  the  round  and 
its  direction  at  separation.  The  components  of  the  dispersion  on  two 
perpendicular  planes  through  this  direction  are  taken  as  the  real  and 
imaginary  parts  of  the  complex  dispersion  Z;  jzj  is  then  the  angular 
dispersion.  *  *  ♦  '  • 

The  dispersion  arising  from  an  initial  yaw  at  separation  is 

denoted  by  and  lies  in  the  plane  containing  the  initial  direction  of 
motion  and  the  missile  axis.  ,  The  dispersion  arising  from  an  ini.tial 
rate  of  turn  of  missile  oxis  Sq  ^t  separation  is  denoted  by  Z2  and  lies 
in  the  plane  containing  the  initial  instantaneous  oscillations  of  the 
axis. 
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The  aerodynamic  lift  and  moment  malalignment  angles 

defined  as  the  angles  which  the  air  flow  makes  with  the  axis 
(at  infinity)  when  the  aerodynamic  lift  and  moment  are  zero, 
are  used  to  denote  the  complex  dispersions  that  arise. 

In  Part  III  the  coasting  equations  of  motion  are  solved  to  give 
Z2,  Zjj  and  when  the  ratio  spin/velocity  is  constant.  We  now  make 
use  of  the  results  of  the  simplified  solution  of  III. 5.  ^or  sufficiently 
large  values  of  the  range  s-s^,  it  is  found  that  certain  transient  terms 
become  negligibly  small  leaving  particularly  simple  expressions  for  the 
dispersions.  The  condition  is  that  s-s^  should  be  greater  than  about  2/6, 
The  error  in  the  solution  can  then  be  about  Q-t  s-s^  =  2/6,  1'^%  at 
s-Sq  =  L/6  and  at  &-Sq  =  G/z  etc. 

8.2  Dispersion  caused  by  unclean  separation 

7/e  now  consider  the  possibility  of  asymmetrical  detachment  of  the 
boosts  producing  a  dispersion.  The  dispersion  arising  from  am  initieil 
yaw  is  shown  in  III. 5.4  to  be  entirely  transient.  The  dispersion  due  to 
am  initial  rate  of  turn  of  the  axis,  on  the  other  hamd,  becomes 


amd  “jjI  are 
of  the  round 
Zjj  and  Zm 


Z 


2 


(43) 


where  is  the  velocity  at  separation,  hence 


IZ2I  =  L/ia:  by  equation  (4I). 


Taking  =,  1,500  ft/  sec,  e  "  1/12  and  L  =  16  to  25  feet,  this  gives 


IZ2I  =  0.6  to  1.0  IsJ  degrees 


(when  is  expressed  in  rads/sec). 

An  alternative  expression  for  this  dispersion  can  be  obtednod  by 
expressing  in  terms  of  the  maximum  yaw  that  it  would  give  rise  to  in  the 
ensuing  motion.  Prom  III.4. (l3),  (l5)  and  (lO) ,  and  III.5.l(l),  the 
solution  for  the  yaw  is 


Vq)  e  sin  p(s-Sq), 


where 


-  {«/2)^ 


111.5.1(2). 
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This  equation  shows  that  the  maximum  yaw  is 

M 


M 


=  (^P  e 


-iteAp 


t  t  '/n 
0  0 


for  vaiues  of  -C  and  n  in  the  ranges  considered  in  paragraphs  7.1  and  7.2. 
Equation  (4-3)  nov/  gives 


IZ2U  1  SmI  ==  1  y W12  e  by  equation  (i+l). 


Taking  s  =  1/12,  •£  =  0.0004-  to  0.0016  ft""^,  and  L  =  16  to  25  feet,  we  find 
that 


Ugl  =  0*08  0.2 


This  shows  that  for  a  dispersion  of  2°  to  arise  the  first  oscillation 
in  yaw  after  separation  would  have  to  hive  an  amplitude  of  10°  to  25° 
corresponding  to  the  range  of  values  of  the  parameters  taken  here. 

8. 3  Dispersion  due  to  lift  maialignment 

The  maximuim  value  that  the  angular  dispersion  |Zj^|  ever  attains 
is  shown  in  III. 5.2  to  be 


ln2  ..  Y^l 
J  (n^  -  Y^)  ^  +  Y^-^— 


where  ^  and  •  n  are  the  lift  and  stability  constants  already  discussed, 
and  Y  is  the  ratio  spin/vclocity. 

As  the  range  increases  after  separation  the  locus  of  in  the 
complex  plane  tends  to  a  circle  whose  radius  is  the  second  term  on  the 
right  hand  side  of  this  equo.tion  and  whose  centre  is  at  a  distance  from 
the  origin  given  by  the  first  term. 

Pig.  20  shows  l%j  L  plotted  against  [yI/'S  for  various  values 

n/C,  All  curves  lie  below  the  curve 


max 
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IzlI  Al  =  se/W  1, 

max 


which  corresponds  to  an  infinitely  stable  round,  i.e,  n  =  «>. 

It  is  seen  from  Fig,  20  that  IZr 1  is  less  than  for  all  values 

'  -^'max  ^ 

of  n  provided  Y  >  2.6,  (i.e.  provided  the  spin  is  greater  than  1,2  to  4.8 
rad/sec,  taking  the  velocity  to  be  1500'/eec.  and  •&  to  lie  between 
0.0004  and  0,0016  f  t“^) .  For  smaller  spins  is  less  than  oCj^  only 

over  a  limited  range  of  values  of  s-s^.  Thus  for  a  non- spinning  round, 
it  can  be  shown  that 


!^l!  "  A”Sq)]  j 


and  from  this  it  follows  that  <a-^  provided 


s-Sq  <  (l  -  ■6Vn^)A 

4=  lA 

for  values  of  ■C  and  n  arising  in  practice,  i.e.  provided  s-Sq  is  loss 
than  625  to  2500  ft,  corresponding  to  -6  from  O.OOI6  to  O.OOO4  ft“^. 

8,4  Dispersion  due  to  moment  malalignment 

The;  locus  of  Zj^  in  the  complex  plane  is  shown  in  Part  III,  para.  5.3 
to  tend  to  a  circle  as  the  range  increases  after  separation.  The  maximum 
Value  attained  by  jZj^l  is  found  to  be 


'/(n^  -  — 


Curves  of  jZjjj^^/aj^  against  IyIA  are  shown  in  Fig.  21  for  various  values 
of  r/C.  The  equation  of  the  envelope  is 

K!  /“m=  [^  +  +  y2]  /  1y1 

max 


4=  1 

for  sufficiently  large  values  of  Ivj  /•&.  It  is  easily  seen  that  for 
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smaller  spins  less  than  ciy^  over  a  limited  range  of  s-s^. 

Thus  for  a  non-rota tir»g  round 


and  so  it  is  required  that 


s-Sq  <  (l  +  ■&  Vn^)/ 1 


-  i/e. 

The  resiilts  of  para.  8,3  and  8.4  can  be  summarised  by  saying  that 
provided  the  spin  is  not  too  small  the  dispersion  is  less  than  the 
malalignment  angle.  For  small  spins  this  is  only  true  over  a  limited 
flight-range,  which  is  about  l/C  if  there  is  no  spin  present  at  all. 

9  Conclusions 

9*11  If  the  aerodynamic  spin  damping  moment*  is  neglected  the  spin  of  a 
round  rotated  by  offset  boost  nozzles  increases  linearly  with  time,  the 
angular  acceleration  being  proportional  to  the  couple  produced  by  the 
nozzles*  With  guided  weapons  conditions  are  usually  such  that  the 
aerodynamic  damping  moment  exceeds  the  boost  couple  tov/ards  the  end 
of  the  boost  period,  causing  the  spin  to  decrease.  Good  agreement 
has-  been  found  between  practical  and  theoretical  spin-forms. 

9.12  When  the  launching ‘ spin  is  zero  the  spin  of  a  given  round  at  each 

instant  of  the  boost  period  can  be  shovm  to  be  proportional  to  the  boost 
couple.  ,, 

9.13  ^^or  the.  purpose  of  evaluating  the  angular  deviation  a  knowledge  of 
the  spin-forrri' over  the  first  part  of.  the  flight  only  is  necessary,  as 

it  is  while  the  velocity  is  lov/  that  the  greatest  part  of  the  deviation, 
arises.  Trajectories  from  502/STV  firings  confirm  this  fact;  the 
angular  deviation  scarcely  changes  after  the  first  500  ft  of  the  3250  ft 
of  boost  range.  The  spin  dan5)ing  properties  of  a  missile  are  therefore 
of  secondary  importance  and  can  be  justifiably  neglected  in  any  evalimtion 
of  the  angular  deviation. 

9.21  Of  the  various  dispersions  arising  during  boost  that  due  to  a  boost 
destabilising  couple  has  been  most  fully  treated,  as  this  is  of  overriding 
importance  for  a  guided  v/eapon.  The  decrease  of  dispersion  when  spin  is 
imparted  by  offset  boost  nozzles  is  sliown  in  Pigs.  5  to  8,  and  when  imparted 
by  the  launcher  is  shovm  in  Pigs.  15  to  17.  For  large  offset  angles  the 
dispersion  decreases  roughly  as  the  square  root  of  the  offset  angle,  and 
for  large  launching  spins  the  dispersion  is  inversely  proportional  to 

the  spin, 

9.22  The  faotor  by  which  it  is  worthwidle  reducing  the  ballistic  dispersion 
is  limited  by  the  existenoo  of  dispersions  that  are  not  reduced  by  spin. 

For  a  Seaslug-like  missile  thero  seems  to  be  no  advantage  in  a  reduction 

of  more  than  3-  According  to  the  theory  reductions  of  this  magnitude 
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can  be  produce:!  in  general  by  ofisetting  the  boost  nozzles  by  less  than  5^. 
Y/ith  a  spiral  launcher  a  reduction  by  3  coulu  be  attained  with  a  launching 
5p>in  of  just  less  than  1  rev/sec. 

9.23  It  shouli  be  noted  that  the  malaligrunent  iispersion  cannot  be  reduced 
by  spin  below  a  certain  value  on  account  of  malaligninent  of  the  principal 
longitudinal  axis  of  inertia^  the  effect  of  wiiicii  increases  with  spin. 

An  estimate  made  at  the  outset  (para.l.f)  showed  that  in  a  typical  case 
the  maximum  reduction  would  occur  at  about  4  revs/scc.  As  this  value  is 
much  higher  than  that  required  for  the  worthwhile  reduction  of  3  mentioned 
above,  the  effect  of  inertia  axis  malaligninent  is  of  no  importance. 

9.24  It  is  appreciated  that  considerable  engineering  difficulties  would  be 
associated  with  the  use  of  a  spiral  launcher.  A  great  advantage  of  such 

a  launcher  is  however  that  the  maximum  spin,  i.e.  the  launching  spin,  would 
be  much  less  than  the  maximum  spin  produced  b^^  offset  nozzles  for  the  same 
reduction  in  the  malalignment  dispersion.  For  Seaslug  the  maximum  spin 
would  be  about  0.8  rcv/sec  compared  with  4  revs/seo.  If  the  maximum  spin 
produced  by  offset  nozzles  v/ere  unacceptably  lai'ge,  it  v/ould  then  be 
necessai-y  to  use  a  spiral  launcher  to  achieve  the  desired  reduction  in 
dispersion. 

9.25  It  shouli  be  pointed  out  that  the  value  5  for  the  reduction  factor 
mentioned  abovo  depends  on  the  relative  ir.iportance  of  malalignment 
dispersions  and  dispersions  due  to  v/ind  error  and  other  errors  not  affected 
by  spin,  and  should  be  revised  in  the  light  of  future  information  about 
these  quantities.  Should  larger  reductions  by  a  factor  of  5  or  more  be 
indicated,  then  a  spiral  launcher  would  be  the  only  means  of  achieving 
them. 

9.31  Y/ith  regard  to  dispersion  of  the  dart  due  to  as3nnmetrical  detachment 
of  the  boosts  at  separation,  it  can  be  shown  that  dispersion  will  be  sm^ill 
so  long  as  the  dart  is  not  set  v/ildly  oscillating  v/ith  large  angles  of  yav/ 
of  more  than  about  10^.  In  view  of  the  successful  operation  of  the 
separating  device  in  test  vehicles  it  is  not  likely  that  such  a  large 
disturbance  to  th^  motion  v;ill  in  general  ecci-ir,  and  so  it  is  concluded 
that  dispersion  from  this  souoce  will  be  negligibly  small. 

9.32  Dispersions  arising  from  the  aerodynamic  asymmetry  of  the  dart 
naturally  depend  rather  critically  on  the  3x)in.  By  assuming  the  spin 
tc  be  constant,  simple  expressions  have  been  obtained  for  the  maximum 
dispersions  due  to  lift  and  moment  malalignrajnt  angles  and 

It  is  found  that  so  long  as  even  a  small  spin  is  present  the  malalignment 
dispersions  will  not  exceed  and  respectively.  This  means  that  if 

the  malalignment  angles  can  be  kept  dov/n  to  the  order  of  0.1^  the  resulting 
dispersions  should  be  negligibly  small. 

9.33  the  spin  is  zero  a  limit  must  be  imposed  on  the  range  of 
uncontrolled  flight  if  large  dispersions  are  to  be  avoided.  From  the 
condition  to  be  satisfied  in  this  case  it  is  found  that  when  aj^  and 

equal  0.1^,  the  dispersion  is  less  than  1^  provided  the  uncontrolled 
range  is  less  than  6,250  ft  to  25,000  ft,  according  to  the  lift  properties 
cf  the  round.  If  the  velocity  is  in  the  region  of  1,500  ft/sec  during 
this  period,  the  time  of  uncontrolled  flight  will  in  no  case  be  required 
to  be  less  than  about  4  secs.  In  view^  of  the  practicability  of  roll- 
stabilising  a  round  and  bringing  it  under  control  w^ithin  1  or  2  secs 
after  separation,  it  can  be  concluded  that  no  large  dispersion  from  this 
cause  is  to  be  expected. 
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Pi^HT  II 

THE  EQUATIONS  OF  MOTION  iJO  THEIR  SOLUTION  DURING  BURNING 


1  Introduction 


During  the  last  15  years  a  considerable  amount  of  work  has  been  done 
in  this  country  on  the  subject  of  rocket  motion  with  a  view  to  applying  it, 
in  the  most  part,  to  small  rockets  of  the  fin  -  and  spin- stabilised  type. 

The  mathematical  part  of  this  work  has  been  published  in  Ref  1:  the  general 
equations  of  motion  of  a  body  losing  irass  are  formulated  and  the  equations 
of  motion  of  a  rocket  in  flight  are  deduced  and  solved,  the  assumptions 
being  clearly  stated  at  each  stage.  The  objeot  of  Part  II  of  this  note 
is  to  adapt  and  extend  this  work  to  apply  to  a  guided  weapon  during 
\ancontrolled  boost. 

The  equations  of  motion  formulated  in  Ref,  1  include  all  conceivable 
effects  to  which  a  guided  weapon  might  bo  subjectol  during  boost,  apart 
from  aerodynamic  asymmetry.  This  v/’ill  nov/  be  introduced,  and  other 
effects  v/hich  are  small  for  contemporary  guided  weapons  v^’ill  be  omitted: 
these  include  l^Iagnus  effects  and  those  malalignment  effects  that  become 
appreciable  only  at  hi^  spins,  such  as  the  principal  axis  of  inertia 
malalignment  effect  (see  1,1, 4)  which  is  possibly  the  most  important  of 
them* 


In  Ref.  1  the  reduction  of  the  equations  of  motion  to  a  form  suitable 
for  solution  proceeds  under  the  assumption  that  the  drag  and  gravity  force 
are  negligible  compared  vd.th  the  thrust.  iil though  this  holds  in  practice 

for  guided  v;‘eapons  diiring  boost,  it  is  not  essential  for  the  manipulation 
of  the  equations  of  motion  along  thi  linos  of  Ref.  1  (provided  the  direction 
of  launch  is  not  near  the  vertical)  ,  The  equations  so  obtained  then  hold 
for  the  motion  after  burning  as  well,  and  no  reformulation  of  the  equations 
of  motion  is  required  for  Part  III. 

A  solution  of  the  equations  oC  motion  in  terms  of  integrals  is 
obtainable  under  any  set  of  assumptions  for  which  the  wavelength  of  yaw 
is  constant:  in  Ref.  1  it  has  been  assumed  that  the  ratio  spin/velocity 
is  constant  -  an  assumption  that  is  not  generally  valid  for  a  guided 
weapon.  If  instead  it  is  assumed  that  the  spin  is  so  small  that  the 
processional  effects  are  negligible,  a  solution  is  obtainable  that  is 
applicable  to  a  slowly  rotating  projectile.  The  method  of  solution  here 
will  therefore  follow  the  lines  of  Ref.  1  with  the  precossional  terms 
omitted,  but  with  the  spin-form  as  an  arbitrary  function. 

The  evaluation  of  these  integrals  for  the  ballistic  deviations  forms 
the  remainder  of  Part  II  of  this  note.  It  is  assumed  that 

(i)  the  linear  acceleration  is  constant, 

(ii)  the  angular  acceleration  is  constant,  and 

(iii)  the  component  of  thrust  perpendicular  to  the  trajectory  is 

large  compared  with  the  lift. 

On  account  of  'the  larger  aerodynamic  surfaces  of  a  guided  weapon,  (ii) 
and  (iii)  are  vnlid  only  at  the  beginning  of  the  boost  period  in  general. 

Any  angular  deviation  of.  the  roiind  will  however  quickly  tend  to  a  limit 
on  account  of  the  increasing  spin  and  velocity,  and  can  still  bo  aocurateljr 
evaluated  under  those  assiamptions. 
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The  integrals  for  the  yaw,  the  axis  inclination  and  the  angular 
deviation  are  evaluated  in  terms  of  the  functions  of  App.A.  These 
functions  have  already  been  used  in  Ref.l  with  positive  arguments,  and 
they  are  now  extended  to  negative  and  imaginary  arguments  in  the 
appendix.  The  evaluation  of  the  yaw  integrals  can.be  carried  out 
exactly.  The  double  integrals  occurring  in  the  angular  deviation  cannot 
be  evaluated  in  general  for  all  values  of  the  angular  acceleration.  Two 
approximate  methods  are  developed  which  cover  all  cases  except  small  or 
zero  launching  spin  together  with  small  or  zero  angular  acceleration. 

The  solution  for  a  non-rotating  round  is  of  course  already  known^*^., 
and  by  interpolation  all  practical  cases  can  be  evaluated  with  good 
accuracy. 

2  Axes  and  notation 

2.1  The  motion  is  referred  to  axes  OX,  Y,  Z  v/hose  directions  are 
fixed  in  space  and  which  are  defined  as  follows: 

0  is  the  C.G.  of  the  missile 

.OZ  lies  in  the  direction  of  the  axis  of  the  launcher,  making 
a  Q.E.  of  a  with  the  horizontal  plane 

OX  lies  in  the  vertical  plane  through  OZ  and  is  perpendicular 
to.  OZ^  in  a  dovvTLwards  direction,  . 

OY  completes  the  right-handed  set  of  mutually  perpendicular  . 
axes  and  thus  lies  horizontally  to  the  left  when  looking 
along  OZ.  •  ^ 

The  ppints  X,  Y,  Z  are  taken  to  lie  on  the  unit  sphere  centre  0. 

2.2  OP  vector  meeting  the  unit  sphere  at  P,  and  'if  Pj^  is 

the  projection  of  '  P  on  to  the  plane  OXY,  then  OP  is  determined  by  z, 
the  complex  number  v^hose  real  and  imaginary  parts  are  the  coordinates 
of  referred  to  axes  OX,  OY.  In  fact,  P  is  located  by  the 
relations  =  XOPjvy  =  arg  z  and  &  (=  ZOP)  =  sin"^  OPj^  =  sin“^  jz]. 

The  following  axes  and  directions  are  defined  by  complex  numbers  in 
this  manner, 

2.21  OT  is  the  direction  of  the  velocity  V  ;  it  is  the  tangent 

— r\ 

to  the  trajectory  of  the  projectile.  The  angular  deviation  is  ZT  and 
is  determined  by  the  complex  number  Z, 

2.22  OA  is  the  missile  axis,  defined  in  3.12;  its  *  inclination' 
ZA  is  determined  by  the  complex  number 

2.23  The  yaw  TA  is  determined  "by  the  complex  number  ',~=3  •  By 
solving  the  spherical  triangle  TAZ,  TEl  can  of  course  be  related  to  Z 
and  ^  . 

The  solution  of  the  motion  in  the  plane  OXI  is  considerably 
simpler  v/hen  the  follov/ing  assLunption  is  made;  a  linean  theory  is 
then  obtained, 

A,1  It  is  assumed  that,  the  sines  and  cosines  of  jzj,  |Sj,  1^1, 
and  all  the  malalignment  angles  are  replaceable  by  jzj,  1, 

and  uniiy  respectively  to  sufficient  accuracy,  and  that  the 
derivatives  of  ' these  quantities  v/ith  respect  to  time  are  of  the  same 
order  of  magnitxade. 
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The  relation  between  Z,  ^  and  is  then 

Z  =  S  -  S  (1) 

2.3  AOA*  is  a  reference  plane  fixed  arbitrarily  in  the  missile, 

2.31  The  total  angle  cr  turned  through  by  the  missile  about  its 
axis  is  the  angle  between  the  planes  AOA’  and  ZOX. 

'i'h  ^ 

2.32  OAp  is  the  direction  of  the  axis  of  the  r  asymmetry,  AA^ 

-  fs 

is  called  the  malalignment  angle  and  is  denoted  by  A’AA^,  denoted 

by  0^,  determines  the  orientation  of  the  asymmetry  in  the  missile. 

Notation  In  the  follov/ing  paragraphs  a  letter  is  underlined  to  denote 
a  vector.  The  same  letter  not  underlined  denotes  the  complex  number 
determining  this  vector  in  the  manner  described  in  2.2. 

3  Aerodynamic  forces  and  couples 

3.11  A, 2  It  is  assumed  that  the  shell  of  a  perfectly  manufactured 
missile  possesses  geometrical  symmetry  of  order  greater  than  2  (see 

Ref, 1.5)  about  a  longitudinal  axis  OA. 

Yfhen  the  missile  lies  in  an  air  stream  flovding  in  the  direction  ^ 
at  a  great  distance  av/ay,  it  follows  from  A.  2  that  the  aerodynamic  force 
R  and  couple  _!  v/ill  both  lie  in  the  direction  OA. 

3.12  An  imperfectly  manuf actui'ed  missile  on  the  other  hand  has 
small  asymmetries  in  its  geometrical  shape  which  prevent  the  axis  from 
being  defined  as  the  axis  of  symmetry.  The  definition  given  here  depends 
on  the  implicit  assumption  that  both  perfect  and  imperfect  missiles  can 
be  launched  from  the  same  laiuncher.  This  means  that  there  is  a  set  of 
points  (of  more  than  2)  on  any  given  missile  (namely  the  points  of  contact 
betv/een  missile  and  launcher)  which  can  be  brought  into  spatial  coincidence 
with  the  same  set  of  points  on  a  perfect  missile.  The  axis  of  the  given 
missile  is  then  defined  as  the  direction  through  its  C.G-,  parallel  to  tlie 
axis  of  the  perfect  missile, 

3.13  Two  directions  OAj,  OA^^  can  be  determined  in  the  body  of  the 

missile  such  that,  v/hen  the  airflow  is  in  the  direction  A^O  at  a  great 
distance  away,  the  resultant  force  is  parallel  to  OA^,  and  when  in  the 
direction  Aj^O  the  resultant  couple  is  parallel  to  For  a. perfect 

'  '  '  ”  r%  n 

missile  both  these  directions  coincide  with  OA.  The  angles  AAj^ 
are  the  lift  and  moment  malalignment  angles  and  are  denoted  by 
respectively, 

3.2  The  components  of  R  and  £  along  the  missile  axis  R^^,  and 
perpendicular  to  the  missile  axis  Rp,  fp  sire  analysed  as  follows: 

3.21  R^  is  termed  the  air  resistance  and  is  oq^al  in  magnitude 
to  the  axial  drag 

3.22  can  be  divided  into  tv/o  parts,  due  to  fins  (or  ^vings) 

having  an  offset,  and  the  restoring  couple  that  arises  v/hen  the 
missile  is  spinning.  If  the  magnitudes  of  are  taken  to  be 

YrV^,  ‘“YpVr  along  0^  then 


r;,  =  “(YRVr  +  YpV2)  oA  , 


(1) 


v/here  r  is  the  spin,  i.e. 


r 


do- 

dt 
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3.23  Rp  arises  from  the  yav/j  it  is  not  quite  the  same  as  the 

lift  which  is  taken  perpendicular  to  the  airflow.  To  the  accuracy  of 
A.l  the  component  of  Rp  in  the  plane  OXT  has  magnitude  [Rpj^  and  so 

when  referred  to  axes  OX,  OY  is  determined  hy  the  complex  number 


(2) 


If  L  is  the  component  of  the  lift  force  in  the  plane  OXY,  referred  to 
axes  OX,  OY,  then 


L  =  Rp-1%l^ 

=  :z:  , 

(4) 

v/here 

! 

11 

(5) 

i^T 

A  missile  with  lift  malalignment  e  will  experience  an  additional 
lift  force  v/hose  component  in  the  OXY  plane,  referred  to  axes  OX,  OY  is 
determined  by  * 


Lr.  =  krY^ 


i(^p+o-) 


(6) 


(2)  gives  the  normal  force  on  an  aerodynamically  symmetrical 
missile.  When  the  definition  of  Rp  is  extended  to  include  (6)  we 
have 


=  L  +  +  l^j  "Zl! 


by  (3).  (7) 


and 


3,24  fp  consists  of 

(i)  a  stabilising  moment  fy,  due  to  the  yaw  and  moment  maleilign- 

ment . 

(ii)  a  damping  moment  Tq  due  to  the  cross-spin. 


41. 


CONFIDENTIAL 


CONFIDENTIAL 


Technical  Note  No,  G.W.  177 


'  The  component  of  fy  plane  OXI  ic  determined  in  magnitude 

and  direction  by 


G'  J  Q 


(8) 


where  d  is  the  distance  of  the  C.-G-.  ahead -of  the  C.P. 

The  component  of  Tq  in  the  plane  OXY  is  determined  in  magnitude 
and  direction  by  ■ 


kQ  .  d  V  ^ 
dt 


(9) 


3*3  To  the  accuracy  of  A,1  the  components  of  along  OZ  have 

magnitudes  -R,  -f  and  in  the  plane  OXY^  referred  to  axes  OX,  OY,  are 
determined  by  -^RS,  where 


E  =  1%!  ,  r  =  I  .  (10) 


Also  the  components  of  Rp,  fp  along  OZ  are  of  the  second'  order,  and 

the  components  in  the  plane  OXY,  referred  to  axes  OX,  OY,  are  determined 
by  Rp,  fp,  where 


>.  ;  fp  =  Ty  +  Tc  •  ■  (11) 

if  Force  and  Couple  produced  by  Boosts 

if.l  Lot  the  boost  thraists  be  reduced  to  a  force  T  acting  at  the  G.&.  and 
a  resultant  couple  G.  Let  Tj^,Tp  and”  Gp  be  the  components  of  ^  and  G 

along  the  axis  and  perpendicular  to"^e*axis  respectively.  Then  to  the 
accuracy  of  A.l,  the  components  of  Tp,  Gp  in  the  plane  OXY,  referred  to 

axes  OX,  OY,  are  expressible  as 


iC^rp+O") 

Tp  =  TAgttj  o 


(1) 


V- 

Gp  =  TepOt^  G 


i(0(j+o-) 


(2) 


where  T,=  1  Tj^\;  ay,  are  the  thrust  and  couple  malalignment  angles; 
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^  are  the  thrust  ajid  couple  malalignment  orientations;  ^3  is  the 
distance  of  each  thrust  axis  from  the  C.G. ,  and  Ag  is  the  radial  offset 
of  each  nozzle. 


A-.2  To  the  accuracy  of  A.l,  the  components  of  G^  along  OZ  have 

magnitude  T,  G  and  in  the  plane  0X1  are  determined  by  T^,  G^  when 
referred  to  axes  OX,  OY;  where  G  =  [G^  [•. 

The  components  of  Gp  along  OZ  are  of  the  second  order,  and  the 

con^jonents  in  the  plane  OXT  are  determined  by  Tp,  Gp,  when  referred  to 
axes  OX,  OY. 

If. 3  Expressed  in  terms  of: 


w 

Q 

^N 

and 

T  and  G 


the  effective  gas  efflux  velocity 

the  total  rate  of  loss  of  mass  from  the  system 

the  distance  of  the  centre  of  each  nozzle  from  the  missile  axis 

the  offset  angle  of  each  boost  nozzle  axis  in  the  tangential 
plane, 

can  be  v/ritten, 


T  =  Qv/  ,  G  =  T  A|^. 


(3) 


5  Jet  dampinp;  couple  J 

3*1  This  is  the  restoring  couple  that  arises  Tvhen  the  missile  is'  rotating, 
due  to  the  additional  sideways  velocity  with  which  the  boost  gases  are 
ejected.  The  magnitude  of  the  axial  component  is  therefore 

J  =  Qkg  r,  (1) 


etnd  acts  in  the  opposite  sense  to  the  spin;  is  the  radius  of  gyration 
of  the  boost  exit  pianos  about  the  missile  axis, 

5.2  The  component  of  along  OZ  is  -J,  and  the  component  in  the 
plane  OXI  is  determined  by 

5.3  .9pmponent  of  the  transverse  jet  daii5)ing  couple  Jp  in  the  plane 
OXI  is  determined  by 


% 


Jp  =5  -  Q 


^GN^ 

dt 


(2) 


where  is  the  distance  of  the  centre  of  each  boost  nozzle  exit  plane 
from  the  centre  of  gravity. 

The  component  of  Jp  along  OZ  is  of  the  second  order. 

6  Equations  of  Linear  Motion 

6.1  These  are  obtained  by  equating  total  force  to  mass  times  accelera¬ 
tion  in  three  fixed  mutiially  perpendicular  directions,  chosen  -here  to  bo 
OX,  OY,  and  OZ^-1. 
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In  the  OZ  direction,  the  forces  acting  are  components  of  thrust,  drag 
and  gravity.  The  equation  of  motion  is  : 


mf  =  T  -  R  -  mg  sin  a, 


(1) 


v/here 


m  =  total  mass  of  missile  at  time  t 


f  =  acceleration  dV/dt. 


In  the  plane  OXY,  the  two  equations  of  motion  are  written  as  one 
in  terms  of  the  complex  quantities  introduced.  The  forces  acting  are 
gravity  and  components  of  Bp  which  have  already  been 

mentioned.  Accordingly  the  equation  of  motion  is 


m  A.  (VZ)  =  (T  -  R)^  +  Tp  +  Rp  +  mg  cos  a.  (z) 

dt 


6.2  The  appropriate  forms  for  T,  R,  etc.,  will  now  be  substituted  in 
(l)  and  (2) .  Writing 


a  =  T/m 


(3) 


(1)  b 


e  comes 


f=a“— V^-g  sin  a, 
m 


ih) 


and  by  equations  (l) ,  2,23(l),  and  3.23(7)  and  3.3(lO)j  equation  (2) 
becomes 


V  dZ/dt  =  alH:  +  (Tp  +  L  +  +  g(cosa+  Z  sincO.  (5)* 


When  ct  <4-5°  the  term  g  sin  d  Z  may  be  neglected  by  A.l.  When  45°  <  “  <90® 
wo  assume  that 


either  (i)  jz]  «  cot  a  A. 3 

or  (ii)  g  <<  a. 


*  The  inclination  of  the  trajectory  to  the  horizontal,  usually  denoted 
by  0,  is  equal  to  •  a  -  g  v/here  2  is  .the  real  part  of  Z.  The  component  of 
gravity  perpendicular  to  the  trajectory  is 

g  cos  0  =5  g(cos  a  +  S  sin  a  +  0(£^)) 
for  small  2.  ,  , 
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A. 3(i)  is  likely  to  hold  in  certain  cases  for  a  near  70°  80°,  but 

for  near-vertical  firing  (ii)  becomes  necessary. 

By  equations  2f.l(l),  3.23(4-),  and  3.23(6),  equation  (5)  becomes 


dZ 

ds 

fe-‘) 

icr 

■  a  e  g  cos  CL 

V2  v2 

(6) 

where  -6 

=  kj/m 

(7) 

(8) 

1-^  ^ 

=  AgOj 

e 

(9) 

and  |i 

Jj 

= 

• 

(10) 

The  first  term  on  the  right-hand  side  of  eqijation  (6)  arises  from 
the  transverse  con:5>onent  of  thrust  when  yawing  and  the  lift.  The 
second  term  arises  from  the  thrust  and  lift  maJ.aligments ,  and  the  last 
term  from  gravity. 

7  Equations  of  Angular  Motion 

7.1  To  the  accuracy  of  A.l,  the  component  of  angular  momentum  of  the 
missile  in  the  direction  OZ  is  Or,  and  the  component  in  the  plane  OXT, 
referred  to  axes  OX,  OY  is  determined  by  Cr^  +  i  A  d^dt;  C,  A  are  the 
moments  of  inertia  about  the  missile  axis  and  any  transverse  axis 
respectively. 

The  equations  of  motion  are  found  by  equating  the  rate  of  change 
of  angialar  momentum  to  the  couples  acting  in  the  three  directions  0X,0Y  and 
OZ^'2  viz. 


—  (Or)  =  G  -  r  -  J  (1) 

dt 


—  +  i  A  L  (G  -  r  -  J)?  +  Gp  +  Tp  +  Jp  .  (2) 

dt  I  dt  j  e  P  P 


7.2  By  equations  3.22(l),  3.3(lO)  and  5»l(l),  (l)  becomes 

d(Cr)/dt  +  (YrV  +  QkQ^)r  =  G  -  YpV^.  •  (3) 
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This  equation  gives  r  when 


V  is  knovm,  we  get 


where 


(4) 


(5) 


The  suffix  o  is  used  to  denote  the  value  of  a  quantity  at  laimch. 

7.3  The  reduction  of  (2)  to  a  convenient  form  is  quite  lengthy.  Firstly 
(2)  can  be  written  by  equations  (l),  3«3(ll)>  3«21.(9)  and  5«3(2)  as 


d^S/dt^  +  Xd^dt  =  (Gp  +  rY)/iA, 


(6) 


where 


AX  =  2Ak  -  i  Cr  +  kQdV 


(7) 


•  2AK  =  Q  +  dA/dt; 


(8) 


by  equations  A.l(2)  and  3.2A(8)>  the  right-hand  side  of  (6)  is 


(Gp  +  rY)/iA  =  a  iIq  0  -  n^V^ 

(9) 

where 

(10) 

iGpl  i^^G  T6b 

'  iia  "  =  -  "  IT  "0  “  =  -  n.  ^  = 

(11) 

=  -  nVaM  a^*«A 

(12) 

and 

n^  =  kpd/A  . 

(13) 
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Next  putting  5=  +  Z  in  (6)  and s eliminating  dZ/dt  using  6.2(6),  we 

obtain 


d^is  ^dX'_Gp  +  rY  ^  dz  .d^Z 
dt2  ^  dt  iA  “  ^  dt  “  dt2 


=  apo 


icr 


-  nV  S  -  \ 


a 6^°”  geos  a- 

-  +  - 

V  V 


<i  /  a  e 
l^p 


dt 


icr 


V 


g  cosa 


V 


2  J 


(15) 


Expressed  in  terms  of  V3I13  as  dependent  variable  and  of  s  as 
independent  variable,  the  left-hand  side  of  (14)  is 


ds  d^Cv^D  dCv.-^;) 

7Z2"+  ^  '7r  =  ^~7~2 —  +  (x  -  -  )  - 

dt  dt  ds‘^  V  d.s 


Xf 

Ll 

-v2 

dt 

w'^J- 

(v:iD  . 


(16) 


Equations  (15)  and  (16)  now  combine  to  give 


d2(V3t) 


2P'(s) 


d^(Y,:H) 

ds 


+  p(s)(v::i:)-=^T(s) 


gcoB  a  /f  \ 

+  ~yz  \Y  “  ’ 


(17) 


where  2P'(s)  =  I  +  X/V  +  (a  -  f)/V^  (l8) 

P(s)  =  n^  +  ^'  +l\/V  +  ?^a  -  f)/V^  +  d[(a  -  f)A2]/ds  (19) 

T(s)  =  T^(s)  +  T2'(s)  (20) 

Il(s)  =  a  e^  (>qip  -  V^^)/v2  (2l) 

T2(s)  =  a  e^'iipA  (22) 


and  dash  denotes  differentiation  with  respect  to  s. 
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Changing  the  dependent  variaLle  in  (17)  to  H  by  the  substitution 


V.'Z.'  =  H  e 


-P(s) 


(23) 


we  obtain 


d^H 

ds2 


G(s) 


V2  U  " 


where 

G(s)  =  P(s)  -  [P'(s)]^  -  P'(s) 
=  n^  -  A^(s)  +  A*  (s)  , 

v/'here 

2/V(s)  =  e  -  a/V  +  (a  -  f)/v2. 


(2O 


(25) 

(26) 


(27) 


Further,  by  eqiiations  6,2(7),  6.2(b),  and  3.23(5),  vi-e  can  write 


o 


(28) 


and 


2A(s) 


>  g  sin  a 
ra  V  v2 


P'(s)  -  A(s)  =  A/V. 


(29) 


(30) 


The  order  of  solving  the  eq\iations  of  motion  is: 

6.2(b)  for  V 
7.2(b)  for  r 

7.3(2i^for  K,  and  thence:^  from  7*3(23) 
and  finally  6.2(6)  for  Z,  ■ 
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8  General  Solution  by  integrals 

8.1  In  this  paragraph  the  equations  of  motion  are  solved  for  the  case 
when  G(s)  is  real  and  constant.  We  begin  by  deriving  conditions  under 
which  this  is  true,  Equation  7*3(26)  is 


G(s)  =  r?  -  A(s)  +  A’(s)  , 


where  by  equations  7*3(29)  and  7*3(7)  a(s)  can  be  e:jqpressed  as 


A(s)  =  +  'Tq  + 

where 


2aq  =  kp/m  -  kQ  d/A 

(2) 

=  -  (k  -  l^r) A  ;  3  =  C/2A 

(3) 

'1^2  =  g  sin  a/2V^ 

(4) 

It  is  clear  that  jT^q|  ’  and  |t2j  can  be  considered  negligible  only  for 
sufficiently  large  values  of  V.  It  might  be  expected  however  that, 
if  the  velocity  of  launch  is  low  and  the  acceleration  is  high,  the  time 
during  which  negligible  vd.ll  be  so  short  that  the 

solution  for  the  ensuing  motion  will  be  unaffected  by  omitting  them. 
This  is  confirmed  in  the  next  paragraph. 


8.2  8.21  Let  Vq  and  V, 


and  be  the  velocities  at  which  jxqj  =  l'’^2l^ 
jTTq’  j  =  1^2*  I  respectively.  It  is  assumed  that  Vq^  and  Vq^  are  imiqucly 
determined,  then  |Trq|  >  |T2  ]  V  >  Vq^  and  |  >  jxg'  |  for  V>  ^02* 


It  is  assumed  that 


B.l 


(“l  *  * 


is  negligible  cai5)ared  to 
for  all  V>  Vc^ 


B.2 


is  negligible  companed  to 
for  all  V  >  Vq^ 


8.22 


When  Max  (Vq^  Vq^)  we 


cjan  then  put 


A^(s)  =  aj^  ,  A’(s)  =  aq' 


(5) 
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in  7.3(26)  giving 

2  P 

Gr(s)  =  n  “  '  ' 

8.23  When  ^  Max  (Vq^  Vq^)  it  is  further  assumed  that 
B.3  g  sin  a  is  negligibly  small  con5)ar0d  to  a, 


Then  g  sina  is  negligible  in  6.2(4).  and  v/ill  not  appear  in  equation  7.3(29). 
In  other  words  the  motion  is. unaffected  by  neglecting  \'T2  \  inA(s),  and  |t2'j 
in  A'(s).  But  when  Vq^  |t2_|  <  Y2\  >  this  range  of  V 

■will  not  affect  the  motion.  Similarly  neither  -will  in  the 

range  • 

8.24  The  assumptions  B.1  and  B.2  put  a  restriction  on  the  magnitude 
of  the  spin;  when'  B.1  and  B.2  hold  the  spin  is  called  'small'. 

8.25  It  is  further  assumed  that 

B.4  kp/m,  kp/m,  kpd/A,  k^d/A  are  constant. 

Then  A(s)  is  constant  and 


G(s)  =  n^  -  =  P^,  say. 


(6) 


8.3  8.31  The  Solution  for  V,  by  6.2(4),  is 

o 

^  2kp(u-s  )/m 

-2kp(s-s  )/m  -2kp(s-s  )/m  e  (a  -  g  sin  a)  du 


V2  =  V  2  DV-  -o^ 


+  2  e 


(7) 


8.32  When  a  is  constant  (7)  gives 


-2kp(c-s  )/m 
f  =  fo  o  ^  ° 


(8) 


and 


(9) 
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2  . 


8.4  8,41  The  general  solution  of  7»3(24),  when  G(s)  equals  p  ,  is 


s. 


H(s)  =K^  cos  p(s-so)  +  K2  sin  p(s-Sq)  -  i.  /  e  T(u)  sin  p(&-u)  du, 


.  (10) 


where  K2_  and  K2  are  constants  depending  on  the  initial  conditions. 

The  contribution  to  H(s)  of  .the  gravity  term  has  been  omitted  as  the 
formulae,  for  gravity  drop  are  given  in  Hef.  1.1. 

Substituting  7*3(20)  into  (lO)  and  integrating  by  parts  we  have 

H(s)=»Kj^  cos  p(s-Sg)+K^  sin  p(s-Sg)-i  j  e^^“^[pT2(u)  cos  p(s-u)-T^(u)  sin  p(s-u)]  du, 

.  (11) 


where 

K3  =  K2  +  T2(sq)/p  (12) 

T^(s)  =  -Ti(s)  +  P'(s)T2(s).  (13) 

Further  it  is  seen  from  (ll)  to  (13)  that 

7.3(23)  (14) 

by  2.23(1),  6.2(6),  7.3(22)  and  (5)  K3  =  (l5) 

and  by  7.3(2l),  7.3(30),  7.3(22)  and  (5)  T3(s)  =  a  A- (16) 

Siibstituting  for  K3,  T2(u)  and  T3(u)  in  equation  (ll) ,  we  obtain 


H(s)  =  ^  cos  [p(s-s^)  +r)]  +  V2a  Ij2  sin  p(s-s^) 


P(u)+io-(u)rM,3  .  ^  n  r  r  \  il 

+  /  e  —  sin  p(s-u)  -  Up  ii  cos  [p(s-Sq)  +  t)  ]  ,  -  > 

Lp  P  J 


(17) 


51. 


CONFIDENTIAL 


CONFIDENTIAL 


Technical  Note  No.  G-.W.  l77 


where 


(18) 


^  ^2  =  Vp 


T)  =  tan  ^  “i/p 


(19) 

(20) 


Hence  aquations  (l7)and  7»3(23)  give 


TiSi(s)  +  +Cp  (s,(Jp)  +  ^c(s, 


(21) 


where 


.-P(s) 


S^(s)  =,  -  -  cos[p(s-s  )  +n] 

1  V//2E  P 


(22) 


e-P(s) 

^2(3)  =  -  sin  p(s“s^) 


(23) 


rP(s)  ^  P(u)+icr(u) 


p(u) 


p(s-u) 


audu. 


u 


(24) 


+  r  \  ^  P(u)+icr(u)  n  r  /  N  1 

5p(s,p)  - - -  /  p(u)  e  -  cos  [p(s-u)  +p]  —  .  (25) 

J  u 


V 


8.42  Differentiating  equations  2.23(1)  and  (21)  to  (25)  we 
obtain  for  d^/ds, 


d^ds  =  dZ/ds  +  d7Z^/d£ 


ds 


P'(s) 


e“P(s) 

^  ^  -7 -  [  P(3-S^)  +  T)  ] 


Y/ile 


p  e"^ r^G  P(u)+io-(u)  ,  \ 

Lg  - -  p  cos  p(s-s^;  +  =- -  /  —  e  cos p(s--iy  — — 


-P(s) 


+  n 


V/V2a 

g--P(s)  ^j-  p(u)+icr(u) 

M-  p  0 


V  /  P 


'u 


r  /  \  T  ^  ® 

sin  [  p(  3-u;  +  TlJ  -  -  - —  fip  ; 
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i. e.  by  equations  6,2(6),  and  l.J>{27) 


4  ^C^~1 

■? 


,  (26) 


where 


e-*P(s) 

^;j_(s)  =  -  cos  p(s-So) 

V/V^a' 


(27) 


e-'P(s)  n 

7(s)  =  -  -  sin  [p(s-3  )  +  Ti] 

^  V/V2a  P  ■ 


(28) 


+  .  N  ^  ^  n  .  r  /  N 

?q(s,p)  =  — -  /  P(u)  e  -  sin  [p(s-u)  +p]-:^ 


u 


(29) 


^  P(u)+ia-(u)  a^du 

%U>P)  = - /  P(u)  e  cos  p(s-u) 


V 


V. 


u 


(30) 


On  substituting  forHJ  (s)  from  equation  (2l)  and  rearranging,  equation 
(26)  becomes 


1  d^  _  T 
p  ds  ^ 


rocn 


f  51(3)  -  5|(3) 


+  L, 


rtt-, 


P 


?2^^)  +  Sq(s) 


p 


(s,Pf)  ?  q  (^jPp) 


i, e.  by  equations  (22)  to  (25)  and  (27)  to  (3O) 


n  d^ 
p2  ds 


(3,Pf) 


,  (31) 


where 


5£(3)  = 


-'P(s) 

e 


“  cos  [p(s-Sjj) 

wt' 


’ll 


(32) 


5;(3.h) 


e-g(3) 

' .( V 


P(u)  +icr(u) 
e 


^  cos  [p(s-u)  -p] 
p  ’u 
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The  solution  for  S  can  now  be  written 

__ — S.  =  ]^S^(s)  +  L2S^(s)  +  Sp(s^p)  +  ^s,  — ^  ,  (34) 


where 


s 


^2(s) 


s 


s 

j  Sp  (ujl^)  du 


Sp(s,n) 


s 


(35) 


(36) 


(37) 


(38) 


8.1.3  The  solution  for  Z  is  novf  obtainable  from  equations  (2l)  and 
(3I-) ,  using  the  relation 


z-z<,=  («-  «o)-C- 


(39) 


9  Assijmptions  required  for  evaluation  of  integrals 

9.1  In  order  to  er\raliiate  the  integrals  of  para.  8,  it  is  neoessaiy 

to  make  certain  additional  assumptions.  E^cpressions  for  the  yaw,  axis 
inclination  and  angular  deviation  due  to  boost  malalignment  can  then 
be  found  in  a  suitable  form  for  numerical  evaluation. 

The  assun^^tions  introduced  refer  mostly  to  the  aerodynamics  of  the 
round  and  are  valid  only  over  a  limited  range  of  velocity.  This 
restriction  is  implicit  in  3.4,  as  the  wavelength  of  yaw  is  constant  for 
a  range  of  subsonic  velocities  only.  The  angular  deviation  can  in 
general  be  evaluated  accurately  under  these  assumptions  without  intro¬ 
ducing  aerodynamic  terms  that  are  of  only  secondary  importance  and  in 
any  case  difficult  to  assess  reliably. 

9.2  It  is  assumed  tl:iat 

C.l  The  acceleration  is  constant  bety/een  launch  and  all-burnt 
and  denoted  by  a. 
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Then 


V  -  =  a(t  - 


(l) 


-  V  ^  =  2a(s  -  s  ) . 
o  ^ 


(2) 


If  the  instant  of  ignition  is  taken  at  t  =  0,  equations  (l)  and  (2)  give 


V 

o 


at. 


(3) 


V 


2 


o 


=  2as^ . 
o 


(4) 


t  is  the  'effective  time  of  launch'  and  s^  the  'effective  launcher  length'; 
they  are  defined  by  (3)  and  (if)  in  terms  of  the  launching  velocity 
and  acceleration  a.  The  'actual  launcher  length'  required  for  the  same 
launching  velocity  is  usually  somewhat  greater  than  so  because  the 
thrust  build-up  is  not  instantaneous  on  ignition.  The  word  'effective' 
is  often  omitted  when  the  context  is  clear.  Prom  (2)  and  (4)  we  have 


^  -  2as  . 


(5) 


9.3  It  is  assiomed  that 

C.2  neglected  in  comparison  with  Yj^Vr  or  G* 

C.3  C  and  G  can  be  taken  constant. 

Prom  7.2(4)  and  7.2(5)  it  follows  that  the  solution  for  the  spin, 
e^ressed  non-dimensionally,  is 

r/rg  =  (r>e)  [e(T)  -  2(To)]  ,  (6) 


where 


T^  =  Yj^s/C 


X  p 

r  u'^ 

E(x)  =  /  e  du. 
o 
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In  order  to  evaluate  the  integrals  containing  o*  in  8./+1  and  8,42  it 
is  assumed  that 

C.4-  the  range  of  integration  can  bo  dissected  in  such  a  way  that  the 
change  inYj^Vr  in  each  interval  is  negligibly  small  con^Jared  with  1  G-  -  Yj^Vr  ]. 

G-  -  Yj^Vr  can  then  be  replaced  by  G  -  YpV^r^  in  the  i-th  interval, 
where  Vj_  and  rj_  are  mean  values  of  V  and  r  over  this  interval. 

7.2(3)  nov;  leads  to 


r  =  cr  =  h^  ,  (10) 

Where 

hi  =  (G  -  YRV.ri)/G  .  (n) 


Integrating  equation  (lO)  we  obtain 


“  V  =  ''2i  “  ^oi)  +  nii  (/s  -  Vi^),  (12) 


where 


n2i  =  h./a 

(13) 

=  ^2/a  (r^i  -  t^^) 

(14-) 

and  3^^  and  t^j^  values  of  cr,  s  and  t  at  the  beginning  of 

the  i-th  interval. 


For  convenience  the  suffix  i  will  now  be  omitted. 
9.4  Finally  it  is  assumed  that 


C.5  lP(s)  I  is  negligible  compared  with  unity. 
C.6  (ij,  (iQ.,  and  are  constant. 


C.7  0^2.  negligibly  small  compared  with  p. 


C*5  is  the  condition  for  P(s)  to  be  neglected  in  the  solution  of 
para.  8.4.  It  is  effectively  a  restriction  on  the  values  of  s  for  which 
these  solutions  are  valid. 


C.6  is  justified  when  the  malalignment  angles  and  orientations  are 
statistical  measures  for  a  number  of  homogeneous  rounds:  for  an  individixal 
round  they  will  vary  in  an  irregular  and  unpredictable  manner* 

0.7  is  not  essential  for  the  evaluation  of  and  undertaken 

in  the  following  paragraphs,  but  will  be  adopted  at  this  stage  for  simpli¬ 
city.  Equation  8.25(6)  then  gives  n  =  p,  and  from  equation  8.2+1(20) 
we  have  n  =  0.  The  symbols  ^  g  and  ^  are  nov/  no  longer  distinct  in 
the  equations  of  para.  8.4. 
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lO  .  Solution  ■  for  the  yaw 

10.1  By  the  assun5)tions  of  para.  9  we  can  write  8.24l(2l)  as, 


(1) 


where 


cos  n(s-so) 

Ts 


(2) 


5,(3)  ^SiLilUiSoI 

/s 


(3) 


^ixr(u) 


cos  n(3-u)  dV, 


u 


irfu)  ,  > 

e  sin  n(s-u)  dV. 


(4) 


(5) 


cCu)  -  0*0  =  n2(u-SQ)  +  n^C/u  -  Vs^) 

2au 

=  Vj^,:n^V2a  = 


by  9.3(12) 

by  9.2  (5) 
by  8.41(18)  and  9.2  (4) 
by  8.11(19)  and  9.2  (4) . 


10,2  In  terms  of  the  yaw  integrals  (s>s^)»  I3(s,Sq)  defined  in  A.4 

and  of  the  functions  a  and  3„  defined  by 

s  s 


etc  =  , 
s  ^  % 


2(n2+n)s  n^ 

+ 


\/^(n2+n) 


2(n2-n)s 

P  =  /  - ^ -  +  — 

^  V  ^  N/2^(n2“"^) 


(6) 
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the  solutions  of  (if)  ajid  (5)  are 


§j(s) =  - 


%(s)  =  - 


icr„ 


2’/'s 


e^^o 


2i7s 


.7^  2 

-lo*  Og 
3  2 


e 


-in(s-s_)  -i^P. 


.  ft  2 


in( s-s. ) 


(7) 


Equations  (7)  can  be  written  by  A.  if.  2(l2f)-(l5)  and  A.4.l(3)“(5)  as 


2?j(s)  = 


e^o 


D(^)  ^d(Ps)' 


1  J 


i[o{s)-c5]p(ccso)  -in(s-SQ)  d(Pso)  ^in(s-so)' 

”■  G  ■f'  -1  -■  -T  --  r  0 

“1  ^1 


(8) 


2i?Q(s)  = 


el  O'er 

i[o-(s)-Og]  r 

7s  L 

-OCI  8i  _ 

D(«So)  -in(s-So)  D(3g^)  _in(&-Sor 
■  ■ '  '  G  ••  G 

8l 


J 


(9) 

This  form  of  the  solution  is  convenient  vfhen  oCg,  3^  are  real  and  positive. 
In  other  cases,  i.e.  when  Og,  3s  negative  or  imaginary,  the  most 
convenient  forms  of  (S)  and  (9)  are  readily  found  by  putting  the 
alternative  expressions  for  given  by  A.2.l(if)  to  (8),  in 

equation  (7).  For  exaiiple,  when  3^^  is  negative  and  s  <  s^, 

(8)  and  (9)  vrould  be  expressed  for  evaluation  as 


2^t(s)« 


XCTr 


Vs 


L-  «! 


Pt'  -J 


^i[o'(s)-crJ  D(as^)  ^-in(&-So)  d(-3*s^) 

^  “1  3l' 


(10) 


elo'of 

P(“s)  F(-3s')n 

_  _  .  JL  .  . 

iKs)-^] 

r 

7s  1 

-  “1  ^1  - 

W  ^ 

3l' 


(11) 


10,3  Constant  spin 

Putting  112  zero  the  equations  (6)  become 
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We  shall  now  find  an  approximate  solution  under  the  following 
assungitions : 


(a)  Vs  <<Vs^ 

(h)  •  >>  1. 

(a)  restricts  the  range  in'which- this  approximate  solution  is  valid; 
the  larger  the  spin  the  greater  the  range  of  Validity.  Equations  ( 
and  (13)  are  then 


“s  “  Ps'  nq/V^  , 


(14) 


The  object  of  introducing  (b)  is  to  allow  the  function  D(u)  to  be 
replaced  by  the  first  term  of  its  asymptotic  expansion  (see  Ref.  I.4) > 
namely  -  • 

D(u)  = 


The  error  in  (15)  'is  about  VJfo  at  u  =  1_,  2^  at  u  =  2  and  quickly 
decreases.  Like  (a),'  (b)  requires  that  the  spin  should  not  be  too 
small.  By  equations  (14)  and  (15)  we  then  have 


u  .  71:  2 
/*  1 

/  e  ^  dx  I  ~i/7ai.  (15) 


-P(“s)  +  P(-P3')  =  0 


(16) 


D(«J  -  = 


2i 

V  “ 


(17) 


Substituting  (I6)  and  (17)  into  (lO)  and  (ll)  we  obtain 


^t(s)  = 


1  e  ^ 


nq  Vs  L 


iIo'(s)-o-J 


-cos  n(s-SQ) 


(18) 


ioi 


1  e 


Cg'(s)  =  - —  sin  n(s-s  )  , 

ni  Vs 


(19) 
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11  Solution  for  the  axis  inclination  ^ 

11,1  Under  the  assumptions  of  para,  9)  equation  8.42(3it-)  becomes 


- =  Tq,  ^i(^)  +  ^  Sq.(s), 


(1) 


where 


^l(s)  =  -  /  ^2^ 


(2) 


Sp(c)  =  /  ?-|(u)  du 


(3) 


^t(s)  =  f  ?g(u)  du 


(4) 


s 

Sg(s)  =  -  J  5t(^) 


(5) 


11*2  ^(s)  will  now  be  evaluated.  Write 


^^(3)  =  ^i(s)  +  i  ^2(3) 


-- 

V  a 


in(u-so)  /  s  /  s 

o  dV^  by  10.l(2)-(3), 


(6) 


then  this  integral  can  be  evaluated  by  putting  =  6^,  where 


6  = 
u 

6^  Vu  + 

(7) 

61^  = 

2n/r. 

(8) 

^o  = 

0, 

(9) 

A.2.1l(l)  and  (3). 

It  is  then  found  that 

1 - 1 

tn 

D 

0 

T,/c  X  in(s-So) 

-  D(Ss)  ® 

(10) 
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Since  8g  is  always  positive,  equation  (lO)  is  in  a  form  convenient  for 
numerical  evaluation  for  all  values  of  s. 


11.3  Sfji(s),  ^q(s)  are  found  by  substituting  10.2(8)-(9),  into  (4)  and 

(5).  In  terms  of  the  integrals  J(j(s,Sq),  Jj3(s,SQ)  of  A. 3  aj^d  S^(s),  we 
have 


2iS^(s) 


-  «1  Pi  - 


^(Ps  ) 

o 

Pi 


(11) 


2  Sp(s)  =  -2 


■Ja(So»s)  JfiCs  s)- 


L  ap 


3l  J 


+  1 


- 


L  «! 


(s) 


10"^ 


•  (12) 


J^(s^,s),  Jp(sQ,s)  can  be  found  by  one  of  two  approximate  methods 
according  as  either  n-]_,  i.e.  spin  is  sufficiently  large  (the 
asymptotic  solution)  or  n2,  i.e.  angular  acceleration  r^  is  sufficiently 

large  (the  approximate  solution).  The  appropriate  formulae  are  given 
in  A. 5  together  with  estimates  of  their  accuracy. 


It  is  found  that  v/hen  (or  V^jj)  is  positive  and  greater  than  7i^ 

there  are  two  expressions  for  the  asymptotic  solution  according  as 
Sq  <  s  <  s^  (or  Sp)  or  G  >  s^  (or  Sp),  For  exaitple,  suppose  Vs^  is 

positive  and  greater  than  /s~q  and  that  is  negative  (this  is  case  (d) 
in  A-.5),  then  by  A. 5. 35(3^)  A, 1.1(3)  we  find  that 


2  e-i°'o  J3(so,3)+iD(3s  )S"(s)=  (-2) 


T5(-3^')  e 


i[a-(u)-cro] 


2  n2  Vii  +  n-j_ 


+  0  ^  (s) 


(13) 


when  s  and  by  A. 5. 35(35),  A. 5. 34(31),  A. 1.1(3)  and  11.2(7) 


2  e  °  Jp(so,s)+iD(3s^)S^(s)  =  2 


^^(-Ps  ') 


2n2  s!  +n^  ^^2^  s  +  n-]_ 


+D(-|3s„«)S  (s) 


+2D(o)  °°^j^2(n2-n)/nni  +  (2i/6i)  D(8 


J 


dO 


when  S(3  ^  s. 

Substituting  in  (ll)  and  (12)  for  Jct(sQ,s)  from  A, 5. 31(13)  for 
J3(^o>s)  fi’om  (13)  and  (lA)  we  obtain 
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Er(u)-crJ 


“(%„)  C(-e  ') 

-1  rr^«(3) 

L  “l  I3i 


(15) 


-1ST 


2e  %(s)  =  2i 


2i< 

‘D(a^)  D(H3^') 

i{o-(u)-cro] 

e 

“1  h'  _ 

2n  Vu+  nn 

2  1 

>  +  i 


D(a D(-3  •) 

-  “1  3i 

. (16) 


when  Sq  <  s  <  ,  and. 


-icr 


2ie  °^(s)=-2i 


iXttg)  D(^') 

i[°'(s)-‘y,] 

^  +  2i 

!>(%  )  ) 

0  ^  .  .0 

“1  Pq' 

T  SLJl. 

2n2^s+ni 

«1  '  Pi' 

2n2^^o+»l 


-  1 


ri^(«s  )  _  d(-3s  ') 

- - 2_gX(3) 

«1  3i’ 


+i2D(o)  “ 


i['^(s3)-o-J  _  in(s-se) 

;2(i^-r)/nn^+(2V5i)  D(6^e  -'^(Ssi^) 


3l'  L 


(17) 


2  e“^'^o?(j(S)  =  2i 


D(a^)  D(3^'y 

“1  "'I?' 


i[c-(s)-o-Q] 


2n2/ s  + 


-2i 


D(as  )  d(-3  ') 

_ O _ 

^1  "  Pi* 


2no  +  n-, 
2  o  1 


+  1 


'D(as^)  IJ(-^  -) 


^(s) 


Pi 


—  S^(s) 


+i2D(o)  - 


i['^(sp)-Oo] 

'h' 


j^2(n2-n)/nni+  (2i/6i)  0(63)®^^^"°^ -iX^sp)  j 


(18) 


when  <  s. 

12  Solution  for  the  angular  deviation  Z 

12,1  By  equations  8.43(39),  10.1(1)  and  11,1(1)  we  have 
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Z  -  Zo  =  Li  Z^Cs)  +  L2  Z2(s)  +  Pj  Zj(s)  +  ■—  Zq.(s)  , 

(1) 

where 

Z^(s)  =  n?2.(s)  -  [5i(s)  “  ^iCsq)] 

(2) 

Z2(s)  =  n  S2(s)  -  52(3) 

(3) 

Z0(s)  =  n  ^0  (s)  -  ^0  (s)  for  6  =  T,  G  . 

(^) 

12.2 

Zero  Launching  spin,  small 

12.21  In  this  case  r^  is  zero  and 

n^L  =  -  2n2fsQ  . 

(5) 

Hence 

eq.uation  A.5.52(A7)  becomes 

=  v/2n2A  (Vu  -  Vi^)  , 


(6) 


showing  that  is  zero  initially  and  positive  for  u  >  s^. 
By  eq^uations  A.  5.52(48)  -  (49),  we  have 


n 


n2  +  n 


^  ^o 

% 


(7) 


n  p  ^  Sq 

n2  -  nV 


(8) 


It  is  seen  from  eq.uations  (7)  and  (8)  that  is  negative  and  is 
positive.  The  arguments  of  D(Y^)  and  D(Yu  -  in  A.5.52(2f5)  are 
therefore  positive,  but  the  argument  of  D(Y,^  "  in  A. 5 .52(1-6)  is 
negative  for  values  of  s  near  Sq.  Provided  -kq  is  greater  than  -O.5 
the  error  in  A. 5.52(A6)  will  be  less  than  20/o  by  A.5.31-. 

.  .  Prom  A.A.1(12)-(13)  we  have 


Vs„  =  n2Vs^(n2  +  n) 

(9) 

=  ^2'^^(n2  -  n)  . 

(10) 

•  63. 

CONFIDENTIAL 

CONFIDENTIAL 


Technical  Note  No.  G.W.  177 


(9)  shows  that  <  Vs^,  v/hich  means  that  >  0  for  s  >  s^.  (lO)  shows 
that  Vs^  is  slightly  greater  than  so  that  in  the  range  -/s^  <  Vs  <  Vs^ 

will  be  negative. 

12.22  Setting  5  ~  we  have  the  following  results 


=  0 

by  10.1  (1)  -  (5) 

(11) 

=  (2i/6i)  D(6so) 

by  11.2  (10) 

(12) 

Q-i^'o  Ja(sg,c«)  =  i(n2  +n)  [D(o)  -  D(-Ka)]  /  ^^'^0 

by  A. 5. 52  (45) 

(13) 

e"^o  jp(s^y«)  =  _i(n2  -n)[D(o)-  D(-Kj3)]  /  2nn2Vo^ 

by  A.  5. 52 (46), 

(14) 

where  are  given  by  (7)  and  (8). 

When  s  is  zero,  and  Kq  are  zero  by  equations  (7)  and  (8), 

and  equations  (13)  and  (14)  becomes 

Ja(o,~)  =  e"^°o  J3(o,«)  =  i/V^2 

by  A.3.32  (4L). 

(15) 

The  solution  for  and  are  obtainable  from  11.3.  (H)  “  (l2)  ; 

they  are 


■io*, 


i  e  °  ?t(“)  =  e 


_  3i  J 


D(«So)d(5so)  d0so)d(5s„) 
L  ai6i  3i5i 


(16) 


-ice. 


p(-So)l>(6so)  D(e3„)I>(6s„)-l 

L  «i  '  3i  J 

L  ai6i  3i6i  J 

,(17) 


where  J3(sq,co)  are  given  by  (13)  and  (14)  or  by  (15),  and  v/here 

IOsq)  can  be  e:xpresse.d  for  the  purpose  of  evaluation  as 


c(3so) 


-i  Tt  3  2 

e  2  (1  +  i) 


D(-Ps,) 
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By  (4)  and  (ll) ,  Z0('»)  is  now  given  by 


Z0(«) 


(18) 


for  0  =  T  and  G. 

12.25  It  in  sho™  in  A. 5.5  that  the  maxinium  error  in  (I3)  and  (lif) 
is  n/2n2  and  that  in  moot  cases  it  is  likely  to  be  considerably  less.  If 
v/e  sacrifice  a  certain  amoiint  of  this  accuracy  by  neglecting  n  compared 
Y/ith  n2,  (16)  and  (17)  can  be  simplified.  More  e^qjlicitly,  if  we  assiame 
that 

(a)  ti/n2  can  be  neglected  compared  v/ith  1 

(b)  (=u)  is  so  small  that  2A(+u)-1  and  2B(+u)-1 

\]n2\l  %  ~  ~ 

are  negligible  compared  with  1, 
then  by  (a)  we  ‘have 


“1  ==  <^1  = 

(19) 

a3  .  -P.  ..  p 

^0  V  n2  V  7^ 

(20) 

/rT 

"a  ’■ 

(21) 

and  hence  by  (b)  and  (20) ,  (2l)  v/e  see  that  ©(-k^j)  and  D(-K|3)  as  v/ell 
as  D(a3^)  D(-Ps^)  can  be  replaced  by  D(o).  Equations  (l3)  and  (l/f)  then 
become 


e‘‘^°‘o  00) 

By  equations  (l9)  and  (22) , 

0"^°°  ^(00)  =  . 


=  e  ’lp(so>®® 
equations  (I6) 


(1  +  i)7C 

2/n^ 


) 


) 


=  i/V2n^  . 

and  (17)  become 


(22) 


(23) 


0 


^2  2Vn2n 


B(^Sq)  • 


(24) 


Equations  (23)  and  (24)  become  exact  as  Condition  (b)  above 

requires  n2  to  bo  fairly  large  because  B(u)  clianges  rapidly  with  u  at 
u  =  0.  Equations  (I6)  and  {17)  on  thOj^othor  hand  would  give  quite  good 
accuracy  even  for  n/n2  eqimil  to  V3  or  V2. 
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12.24  V»hen  cr^  is  taken  as  zero,  the  solutions  of  Zrp(oo)  and  Zq(«<>) 
are  by  eq.uations  (l8),  (23)  and  (24), 


R1  Zj(c)  =  Im  Zt(co)  =  ly^  4(63^) 


(25) 


R1 


Zg('^) 


n  2Vi^n 


2=r  B(6  ) 


(26) 


Im 


Zpfco) 


B(6^  )  -  —  . 

2^1  "o  ^2 


n  ZVnIsi  "o'  n, 


(27) 


For  large  n2,  the  real  and  imaginary  parts  of  Zq.(«’)  are  equal;  we  have 


1Zt(c)i  =7./^  4(633) 
'  2nr 


(28) 


IZsC-)! 


% 


n 


B(SsJ 


(29) 


12.3  Constant  Spin 

12,31  V/hen  <  s  <  Gj3  the  solutions  for  Z^(s)  and  Zq(s)  follow  from 
(4),  10.2(l0)-(ll),  and  11.3(l5)-(l6)  v/ith  n2  =  0. 


i[o<s)-o-o] 


F(aso)  N(-^s^')■ 

L  «1  3l'  . 


1  i  sin  n(s-SQ) 


/s 


+  i  n  ^(s) 


N(aso)  B(-3s^’) 
^l' 


*—  a 


cos  n(s-So) 


/s 


-  n  ^{s) 


(30) 
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2e  ^(s)=  i 


1\  ^E(-Ps')A  1  y|^l[0-(3)-0-J 

"Ov/s  n~J  ■%'  V/s 


+  1 


_  “1  ""  h'  _ 


1  i  sin  n(s-sJ 


+  1 


I)(asJ  D(-Pp  ') 


_  “i  Pi  j  L  ^ 


QOS  n(s-So) 


+  n  4  (s) 


(31) 


where 


0^  =  =  -^2^% 


V^=  -  =  -  n^/2n 


The  functions  Pg'  are  defined  by  10.3(l2)-(l3) . 

12,32  For  values  of  s  somev/hat  less  than  sp  it  is  possible  to 
simplify  (30)  and  (3I)  by  using  the  asymptotic  expansion  of  D(u)  ,  namely 


D(u)  = 


2  3 

%  u 


+  ^  f  1  - 

TOi  I  7t2uZt- 


(See  Ref. 1.4) 


1.  e. 


D(u)  =?  iAu  v/hen  %x?  »  1  . 


If  v/e  assume  that 

(a)  s  lies  in  the  range  for  v/hich 


o 

2n  (■'^Sp  -  /s)  >>  1  (say>  4), 


the  asymptotic  expansions  of  D(-p'  )  and  DCdg)  are  valid,  and  we  have 


(/s  {j^ 


- 

V^/  ni 


i _  /2n 

VsJ  7; 


(32) 


A  f\arther  simplification  is  obtained  by  assxoming  that 


(b)  , 
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then 


-i-jf 


(33) 


by  the  asymptotic  expansion. 

Substituting  (32)  and  (33)  into  (30)  and  (31)  vfe  get 


Z,(s)  =  ^ 


cos  n(s-So) 

r^s - Vs - " 


(34) 


e“^^o 


Z^i(s)  = 


1  i  sin  n(s-so)  . 

vr - v^ — ^  ^  ^  ^  ^2(2) 


L''^P 


(33) 


Putting  cr^  =  0  and  dividing  (34)  and  (35)  into  real  and  imaginary  parts 
we  have 


R1  Zj(s) 


sincr(  s) 
npVs 


(36) 


Im  Zrn(s)  =  — 

n. 


cos  n(s-So)  -coso(s)  ^ 
Vs 


Si(s) 


(37) 


R1 


2iG(s)  2 


n  ^1' 


(38) 


Im 


Zg(s) 


n 


^1 


^2(3)  - 


sin  n(s-SQ) 
nVs 


(39) 


where  ^(s),  ^(s)  are  given  by  11.2(l0),  i.e.  by 


Sl(s)  =  V2t^  [-A(6s^)  +A(6s)  cos  n(s-S|^)  +  B(6g)  sin  n(s-SQ)]  (40) 


^(s)  =  ■/2n/n  [b{5s^)  +  A(5^)  sin  n  (s-s^)  -  B(6s)  cos  n(s-3j^)  ].  (4L) 
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Equations  (39)  and  (ll)  combine  to  give 


im  -  ''ii 

n  V  n 


n-, 


(42) 


where 


8(63^) -13^(63)  Din|  (832-  637)-  8(63)  003  ^  (63^-837)  (w) 


and 


A,  (u)  =  -i-  -  A(u)  . 

7CU 


12.33  It  can  be  shown  that  for  sufficiently  large  s 

(44) 


12g(s)1  .  Zg(s) 

- —  T  Im - 

n  n 


This  follows  because  A2(u),  B(u)  are  rapidly  decreasing  functions  such 
that,  when  5^  has  increased  to  about  1  or  2,  A2(63),  B(6g)  are  negligibly 
small  compared  to  ^(63^)  (which  is  or  a  little  less).  We  then  have 


and  hence 


n 


n 


% 


^lv2n 


Rl 


Zg(s) 


/  Im 


n 


Zq(s)  _  2  fin 
n  “  ni 
^/2/N 


because  n]L  satisfies  A. 5. 41  (36),  where  N  is  large. 
Hence 


Im 


Zg(s)  12g(s) 


n 


<  Im 


Zg(s) 


n 


n  \  n2  y 


which  proves  (44) . 
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12.34-  When  sg  <  s,  the  solutions  for  Zj(s)  and  Zq(s)  can  be  obtained 
from  (if)  by  substituting  equations  10.2(l0)-(ll)  and  11. 3(l7)“(l8)  ir.  it. 


-lOX 


2  e  Z^s)  = 


1  \  dOs')/1  1  -n  i[q'(s)-cro] 


_  “l  v/s^  7^3  y  Pq'  \/s  '^rJ_ 


p(%)  ^(H3so')n 


L  aq  3^,  _j 


1  i  sin  n(s-S3) 


+  1  n 


/s 


%(=) 


B(“so)  ^(-gs^'nrios  n(s-So) 

_  “q  3  t  JL  /g 


-  n  K^(b) 


3q'  11/^  6i 


^in(u-Sj3) 


(4-5) 


-Kr, 


2  e  °  ZJs)  =  i 


)(ag)  /I  1  \  dOs')/1 


“l  vs  \/s  Vs^/_ 


i[a<s)-cr^] 


+  1 


I)(as^)  D(-fe^')" 


—  <^1  Pi'  -J  I—  Vip  7s 


1  i  sin  n(s-s^) 

+ - ^  -in  ^(s) 


.+  1 


I)(aso)  D(-Pso’)' 

-  Pi-  J 


COS  n(s-So) 

Vs 


+  n  Sj^(s) 


2xE^  ^iKspjxrjf  1  j  ^in(it.sp)T  _ 

Pi'  ILVu  &!_  J  J 

^  P 
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12.35  As  s  tends  to  infinity  (k5)  and  (2f6)  give 


2  e^^^o  Zt(«)  = 


E(“sJ  i^(--Pso'r 


L  <xi 


3l' 


+  i  n 


"D(“3o)  D(-3so')' 


L_  a- 


n 


Si(«) 


Pi' 


+  i  D(  63^ 


(A7) 


2  e“^°'o  Zq.(co)  s  -  i 


T^(<^so)  D(-feo'^ri 
-  “1  Pi' 


+  i  n  S2(") 


+  1 


D(‘"so)  ^  D(-Pso') 

-  «1  3i’  J 


n 


-^D(o) 


Pi 


+  i  ^(63^ 


l7-.6s^  • 


(A8) 


Substituting  ^ 2n/%  and  assiaming  as  before  that 

<<  equations  (A7)  and  (48)  can  bo  simplified  to  give 


e  o 


Zj(-)  =  i£^-(l-i)  5  [Aj^(S^^)  *  iBCS^ij)] 


(49) 


^  ^  A),  a..)  i  .  ^(a33)]e^"<^3>-o] .  (30) 

X 
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It  is  now  shov/n  that  when  (49),  (50)  are  vail  id,  i,e,  when 
4.5.41(56)  holds,  the  A]_(£^p)-,  3(6^  )-teriiis  are  small  compared  with 

tlie  A(6sq)-,  B(63^)-terras.  By  condition  A. 5. 41(56)  we  have 


^  >  —  >  2  for  N  >5, 


and  hence  B(63  )  =  I'Y  the  asymptotic  expansion  (Ref.  1,4)^ 


Therefore 


and  hence  is  very  much  less  than  4(63  )  B(63  )  (which  are  usually  just 

less  than  ^) .  Likewise,  for  63p>  2,  A]_(63p)  is  less  than  B(6sp)  and 

so  is  very  much  less  than  A(6s  )  B(63  ),  Equations  (49)  and  (50) 

o 

then  give,  on  taking  cy  =0, 


R1 


Im 


Z(^  (00) 

=  0 

Z(ji  (00  ^ 

=  ■/2'TOi 

4(63  ) 
0 

^1 

Zg(~) 

n 

2 

^1^ 

(small) 

B(6s^) 

n 

v  n 

ni 

(51) 

(52) 

(55) 

(54) 


<  sp  not  satisfying  condition  12 *32 (a)  it 
\  and  the  terms  in  ^ol  and  in 


12.36  For  values  of  s 

is  usually  found  that  in  (30)  and  (31 )  the  terms  in  e 
■^l(^s)>  small  compared  with  A( 63^)-,  ^(63  )-terms. 

terms  are  neglected  the  solutions  reduce  to  (51)  -  ^54)*  The  same  is 
true  of  (45)  Q-nd  (46)  when  s>  s^;  and  so  (36)  -  (39)  can  be  taken  to  hold 
for  all  s.  The  error  in  these  equations  arises  mostly  when  s  is 

near  Sq  since  this  is  \hen  the  approximate  terms  are  largest.  For 

such  values  of  /s  (i.e.  «  Vs^)  it  is  shown  in  A. 5. 41  that  the  error  is 
less  than  i.e.  on/v^ . 


and  in 
When  these 
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Pi\RT  III 

A  SOLUTION  FOR  THE  MOTION  i\FT£R  BURNING 
1  Introduetion 


Part  III  of  this  note  is  eoncerned  with  the  flight  of  a  missile 
in  the  interval  between  separation  of  its  expended  boost  and  the  first 
application  of  control.  IXiring  this  time  dispersions  may  arise  from 

(i)  faiilty  detachment  of  the  boosts,  and 

(ii)  aerodynamic  asymmetry  of  the  round. 

The  effeet  of  (i)  ean  be  evaluated  in  terms  of  the  initial  conditions, 
i.e,  the  initial  yaw  and  initial  angular  velocity  of  the  missile;  the 
effect  of  (ii)  can  be  evaluated  in  terms  of  the  lift  and  moment 
malalignment  angles.  These  causes  of  dispersion  have  already  been 
introduced  in  Part  II,  and  so  little  further  v/ork  is  necessary  to 
obtain  the  equations  of  motion.  In  fact,  under  assumptions  for  which 
the  wavelength  of  yaw  is  constant,  the  solution  by  integrals  given  in 
11,8  can  be  taken  over  at  once  with  the  thrust  equal  to  zero. 

V/hen  the  spin/veloeity  ratio  is  constant  the  yaw  and  angular 
deviation  ean  be  evaluated.  Finally,  the  non-transient  terms  in  the 
angular  deviation  are  considered  under  simplifying  assumptions. 

2  The  equations  of  motion 

The  equations  of  motion  after  burning  are  found  by  putting  T  and 
hence  a  equal  to  zero  in  the  equations  of  motion  during  boost,  namely 
in  equations  11.6.2(0  and  (6),  11.7.2(0  and  (5)  and  11.7.2(20.  It 
is  understood  that  the  parameters  in  these  equations  now  refer  to  the  dart 
alone.  We  obtain 


f 


k-r\  p 

—  “  S  sin  a 

m 


(1) 


ds 


^j[^iL+^(s)] 


+ 


g  cos  a 


v2 


(2) 


I(t)  r  -  r 

o 


J. 

o 


l(  t)  -I -  dT 


(3) 


where 


i(t) 


Y  Yrr- 


dT 


C 
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G(3) 


H  = 


g  COS  CL 

_  v2 


(4) 


The  suffix  o  now  denotes  the  value  of  a  quantity  at  the  end  of  the 
separation,  i.e.  when  contact  between  all  boosts  and  the  dart  has  broken. 
By  11.7.3(26),  (29)  and  (30) 


G(s)  =  n^  -  A^(s)  +  A'  (s) 

(5) 

2a(s)  =  kp/m  -  4/v  +  g  sin  a/V^ 

(6) 

2P'(s)  =  kp/m  +  X/V  +  g  sin  oc/V^, 

(7) 

and  by  11.7.3(7)  and  (8), 


\  -  i(c/A)r  +  k^dV/A 


(8) 


We  shall  also-v/rite 


2cc^  =  kp/m  -  k^d/A  as  before  in  II. 8.1  (9) 


2a2  =  kp/m  +  k^d/A  .  (lO) 


3  General  Solution  by  Intepirals 

3.1  The  procediire  here  is  similar  to  II. 8;  assumptions  are  introduced 
which  permit  G(s)  to  bo  taken  constant;  then  2(a)  can  be  integrated  and 
the  yaw  obtained.  It  is  assumed  that 

D.I  kp/m,  kp/m,  kpd/A,  kQd/A  are  constant,  and 

D,2  A(s)  is  constant  and  denoted  by  A* 

Then  by  2(5)  we  have 


G(s)  =  n^  - 


p 

=  p'^  a  constant. 
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It  is  readily  seen  from  2(6)  and  2(8)  that  1.2  ijiiplies  the  follov/ing 
alternative  assumptions 

(a)  cither  (i)  spin/velocity  ratio  constant, 

or  (ii)  spin  sufficiently  small. 

(b)  either  (i)  velocity  constant, 

or  (ii)  sinCt  sufficiently  small. 

In  practice  (a)(i)  and(bKi)  hold  for  a  short  interval  after  separation, 
while  (a)(ii)  and  (b)(ii)  will  almost  alv/ays  be  satisfied  even  when 
sin a  =  1. 

3.2  In  view  of  D.l  the  following  solution  for  V  is  obtained  from  2(l) 


V‘ 


■2  =  V  2  _  [i  _  Q-2kD(&-So)/m  j 


mg  sina/kp 


(2) 


Hence 


f 


=  f 


o 


(3) 


and 


V2  -  [1  -  .  (4.) 


3.3  The  solution  for  H  is  obtained  by  the  method  of  II. 8.4,  v/hich 
will  not'  be  repeated  here.  The  only  malalignments  now  present  are  due 
to  the  aerodynamic  asymmetry. 

Prom  equations  II.8.U  (l8)  -  (25),  11.6.2  (8)  -  (lO)  and  II. 7. 3 
(lO)  -  (l2) ,  we  have 


ivl 


C(s)  +  -^  g  (s)  +  C.  ^  (s),  (5) 


where 


£,^3)  H  =03[p{3-S„) 

V  P 


(6) 


?2(s)  =  -2 -  sin  p(s-s  ) 


(7) 
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o-P(=) 

V 


s 


t 


P(u)+icr(u) 

e 


—  003  [p(s-u)+'n]V^  du 
o 


(8) 


%(s) 


V 


/ 


P(u)+icr(u)  ,  , 

e  sin  p  (  s-u)  du 


So 


(9) 


(10) 


(11) 


N 


tan“^  a/p  . 


(12) 


The  angular  deviation  can  be  obtained  by.  integrating  equation  2(2), 
where  the  ya\r  is  given  by  equation  (5) . 

If  *Spin  proportional  to  velocity^  solution 

4.1  Wo  nov/  make  the  following  assumptions  under  which  it  will  he  proved 
that  the  spin/Velocity  is  always  constant. 

E.l  kym  +  g  sin' a/V^  constant 


fi.  2  constant 


E.3 

At  the  end  of  separation  -i. 

(  ^')  =  0,  i.e. 

=io . 

(1) 

ds 

\V 

r 

V 

0 

0 

4.11 

E.l  implies  either  that  the 

change  in  V 

is 

small 

or  that  the 

drag  kpV^  is  very  lauch  larger  than  the  component  of  gravity  mg  sin  a. 
2(1)  gives 


V  =  V„  (2) 

where 

6  =  kp/m  +  g  sinaA^  .  (3) 
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4.12  By  £.2  and  (2),  2(3)  gives  the  solution  for  the  ^in, 


[r(YR  -  06)  +  Yj,v]  =  [r„(YR  -  06)  +  YpVj  , 


and  this  equation  can  be  written  as 


vVv  r/ 


ih) 


=  0  by  equation  (l) 


Therefore 


r/V  =  r^y^  =  Y  say 


and 


<^(3)  -  O'o  =  Y(s-So) 


(5) 


4.2  Comparing  2(6)  and  2(7)  it  follows  by  D.2  and  2(8)  that  P'(s)  is 
also  constant.  Write 


P'(s)  =  Pq  +  ip2  (^) 


where 


Pi  =  02  +  g 

sin  a/2V^^ 

(7) 

P2  =  -3Y 

,  P  =  C/2A. 

(8) 

By  (2)  and  (6),  we  have 


(9) 
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v/here 


B(3)  =  (p3_'  +  ip2)(s-SQ)  (lO) 

Pi'  =  Pi  “  S  (ll) 


where 


“2'  ”  “2  “ 

=  (kp  -  k33)/2m  +  k^^d^A  by  11.3.23(5)  (l2) 

4.3  The  solution  for  ;:zr  ,  3-3(5),  can  be  written  using  (9)  as 


^2(3)  +  Cl  gj;(s) 


P 


(13) 


where 


-B(s) 


H  cos 
P 


[p(s-Sq) 


+  p] 


(14) 


2(3)  =  sin  p(s-Sq) 


(15) 


gj(s)  =  -  e 


rB(3) 


s 

f  e^^  )+i‘^(^)  n  j,Qg  [p(s-u)  +p]du 
J  P 


(16) 


%(s)  =  -  e 


■B(3) 


f  B(u)+io-(u)  .  /  \  n  /t,\ 

/  e  sin  p(s-u)  du  (.17) 


and  o'(u)  is  given  by  (5) . 
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The  functions  5^(s)  and  5jj(s)  will  now  be  evaluated.  Write 

Hq  =  Pi'  +  ih^  =  pi'  +  i(p2  +  Y)  (18) 
Hi  =  pi'  +  ihi  =  pi'  +  i(p2  +  Y  -  p)  (19) 
H2  =  Pi'  +  ihg  =  Pi'  +  iCPg  +  +  P)  (2^)) 
gl  =  p  -  iA  =  n  e“^^  (2l) 
g2  =  p  +  iA  =  n  ,  (22) 


and  let  e(X) ,  E(X) 


X(s-Sn) 

denote  e  and 


X(u-Sq) 

o  du  respectively,  so 


that  in  general. 


s 


o 


e(Xi)  e(7.2)  =  e(Xi  +  Xg)  (23) 

and 

e(X)  =  [e(X)  -  1]A  .  (24) 


In  this  notation  (16)  can  be  written 


=  -  e(-H^)[g2  e(ip)  £(%)  +  gi  e(-ip)  E(H2)] 


=e(-Ho) 


^  e(ip)  +1^  e(-ip) 

LHi  ^2 


^  o(ip)  e(H^+^  a(^)e$y 


L% 


^  e(-Hi)  +  —  e(-H2) 


l_Hi 


Hn 


12  + Si 

Lhi 

«2 

. by  (24) 

by  (18)  -  (20), 
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1.  s. 


2p  e 


©(ip)  +  ef-ip) 

-g2  _^gl 

%  J 

3 

(25) 


Similarly  (17)  becomes 


-B(s) 

-e(ip)  _  e(-ip)-’ 

-1  1- 

e 

Lhi  H2  J 

-h\- 

^i[o-(s)-crJ 


if.4  Prom  2(2),  Z  is  given  by 

Z  -  Zq  =  ^  J  (u)  du  +  I 

C!  J 


du  ; 


(26) 


(27) 


the  term  giving  the  gravity  drop  has  been  omitted.  Equation  (27)  can 
be  written 


Y/here 


Z  -  =  P.2^Zq_(s)  +  1^222(3)  +  Cj^Zj^(s)  +  —  Z-^(s),  (28) 

P 


Pi  3; 

(29) 

^^2  =  ^'/P 

(30) 

s 

Zi(s)  =  1  ^(u)  du 

(31) 

^0 

^2(3)  =  j  du 

(32) 

«.  J 

So 

Zl(s)  =  J['e?L^u)  +  e^°*^^^Jdu 

(33) 

■  ■  i 

s 

^l(s)  =  f 

(34) 
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To  evaluate  (31)  -  (31-)  write 


=  Pi'  +  ikg  =  Pi'  +  ip2 

%  =  Pi'  +  iki  =  Pi'  +  i(p2  -  p) 

=  pi'  +  ik2  =  Pi'  +  i(p2  +  p), 


(35) 


(36) 


(37) 


Then 


2pZi(s)  =  g2E(-Ki)  +  giE(-K2) 


(38) 


2iZ2(s)  =  E(-K^)  -  E(-K2) 


(39) 


2p  Zj^(s) 


f  B(-%) 


E(iY)  (Z*X)) 


L-i_ 

LHi'h2J 


E(iY)  .  (U) 


l*-.5  To  evaluate  Zi(s),  Z}.-(s)  write 


where 


=  (Pl 

-  A)  +  iho  -  ig2  =  go  -  ig2 

(42) 

1 - 1 

n 

CM 

-  A)+  ih^  +  ig^  =  go  +  igl 

(if3) 

H 

II 

H 

-  a)  +  ik^  -  ig2  =  (go  "  "  ^^2 

(44) 

K2  -  (Pi' 

-  A) +  ik^  +  ig^  =  (g^  -  iY)  +  igi  , 

(45) 

So 

=  (Pl'  -  a)  +  i(p2  +  Y)  . 

(46) 

81. 


CONFIDENTIAL 


CONFIDENTIAL 


Technical  Note  No.  G.Vrf’.  177 


The  following  relations  are  readily  obtained  from  (42)  -  (45) > 


62^2  ^  SA  =  ^o^^l  ^  22^ 


(47) 


^^o  ~  ^"^^o  ~  £2_42^  ^°1  ^2^ 


(48) 


Ej_H2  =  8^2  ^  ig^(g^  _  g^)  ^ 


(49) 


A2  ""  "  ^''')(gl  “  S2)  +  §1  £2  ’ 


and  from  (2I)  and  (22)  we  obtain 


^2  =  2P 


(51) 


gl  -  g2  -  “2iA 


(52) 


and 


®1®2  = 


-  V.2 


n2  by  3.1(1) 


(55) 

(54) 


Then  evaluating  (40)  hy  (24)  we  obtain 


2^(s)  ^  (c^  -  Cg)  -  [c^  e(-K^)  +  C2  ©(-Kg)]  +  Cg  e(iY),  (55) 


v/here 


=  °1 


4*  Or 


4_  /  ^2  ^  % 


2p  \l3g_K^  HgKg 


(56) 
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C 


2 


2p  VHi  Hs/J 


/  iY 


(57) 


c.  =± 

2p 


(58) 


- 

^2p 


(59) 


By  (^7)  to  (54)  and  C2  are  found  to  be 


+  iYg^  -  g^2)  /  (60) 

C2  =  (n^  +  iYg^)  /  iY  H^H2  ,  (6I) 

where 

%^2  "  ^  +  So^  By  (49), (52)  &  (54) 

.  (62) 

V2  =  +  2A(g^  -  iY)  +  (g^  -  iY)^  by  (50),(52)  &  (54) 

.  (63) 

and 

-  iY  =  l<^d/A  -  2i8Y-  6  by  (46),(ll),2(6)  &  2(7). 

.  ( 64) 


Likewise  evaluating  (41)  by  (24)  we  have 


-icr, 


°  Zfi(s)  =  (Di  -  D2)  -[d^  e(-Ki)  -  d2  e(-K2)]  +  ^2  e(iY),  (65) 
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where 


D-,  =  d-,  - 


1  “ 


1  1 


2iVHiKi  H2K2 


(66) 


2Y  VHjl  Kg 


(67) 


d.  = 


2i  i 


% 


(68) 


 e 


2i  R 


H2K2 


(69) 


Equations  (66)  and  (67)  can  be  evaluated  in  terms  of  gg,  Y  and  A  as 


=  P-6  [2go  -  iY  +  2A]/Hj_H2K^K2  (70) 

=  i  /  Y  Rqlig  ’ 


v/here  H^Kg,  Kj^^g  are  given  by  (62)  and  (63). 

4.6  .  The  solutions  for  2^(0),  Zg(s)  are,  by  (38),  (39)  and  (24), 


2P  Z3_(s) 


+  gl 


e(-Kg)  - 

^  - 


(72) 


2i  22(3) 


1_ 

Kl  "  K2 


K3(-Ki)  e(-Kg) 

_  Kl  “  K2  _ 


84. 


(73) 


CONFIDENTIAL 


CONFIDENTIAL 


TechnicaJL  Note  No,  G.W,  177 


where  by  (I4),  (45)  and  (51) 


^  +  =  2p(go  -  iy)  /  K1K2 

At  Ao 


1  "^2 


(74) 


(75) 


5  Simplified  solution  for  the  angular  (ieviation 

5*1  Wo  nov/  assume  that 

F.l  The  velocity  is  constant 

F.2  The  gyroscopic  effects  are  negligibly  small,  i.e,  small  spin 
F.3  The  cross-spin  damping  is  negligibly  small  • 

It  then  follows  from  these  assumptions  that 


6^  sin  3,  and  kQd/A 


may  be  put  zero  in  the  preceding  equations:  the  following  sin^lifying 
relations  are  then  found  to  hold. 


A  =  a(s)  =  P*(s)  =  kp/2m  by  2(6)  -  (8) 


=  ai  by  2(9) 


*2  by  2(10) 


=  Pi  by  4.2(7) 


(1) 


=  Pi'  by  4.2(11) 


=  since  is  zero 


=  t/2  , 


and 


p2  =  n2  -  (e/2)2 


P2  =  0 


by  4.2(8) 
by  5.1(1) 
by  4.5(64) 


(2) 

(3) 
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Hence  by  4.5(62)  and  4.5(63) 

Ej^H^  =  (n2  -  Y^)  +  iY^ 

(4) 

=  r? 

(5) 

and  by  2f.5(60),  2f.5(6l)  and  (5) 


0l  =  -  i/HlHj 

(6) 

O2  =  (n^  -  y2)  /  iY  . 

(7) 

By  4.5(70)  and  (71), 


=  p6(e  +  iY)  /  n^ 


(8) 


D2  =  ip.e  /  Y  H^H2 


(9) 


and  by  4,6(74)  and  (75) 


Sp 

=  .0 


(10) 


K. 


=  2ip/n^. 


(11) 


5,2  It  is  seen  from  4.5(55)  that  2t(s)  is  the  sum  of  three  terms; 
a  constant^  a  transient  and  an  oscillatory  term*  The  damping  factor  in 

the  transient  term  is  ^  and  for  sufficiently  large  s-s^ 

the  contribution  from  this  term  is  negligibly  small.  The  locus 
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then  becomes  a  circle  of  radius 
I^lI  1^1  ^2  I  f^om  the  origin. 


IClI  I  Coj  whose  centre  is  at  a  distance 
By  (k),  (6)  and  (?)  we  find  that 


I^lI  IC2I  =  jn^  “  Y^l  /  IyI  Jir?  -  Y^)  ^  (12) 


I^lI  1^1  "  ^2  1  “  “l  V1*^1 


(13) 


In  the  special  case  Y  =  0,  4.5(55)  becomes 


Zl(s)  =  (Ci)y^q  +  (i'^‘^2)Y=0 


(14) 


for  sufficiently  large  s-s^.  Hence 


^L- 


"'Lo 


"•  (Xt 


(15) 


In  this  case  as  s  increases  the  angular  deviation  will  increase 
indefinitely.  The  negative  sign  in  (I3)  appears  on  accoiint  of  the 
v/ay  the  lift  malalignment  is  defined  in  11,3.13. 


5.3  Likewise  the  angular  deviation  locus  due  to  moment  malalignment 
is  seen  by  equation  4.5(65)  to  be  a  circle  for  sufficiently  large  values 


cf  S-Sq. 

distance 

that 


The  radius  of  the  circle  is 

Cm 


"  ^21 


Cm 


from  the  origin. 


lD2land  its  centre  is  at  a 
By  (4),  (8)  and  (9)  yro  find 


p 


D, 


'2  1  =  “ivl 


/  1y|  y (n^  -  Y^)^  +  Y^.C^ 


(16) 


% 

p. 


1^1  -  ^2! 


(17) 


For  a  non-spinning  round  Y  is  zero  and  4.5(65)  hecomes 


pe 


n*" 


n*" 


(s-s  ) 
o 
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The  angular  deviation  is  then  determined  hy 


%  " 


%Io 


if^T.+cr.) 

e  ^  d 


n 


(18) 


5.4-  The  angular  deviation  due  to  an  initial  ya^Y  is,  hy  equations 
4.4(28)  and  4.4(29) 


=  0  for  sufficiently  large  s-Sq 


hy  equations  4.6(72)  and  (lO) .  This  shows  that  the  angular  deviation 
is  entirely  transient. 

5.5  The  angular  deviation  due  to  an  initial  rate  of  turn  of  axis  is 


^  Z2(=) 


by  4.4(28)  and  4.4(30). 


For  sufficiently  large  s-s^,  equations  4.6(73)  and  (ll)  give 


=  p/K^Kj 


hence 


Zj  -  Zju 


n 


V 


(19) 
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APPENDIX  A 

EVALUATION  OF  DTTEGRALS 


1  Fresnel  functions 
1.1  Write 


e(u)  =  C(u)  +  i  S(u)  =  /  e^^  dx 


for  real,  u,  and  define  d(u)  by  the  equation 


-XTrn‘- 

d(u)  =  b(u)  +  i  A(u)  =  e  [(l  +  i)  /  2  -  s(u)] 


(1) 


(2) 


for  all  real  u.  It  can  then  be  shown  that 


(a)  B(u)  and  A(u)  are  monotonic  decreasing  functions  for  positive 
values  of  u  (Refs.  1.3  and  4.3)  and  so  can  be  conveniently 
tabulated. 


(b) 


D(-u) 


e  ^2  ^  (l  +  i)  -  D(u)  =  2  e  ^2  ^  d(o) 


and 


TC  9 


D(-u)  =  e-^2  "  (1  -  i)  -  D(u)  =  2 


-I 


u^,  , 

D(o) 


D(u) 

D(u), 


where  the  bar  denotes  the  complex  conjugate.  Eqxiations  (3)  show 


that  B(u)  and  A(u)  are  oscillatory  for  negative  u. 

(c)  D'(u)  =  -[l  +  i7i:u  D(u)  ]  (4) 

[uD^^“^^(u)  +  (n  -  1)  D^^“^\u)  ]  n  >  1  (5) 

(d)  ujD’(u)j<  |d(u) 1  for  all  finite  positive  u  (6) 

(e)  |d' (u)  j  <  ll’(u)  1  for  all  positive  u  (7) 


1.21  Prom  now  on  we  shall  be  concerned  mainly  with  fuctions  and  integrals 
of  d(u)  where  u  is  of  the  form 


u  =  Xg  =  7^  Vs  +  , 


(8) 


p 

where  ^^o  real  constants,  and  denotes  the  positive  square 

root  of  Xq^.  Then,  X^  is  real  or  imaginary  according  as  X^^  is  real 
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On 

or  imaginary;  and  is  real  or  imaginary  according  as  is  positive 
or  negative. 

V/hen  7.^  is  imaginary,  write 


7g  =  i  7.g'  ,  -(9) 

where 

/s  (10) 


and 


i^l'  =  7-1  .  =  7.^  .  (11) 


1.22  Let  (1)  and  (2)  define  e  (v)  and  D(v)  for  imaginary  arguments'  v 
equal  to  iu.  It  then  follows  that 


s(iu)  =  is (u) 


(12) 


and 


D(iu)  =  iD(u)  . 


(13) 


2  Integrals  occurring:  in  ^‘aw  f^jinctions 

2 , 11  The  integral 


dV^  = 


where 


(1) 


2  a  u  , 


(2) 


can  be  evaluated  in  terms  of  Piresnel  functions. 


92. 


CONFIDENTIAL 


CONFIDENTIAL 


Technical  Note  No.  G,  W,  177 


By  1.21(8),  (1)  and  (2)  give 


.  p 


i^  ^  y  2~’ 

D(X^)e^  -D(X^)e^2*'3 


by  l.l(l)-(2). 


.  (3) 

2.12  We  define  s\  to  be  the  value  of  s  for  which  is  zero,  then 

Vmen  <V^  it  becomes  convenient  to  apply  l.l(3)  to 
(3)  to  obtain  the  alternative  form 


.%  y  2 

XTT 


D(-’^sJ  c  2  -o  _  D(_a  ) 


•  ^  i  2' 


io<  s  <  (4) 


V^o^s)  =  - 


.  Tt  .v  2 
^O 


D(\J  +  (1  +  i) 


So<  s^<s, 


(5) 


2.13  When  is  imaginary,  is  real  aind  is  negative  for  s  <  s^. 
The  alternative  forms  of  equation  (3)  in  this  case  are 


-i-  '2 

HK'  )  e  2  ,)  Q 


•  ^  a  '2' 
-^2\ 


a^<.  Sq  <  s  (6) 


I\(so,s)  =- 


T  ^  "X  ^P  .  ^  -  Ip 

D(-V)  ®  ^  °  e-^2^ 


<  S-x  (7) 


=  - 


\l' 


'2  i  5  ■> '2 

5(-Xs^)  e  2^0  -  d(X3-)  e  +(l-i) 


(8) 
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3  Inte.crals  occurring  in  lispersion  functions 
3.11  yie  nov/  evaluate  the  integral 


^i[a-(u)-<rj 


D(\j)  dV^  =/2a  e  J^Cs^js), 


(1) 


where 


^■(u)  -  0-^  =  n2(u  -  Sq)  +  n^C/u  -  Vs^)  .  (2) 


It  appears  very  unlikely  that  this  integral  can  be  evaluated  exactly 
for  general  values  of  the  parameters.  For  the  present  paragraph  it  is 
supposed  that  is  real  and  positive  in  (s^,s);  the  other  oases  are 
dealt  with  in  paras.  3.21  -  3.26. 

312  Equation  (l)  can  be  expressed  as 


2Y. 


i^  (x-Xo) 


by  the  substitution 


v/here 


Yu  =  Yi  VU  + 

=  2n2A 

YlYo  =  ^lA  , 


(4) 

(5) 

(6) 


so  that 


cr(u) 


). 
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3.13  Caso_£_Yg  positive  for  3  >  i.e,  n2  positive 
Partial  integration  of  (3)  gives 


ioV, 


2Yi  ’ 


R 


(x-Xq) 


p=0 


p+1 


2i 

IT 


dxP 


/x 


as 


00 


,=Y„  2 


provided  Yj^  +0. 


(7) 


This  series  does  not  converge  as  R  ©o,  being  an  asymptotie  series. 
As  far  as  the  seeond  term  equation  (7)  is 


J.(s, 


o? 


3)  = 


2Yl 


2i 

7; 


+  T, 


u 


u 


(8) 


whore  stands  for  the  second  term,  namely 


Yi 


/ 


(9) 


The  condition  that  |tJ  shoiald  be  loss  than  l/N-th  of  |2D(7"^)/7;Y^1  ,  the 
first  tern  in  equation  (8) ,  is  that 


1 


/7rf^2  <  i/N. 


(10) 


By  1.1(7),  (10)  holds  if 


and 


(11) 


Yg^^  -  (/2  1%  Ysg  -  NA  >  0 


95. 


CONFIDENTIAL 


CONFIDENTIAL 


Technical  Note  No.  G.Xi.  177 


3.14 


Case  :  Y 


nep;ative 


for  all  s  ^  i.e.  r\2  negative 


Write 


Y„ 


•2 


and  substitute  y  =  Y  '2  and  y  =  Yo  '2  in  (l) 

J  u  o  ^0  ^  ^  ’ 


-io-  ,  1  -i^(y-yo)  , 

°  =—  i  ,2®  ^  ^^Vy)^ 


(12) 


The  first  tv/o  terms  of  the  asymptotic  expansion  give 


2Y^' 


iQ-(u) 


2i 

7^ 


c(\i) 


Y. 


u 


1  as  Ys^'2.>«., 
^So 


(13) 


v/hero 


T 


! 


u 


Y. 


i  D'  (X  ) 

1 


/  2Y 

u 


(14) 


The  condition  equivalent  to  (ll)  is 


/^>  /s" 

vj  Y 


and 


(15) 


Y3  '2  _  (/2  N\  /r.Y  ')  Yg  '  -  NA  >  0 
*^0  ^o 


3.15  evaluations  of  (3)  given  in  3-13  and  3.14  can  be  eoqiressed 

as  follov/s.  If 


2(u) 


1  G 


ia-{xx) 


/t^YhY. 


I'u 


(16) 
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then,  provided  +  0  and 


1Ys^2i  _  (/2  N?^  /  7;1yJ)  jY^^l  -  N  A  >  0, 


(17) 


the  solution  of  (3)  can  be  taken  as 


J^(so,s)  =  D(X3^)  S(so)  -  0(^3)  2(s)  (18) 

with  an  error  of  less  than  1/N. 

3.16  Tfi/hen  =  0,  it  can  be  shown  that  for  sufficiently  large  n^ 


J^(so,s)  =  D(Xso)  2(^o)  -  ^(’^s)  2(s)  ,  (19) 


where 


2(s) 


.  io'(s) 

1  e 


/  . 


(20) 


The  condition  of  validity  for  (l9) ,  corresponding  to  (l7)  is  that 


|n-]_|  >  /2  NX2_  . 


(21) 


It  is  seen  tliat  (20)  and  (2i)  agree  with  (l6)  and  (l7)  with  Y^^  =  0, 

YiYo  =  niA. 

3.21  So  far  in  this  paragraph,  Xu  has  been  assumed  real  and  positive. 
This  is  case  (a)  below.  The  other  three  cases  have  still  to  be 
considered. 


(a) 

real  and  negative  or 

less 

tlian 

^0 

(b) 

Xu  imaginary  and  Vs^  negative  or 

less 

than  ' 

(c) 

(i) 

^u 

real,  s^^  <  and 

s  <  s« 

k 

(ii) 

real,  s^  <  s^  and 

! 

3 

(d) 

(i) 

imaginary,  s^  <  s^ 

and 

®o  < 

S  <  S, 

(ii) 

K 

imaginary,  s^  <  s^ 

and 

^0  < 
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3 . 22  Case  (b) 

By  1.22(13  )v/G  can  v/rite  (l)  as 

72a  JxCso.a)  =  i  f  (1^.)  dV^  (22) 

and  carry  out  tho  v/ork  of  3.13  -  3.l6  v^ith  D(x^)  replaced  by 
The  result  v/ould  then  be 


S(sq)  -  0(^3')  s(s)] 


(23) 


provided 


l^s 


VF  >  /T, 

o  Y 


(V2  NN^'  /%  IyJ  )  1 


-  N/9t;  >  0 


(24) 


In/hen  n2  =  0,  this  condition  becomes 


j  n^  j  >  V2  N  N  ^ ' 


(25) 


3.23  Case  (c)(i) 

By  1.1(3)  we  can  vrrito  (l)  ac 


V^e  ^  J'^(sq^s)  =2B(o)  e 


0  ^  n  —  Ai^_m  ^ 


f  ^2  ^^u  “^s  ) 


o  dV^-  e 


(26) 


where 


2  2  ..  7..  2  a,  2n  2  ^  2^ 


.u  -  -  ^s,")  -  (V  -  \/)  1 


^  •  (27) 


=  P^U  +  fiQ 
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The  first  integral  on  the  right  hand  side  of  (26)  can  be  evaluated  by 
2,11(1),  and  the  second  integral  by  the  method  of  3.11  to  3. IS,  Vfe 
obtain 


D{-\)  S  (u)  -2D(o) 


(28) 


provided  (17)  holds. 

3. 2if  Case  (c)(ii) 

Put  s  =  s^  in  (28)  and  s^  =  in  (l8)  and  add  - 


=  K(so,3^)  -  [i^(-N3^)  L(so)  +  D(\s)  E(s)]  ,  (29) 


where . 


K(s„,S7,)=2D(o)[2(s^)  +  f^ 


D(|^3  )  e  2  °  -D(p.,  )  o  •" 


\  Q 


(50) 


provided  (17)  holds. 

3.25  Case  (d)(1) 

By  1.1(3)  we  can  write  (22)  as 


VTa  e""^' 


•  ^ 


Xs  '2 


,[cr(u)-cr] 


Jp.:)dV, 


u 


where,  by  (27)  and  1.2l(9), 


(31) 


(32) 
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(31)  can  be  evaluated  to  give 


=  i 


I(u)  -  2  D(o) 

1*1 


i[o-(u)  +  |Xu'^ 


(33) 


provided  (2if)  holds. 
3.26  Case  (d) (ii) 


J^(Sq,s)  =  K(3o,s^)  -  i[D(-?^^')  Z(So)  +  DO.3')  2(s)]  (3^ 

provided  (24)  holds, 

3.3  The  approximate  solution 

3.31  An  approximate  evaluation  of  J)^(so,s)  can  be  obtained  #ien  either 
2  2 

(a)  are  positive  and  is  nearly  equal  to  Y^, 

2  2 

or  (b)  positive  and  is  nearly  equal  to  Y^* . 

In  case:  (a)  never  negative  and  in  (b)  are  never 

positive.  In  (a)  Y^^,  may  be  positive  or  negative,  and  in  (b) 

Y^*^  X^*  may  be  positive  or  negative  according  as  Vs  is  greater  or 
less  than  Vsy,  Vs^. 

(a)  and  (b)  are  divided  into  tv,ro  cases 

(i)  ^3,  alv/ays  positive  i.e.  ^^s^  <  Vs^. 

(ii)  X^  alv^rays  or  sometimes  negative  i.e.  Vs^  ^ 

3*32  case  (a)(i) 

Write 


\  -  ^(1  - 

Yu  =  X  + 

where  are  chosen  so  that 


e  =  1  -Xl/Tj. 

(35) 

■X  ®  “  “  ("''oh.  ’ 

(36) 
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then  (l)  can  be  written 


-icr 
0  ^ 


/  \  1  (x+’^)^-(xq+k)2] 

J7^(so,s)  ®  ^  ■^D(x(l-e))  dx  . 


(37) 


Yoo- 


By  the  mean  value  theorem  applied  to  D(u)  wo  have 


B(x(l-e))  -  b(x)  =  -s  X  B'(x)  +  0(  e^) 
A(x(l-e))  -  A(x)  =  -  e  X  A’(x)  +  O(s^) 


and  hence  that 


1d(x(1-s))  _  D(x)l  <ex  1d*(x)1  +  lo(s2)l 

<  s  |d(x)1  +  l0(s^)  1  (38) 


by  1.1(6)  since  x  is  positive.  • 

If  D(x(l-e) )  is  replaced  by  D(x)  in  (37) j  (38)  shows  that  the 
fractional  error  is  not  greater  than  s  for  s  sufficiently  small. 
liVith  this  approximation  made,  it  is  easily  checked  using  1.1(A)  that 
(37)  can  be  uvaluaxed  to  give 


TCYp 


■  io-(u) 
e 


•i)(Yu)  "  D(Yu“«?0 


(39) 


In  the  special  case  =  0,  equation  (39)  is  exact.  By  (35)  >(36)  wc 
then  have  \  ~  written 


’\(«o»®) 


io(u)  D(Y^)  -  D(A^) 


(AO) 
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In  the  special  case  when  =  0, 


(40)  becomes 


TCY"- 


icr(u) 

o  C'(YJ 


1  •- 


(41) 


In  the  case  k\  ^  0,  +  0,  v/e  see  from  (39)  that  the  solution  is  given 

by  (41),  but  now  only  approximately. 

3.33  Case  (b)(i) 

V/rite 


=  y(i“®4') 

(42) 

V  =  Y  +  ^^x' 

(43) 

where  e^' ,  k^'  are  chosen  so  that 


H  e'  =  1  -  =  e  (14) 


then,  by  1.22(l3),  (l)  can  be  written 


Y  *-K  * 


-icr 

e.  °  Jx(s 


“  Vi'  / 


-i|[  (y+K')^-(yo+K')^]_ 


D(y(l-  e'  ) )  cly.  (46) 


Y  '  -k' 

^o 


Since  y  is  positive  in  the  range  of  integration,  (38)  holds  with  x 
replaced  by  y,  and  we  can  evaliiate  (,46)  approximately  to  give 


■\(so>s) 


1 


9iYi 


"iO-(u)  d(Y^')-d(y^'^k^') 
0  - .  . . . .  . 

^X 


(47) 


with  an  error  not  greater  than  e . 
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By  1.21(11),  1.22(13)  and  (45)  it  is  seen  that  (47)  reduces  to  (39). 
3.34  Case  (a)(ii) 

When  7-2  is  negative  it  is  possible  to  evaluate  (37)  tiy  this  method 
only  when  the  minimum  value  of  i.e.  7^  ,  is  fairly  near  zero.  This 
is  because  the  inequality  ° 


lD(u(l-e))  -  D(u)l  <  e  1d(u)1  ,  (48) 


which  is  used  to  evaluate  (37)  ,  does  not  hold  for  all  negative  arguments. 
Let  u  =  U2(e)  Be  the  value  of  u  (negative)  for  which  the  equality  holds 
in  (48),  then  if  Tqq  >  the  method  of  3-32  is  still  valid.  The  value 

of  u.Q(i)  is  about  -0.5. 

3.35  Case  (b)(ii) 

If  X3  *  >  Uq,  (^4.7)  still  holds  vriLth  a  maximiim  error  e  for  suffi¬ 
ciently  small  e  • 

3,4  approximate  solution  of  J^Csq  ,s)  can  be  obtained  when  7^  is 

small  compared  v/ith  Y^.  It  can  easily  be  verified  that 


Jx(sq,s)  +  ^^d^Cs^^s)  =  - 


E(X„)  D(4„) 


(4-9) 


where  is  given  by  (2/).  'iVhen  is  nearly  equal  to  Yp,  JJ_^(sQ,s)  is 
given  by  (39);  (49)  then  gives  J^(sq,s). 

A  similar  result  holds  when  is  small  compared  with  , 

The  maximum  error  is 


-  or  X^»2/2Y^'2  . 


4  The  Yaw  Integrals  I(;f(sQ,s),  I^(sq,s) 


4.1  These  are  denoted 'by 


/2”a 


%  2 
“so 


in(u-So)  iKu)-o-o] 


dV, 


u 


(1) 
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.7^0  2 

-ip  PSq 

/2iL  e  ^  ° 


^(so,s)  = 


-in(u-So)  i[0<u)-crQ] 

0  e  ’  ^  V . 


(2) 


where 


and 


cr(u)  -  =  n2(u-SQ)  +  (Vu  -  V^) 


=  2au 


(3) 


Write 


and 


I  (“u"  -  -o)  +  -l(^-  -  ^o) 


(0 


2  =  ^^2  ~  n)(u  -  Sq)  +  n3_(Vu  -  V^)  (5) 


=  a  ^  Vu  + 


(6) 


3^  =  ^  /u  +  3o  . 


(7) 


then 


=  2(n2  +  n)  /?; 


«l“o  = 


/  % 


(8) 

(9) 


3^2  -  2(n2  -  n)  A 


(10) 


^^o  = 


(11) 
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If  is  the  value  of  u  for  which  =  0,  then 


^  -  i"‘i/2(n2  +  n)  . 


(12) 


Similarly 


=  -  n^/2(n2  -  n) 


(13) 


4.2  The  evaluation  of  (l)  and  (2)  is,  by  2.1l(3), 


la(so.s)  = 


4.%  „  2 

2 

e  -D(as)  e 


•  „  2-1 


(14) 


1  r  i5.3  2 

4(^=0-=)  =  K-  ^  °  -i:(Ps) 


2-, 


(15) 


When  a  3  are  negative  or  imaginary,  (14)  and  (15)  can  be  transformed 
by  2.1(^)  -2.1(8)  SO  that  all  arguments  are  positive. 

5  The  integrals  Jct(sQ,  s) ,  J3(so,  s)  ' 

5.1  These  integrals  are 


V2a  e  °  J£x(sq,s) 


i[o-(u)-o-J 


D(au)  dV^  (1) 


So 


V2a  e  ^  o  J3(sq,s) 


i[o-(u)-<rJ 


E(3u)  ^^u»  (2) 


where  cr(u)  and  are  given  by  4.l(3)^  and  and  3^  by  4.1(6)  -  4.l(ll) 
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5.2  Conditions  for  the  validity  of  the  Asymptotic  Solution 

The  conditions  under  which  the  asymptotic  solution  holds  with  an 
error  of  less  than  1/N  are  given  by  3..15(l7)  ot"  3.22 (2if);  they  can  both 
be  written  as 


o  Y 


lYs/l  -  (/2  NjN  1  A  1y.  1)  lYg  1  -  nA  >  0 


(3) 


where  now  stands  for  either  or  3;[_. 
By  3.12(1.)  to  3.12(6)  we  have 


lYgl  =  A  1  +  n^  /  2x12) 


(4) 


when  s  >  Sy 


and 


^  ^  /  2n2  . 


(5) 


By  11.9.3(13)  -  (II-)  and  11.9.2(3)  -  (k) ,  equations  (k)  and  (5)  give 


Ys  1  =  jr^l  /  H 


(6) 


Vi;:  -  VT 


=  “  t'o  /  n2 


(7) 


Equation  (3)  then  requires  that 


ryn2>  0 


(8) 


and 


n2  Va  jn2jNa 


(9) 
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If  we  choose  to  be  positive,  then  (8)  requires  n2  to  be  positive  and 
(9)  gives 


2r^  /Van2  ^  1  ^ 


(10) 


Of  the  two  values  for  |7\2.  1j  '/2(n2+n) /%  and  ^2  |n2-n  |/'^,  the  former  is  the 
larger  since  n  is  positive;  substituting  it  into  (lO)  we  have 


1 


In  the  particular  case  when  n/n2  is  small  and  N  is  large,  (ll)  gives 


2  a  n2 


(12) 


There  has  been  no  loss  of  generality  in  taking  r  to  be  positive: 
if  it  were  assumed  negative^  (8)  and  (9)  would  lead  to  the  same  equation 
as  (11)  with  r^  and  n^  replaced  by  -r^  and  -n^. 

5.3  The  Asymptotic  Solution  -  Expressions  for  ^ * 

ViOien  the  conditions  of  5*2  hold  the  evaluation  of  for 

X  =  a  or  3,  is  given  by  the  results  of  3.11  to  3.26.  These  solutions 
will  now  be  written  out.  There  are  four  cases  (see  3.2l). 

5.31  Case  (a)  ;  ot^2  (3^2)  positive,  and  (or  Pg)  positive 
for  all  s  ^  s^. 

From  3.15  (l6)  ,  3.15(18)  and  3.12  (4)-(6) ,  wc  have 


J7.(so,s)  =  -  1 


D(\i) 


i[oi;u)-o-  ] 


2n2  Vu  +  np 


(13) 


for  7,  =  a  or  ^  . 

5.32  Case  (b)  !  (or  P. ^)  negative,  and  ot  '  (or  p  ')  positive 


for  all  s  >  Sq. 


By  3.22(23) 


-icr, 


D(X^') 


for  X  =  a  or  3. 


2n2  Vu  +  n^ 
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5.33  In  cases  (c)  and  (d)  the  results  can  be  somewhat  simplified. 
Prom  3.23(27)  we  have 


Pl  - 


(13) 


^'l^^o  =  '^I'^o  -  H’^o 


(16) 


and  hence  from  3.12(5)  ~  (6)  said  4.1(8)  -  (ll)  equations  (15)  and  (16) 
beconie 


=  -  2n/7^  (X  =  a) 


2-0/7.  (X  =  p) 


(17) 


=  0  ( X  =  a  or  8  )  . 


(18) 


If  we  write 


=  +  ^2n/% 


(19) 


63  = 


(20) 


and 


2(f)  =  2''(=o.s)  =  (2i/6l) 


15(63  )  -  B(6  )  o 

L.  O 


in( s-Sq) 


(21) 


we  find  using  1.22(13),  and  (17)  -  (2l)  that 


% 


-ins„  — 


D(Hg)  e^^^®  e  °  (s)  /  2i  ( X  =  a)  (22) 


=  e^""°^^(s)  /2i  (X  =  3)  .  (23) 
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Hence  since  cr(u)  -  /  2  is  independent  of  u 


i[cr(u)-|\u'"] 


(-i/2)  S^s) 


.  r 


(\  =  a)  (24) 


=  (i/2) 


.r_  7^  Q  2, 

//  N  V2  ^O  ] 
4  (,s)  e 


(4  =  3)  (25) 


and  in  the  particular  case  Virhen  s  = 


1  1*^  u  2 


l-p_  2, 

0  2^s;\.  ] 


l[ots^)-|Ps^] 


>(1/2)^  (s„)  e 


p  _  T 

^  ^0~2 


(4“a)  (26) 


=(-l/2)S'‘(sp) 


(4-  3).  (27) 


5.34  Gaso(c)  (i)  ;  (or  3^^^)  positive  and  ttg  (or  3g)  negative. 

Prom  3.23  (28)  and  (24),  (25) 


-10"^ 

Q  Ja(so^s)  =  i 


^(-Oa) 


i[o-(u)-crJ 


2n2  /u 


+  n-i 


+  iD(o)  ?’(s)  (28) 


°■J^(so.s)  =  i 


D(-e„) 


2n2  +  np 


.  7C  3_^  2 

-  i  D(o)  ^''(s)  e""*^  .  (29) 
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Case  (c)(ii);  (or  positive  and  Og  (or  positive. 

Prom  3.24(30)  and  (26) ,  (27)  we  have 


-iCTp 


K(j 


o> 


i[cr(scj)-o-o] 

ST)  =  2D(o)  — 

n  Hi 


■il/(n2+n) 


-I2  «£ 


(T  =  a)  (30) 


=  -  2D(o) 


.  i[o-(sQ)-o-J 
1  e  ^ 

n  n^/(n2-n) 


+  I  ^  (sp) 


J 


(7^=  P)  ,(31) 


and  from  3.22f(29) 


3"^°°  Jx(3„,3)  =  K(; 


o> 


3^- 


P(-N.o) 

2n2''^+np  2n2'^s+np 


D(X  )  e 
s 


i[o<s)-cro] 


(32) 


for  T.  =  a  or  p.  •  ■■ 

5.35  Case  (d) (i) ;  a-^  (or  negative  and  (or  P^)  negative. 

Prom  3.25(33)  and  (24-),  (25) 


o  Ja(So,s)  = 


D(-V) 


i[o-(u)-o'Q] 


2x1'^ u+n^ 


a.  »2 


D(o)  ^(s)  e 


(33) 


-ICTn 


i  [o<u)-o*o] 


2n2'^+ni 


X/  N  "2 


+  D(o)  S  (s)  e 


i:?Ps'2 


(3if) 
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Case  (d)(ii);  (or  negative  and  (or  Pg)  positive. 

Prom  3.26(34) 


Jx(sq,s)  =  e  °  K(sq,s^.)  + 


F(X 


_ 

2n2'''^+ni  2n2Vs+np 


(35) 


for  \  =  aor  P,  where  e  K(sq,s^)  is  given  by  (30)  or  (3I). 

5.4  Asymptotic  Solution  in  Constant  Spin  Case 

5.41  When  riy  is  zero,  the  condition  corresponding  to  (9)  is, 
by  3.16(21)  and  3.?2(25), 


jn^l  >  V2  N  14^1  (36) 

where  |7t^[  now  equals  J2n/K ,  and  n^  equals  %/2/a  r^. 

It  should  be  noted  that  when  4^  stands  for  and  n^  is  positive, 
or  when  ?tg  stands  for  Pg  and  n-j^  is  negative,  the  error  is  usually  much 
less  than  i/Nj  this  is  shown  in  the  foUov^ing. 

By  A.  1(8)  -  (9) 


>  n^/V2n  -K 


since  is  supposed  positive 


>  V2  N/t: 


(37) 


by  (36). 

For  Values  of  of  this  size  (i.e.  > l)  the  asymptotic  expansion  of 

D(u)  holds  with  an  error  of  less  than  10^,  and  gives  (see  Ref.  I.4) 


I  r'(u)  /  D(u)|  =  l/u  .  (38) 
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Hence  the  error  in  the  asymptotic  solution  3'.l6(l9),  which  in  generaJ.  is 
the  left  hand  side  of  equation  3.13(lO)>  is  seen  in  this  case  to  be 


D'(0 


E(“s) 


,  which  is  less  than  %/  2  N^  by  (36),  (37)  aui  (38);  i.e. 


the  error  is  very  much  less  than  1/N. 

Similarly  when  stands  for  and  n^  is  negative  we  have 


by  4.1(10)  -  (11), 


^1 


>  -  n^  /  'J2n% 


since  nq  is  supposed  negative. 

By  the  same  argument  as  above  the  asymptotic  solution  3.22(23) 
has  an  error  less  than  %/2  N‘^. 

In  the  remaining  cases:  v/hen  stands  for  and  nq  is  negative, 
and  when  \s  stands  for  and  nq  is  positive,  we  can  say  that 


and 


-  Oo  >  72  N/t^,  when  n-,  <  0,  provided  /s  «  -n-|/2n  = 

X  X  nr 


-  p^'  >  7’2  N/t^,  when  nq  >0,  provided  Vs  «  n^/2n  =  7s^. 


Hence  only  when  7s  is  very  much  lass  than  7s^  (or  Tso)  will  the  error 
be  as  small  as  7^/2  N^.  As  s  approaches  s^^  (or  s^  the  error  increases 
to  its  maximum  value  1/N,  since  at  s  =  s_  (or  so)  the  equality  sign  holds 
in  1.1(7).  As  s  increases  above  2sq^  {or  2sp),  (or  P^)  exceeds 
72  N/t^  and  the  error  becomes  less  than  V2  N^  again.  The  error  in  the 
asymptotic  solution  is  therefore  less  than  7^/2  N^  except  when  s^  and  s 
are  near  s^j  (or  s^ ,  where  the  error  is  1/N.  It  is  usually  found  that  _q 
v/hen  (36)  holds  /^<<Vs^  (or  V^)  ;  for  example,  if  N  =  5  aind  Ttn  =  0,02ft, 
(or  Sp)  ^  1200  ft. 

5.42  Prom  4.1(8)  and  (lO)  it  is  seen  that  aq^  is  positive  and 
p]_^  is  negative;  therefore  if  rQ  is  chosen  in  a  positive  sense  nq  is 
positive,  and  hence  is  positive  for  all  s  >  Sq,  v/hile  Pg'  is  negative 
for  s  <  s  <  ^  and  positive  for  s  >  Sq.  Hence  J^(sq,s)  falls  under 
case  (a)  and  J^(sq,s)  under  case  (d) .  By  (13)  we  have 


i[o-(u)-o-J‘ 


(39) 
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and  by  (-34)  and  (35) 


S  yX  '2 

+  D(n)S  (s)  e  5*^0 


(40) 


L  ^  J  P  ^ 


(41) 


5*5  The  Approximate  Solution 

5.51  This  solution  holds  in  the  case  when  or  4^'Ai' 

(Xl  =  or  near  unity.  The  maximuiii  error  in  the  solution  is 

then 


lel=  ll  -  /Yil  =  ll  -Xi«  /Yi'  1  . 


By  4.1(8)  and  (lO),  and  3.12(5),  Xj/Y]^  ~Ji+  n/n2  in  the  two  cases, 

showing  that  this  solution  holds  when  n/n2  is  sinall.  The  maximum 
error  is  then 


je  I  =?  n/2n2  . 


It  should  be  noted  that  the  error  will  usually  be  much  less  than  this  value. 
The  reason  is  that  when  the  integrands  of  J^(3,Sq),  Jp(s,SQ)  are  largest 
(3^  will  be  near  zero  -  v/here  the  approximation  to  D(x(l-e))  in 
para.  3,32is  most  accurate.  The  error  in  the  integrands  only  approaches 
|e|  when  tend  to  inf  ini  ty  by  l.l(6);  i.e,  when  the  integrands 

of  J(x(s,Sq),  Jp^s,Sq)  approach  zero. 

There  is  a  further  condition  which  is  required  to  hold  when 
13^  or  Pg*  *are  negative.  In  3.3^  this  is  discussed  as  case  (ii). 

In  most  practiceil  examples  it  is  almost  certain  to  be  satisfied.  The 
condition  is  that  etc,,  should  not  be  *  too  negative*,  but  greater 
than  about  -0.5.  V/hen  n2  is  positive  it  is  seen  from  i4-.l(8)  to  (ll) , 

11.9.3(13)  -  (14)  and  II,  9.2(3)'-(4)  tliat  oc^  and  are  greater  than  -0.5 
provided 


2nVs^  +  V27a  r^  >  -- J'K  (112  +  n)/2 

(42) 

-2n/s^  +  'fz/a.  >  -jT^{n2  -  n)/2  . 

(43) 
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If  Tq  is  zero  or  positive  (42)  will  always  be  satisfied.  When  r^  is 
zero  (this  is  the  particular  case  for  vfhich  this  method  is  used  to  obtain 
the  results  of  Pt.l),  (4-3)  req^uires  that 


s 


< 

o 


(Mf) 


Hence  if  n  =  0.01  ft  ,  s  =  0.1  (44-)  requires  that  s^  <  200  ft.  V/hen  n2 
is  negative,  the  required  conditions  are  those  obtained  by  replacing 
rQ,  n2  by  -r^,  -n2  respectively  in  (42)  and  (43). 

5.52  Choosing  n2  positive,  the  solutions  for  Jgj(s,SQ),  Jq(s,  Sq) 
are  given  by  3.32(39).  It  is  found  that 


i(n2  -I-  n) 
n 


:  [d(Yu)  -  d(Yu  -  Kcc)  ] 


^ifr(u)-.c^  “I 


(4-5) 


"i°'o  /  V  4(^0  -  n)  r  r  ,  .  ,  ST  ircr(u)-ov,]  (  , 

e  °  J3(sq,s)  = - ? -  [  D(Yj  -  D(y^  -  Kp)]  e  *■  ^  \  ,(46) 

“  n,  4  >3^5 


n 


where  by  3.12(4)  -  (6)  and  3.32(36) 


Y^  =  (2n2  Vu  -f  n^)  /  '^2%n2 


(47) 


Kgf^  =  n  /  V2  7tn2  (n2  +  n) 


(48) 


Kp  =  -  n  n^^  /  V 2  xn2  (n^  -  n) 


(49) 


5.53  If  n  is  negligible  in  comparison  ¥n.th.  equations 
4.1(8)  -  (11)  and  3.12(4-)  -  (6)  show  that 


(50) 


hence  Jjj(s^,s)  and  Jp(s^,s)  are  equal  and  independent  of  n. 
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Prom  equations  (47)  to  (49)  it  is  found  that  both  -  Kj^and 
Y^  -  Kp  tend  to  Y^  as  n/n2  tends  to  zero,  so  that  in  the  limit  equations 
(45)  and  (46)  become 


i 


i[o-(u)-o'J 
e  ^ 


D'(’^u) 


(X  =  a,  p). 


(51) 


5.54  If  the  launching  spin  is  Sufficiently  large  P(Y^)  can  be 

replaced  by  the  first  term  of  its  asymptotic  expansion  (Ref.  1.4).  The 
required  condition  is 


7CY3  2  » 1 


i. 


e. 


by  5.2(6). 
Then 


and  from  (5I) 


4 


a  nr 


>  > 


1 


D(r  J  =  iAy^ 


=  -  3 


i[o‘(u)-o-o] 


Y, 


u 


(X  =  a,  3). 


(52) 


In  this  case  however,  v/hen  both  n2/n  and  Va  n2  are  very  muoh 

greater  than  unity,  the  asmptotic  holds  in  view  of  the  condition  5.2(l2)* 
By  equations  (50),  3.15(l6;  and  3«15(l8)  we  obtain 


-ICr 


A.  =  " 


% 


Yl 


^i[cr(u)-c7j  D(f  u) 

Y 


u 


(X  =  a, 3), 


which  agrees  v;-ith  (52)  asynpitotioally. 
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KEY  FIG.  I 

T  a  t/-ll  G  -  BOOST  COUPLE 

’g  «  -dG/C  C  a  M.  OF  I.  IN  ROLL 

■ll*a2c/airR  a  a  acceleration 

r  a  SPIN  Tr  *  SPIN  DAMPING  CONSTANT 

t  -  TIME 


A  a  SPIN  PRODUCED  BY  OFF-SET  NOZZLES 
B-SPIN  PRODUCED  BY  SPIN  AT  LAUNCH 
C  a  A  -t-  B 


FIG.I  SPIN -FORMS. 

(plotted  non- dimensionally) 
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FIG.  2 
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FIG. 2  DISPERSION  FUNCTIONS  FOR  BOOST 
DESTABILISING  COUPLE.  ZERO  LAUNCHING 
SPIN.  CONSTANT  ANGULAR  ACCELERATION. 
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FIG.  3 
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•  FIG.  3  DISPERSION  FUNCTIONS  FOR  BOOST 
DESTABILISING  COUPLE.  ZERO  LAUNCHING 
SPIN.  CONSTANT  ANGULAR  ACCELERATION. 
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FIG.4 
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So  »  EFFECTIVE  LAUNCHER  LENGTH 
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FIG.4  DISPERSION  FUNCTIONS  FOR  BOOST 
DESTABILISING  COUPLE.  ZERO  LAUNCHING 
SPIN.  CONSTANT  ANGULAR  ACCELERATION. 
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FIG.  5 


FIG.5  DISPERSION  FUNCTIONS  FOR  BOOST 
•  DESTABILISING  COUPLE.  ZERO  LAUNCHING 
SPIN.  CONSTANT  ANGULAR  ACCELERATION. 
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FIG.  6 
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FIG.6  DISPERSION  FUNCTIONS  FOR  BOOST 
DESTABILISING  COUPLE.  ZERO  LAUNCHING 
SPIN.  CONSTANT  ANGULAR  ACCELERATION. 
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FIG.7&8 


FIG.  7 
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FIG.8 

FIG.7&8  DISPERSION  FUNCTIONS  FOR  BOOST 
DESTABILISING  COUPLE.  ZERO  LAUNCHING 
SPIN.  CONSTANT  ANGULAR  ACCELERATION. 
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FIG.  lO. 
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FIG.IO.  YAW  LOCUS  OF  MOTION  IN  FIG.  9. 
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FIG  II 


FIQI.DISPERSION  FUNCTIONS  FOR  TRANSVERSE 
BOOST  FORCE,  ZERO  LAUNCHING  SPIN, 
CONSTANT  ANGULAR  ACCELERATION. 
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FIG  12 
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FIG  12,  DISPERSION  FUNCTIONS  FOR  BOOST 
DESTABILISING  COUPLE,  CONSTANT  SPIN. 
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FIGI3. 


FIG  13.  DISPERSION  FUNCTIONS  FOR  BOOST 
DESTABILISING  COUPLE.  CONSTANT  SPIN. 
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FIG.I4&IS 


DESTABILISING  COUPLE.  CONSTANT  SPIN. 
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TVi  «  Jzla.  »  1-5 
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FIG.  16  DISPERSION 
FUNCTION  LOCUS. 
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FIG.l6ai7 
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FIG.  17  YAW  LOCUS. 


FIG.I6&I7  CONSTANT  SPIN.  CONSTANT  BOOST 
DESTABILISING  COUPLE  WHOSE  AXIS  IS 
INITIALLY  IN  DIRECTION  OF  IMAGINARY  AXIS. 
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FIG.  18.  FACTOR  R  BY  WHICH  DISPERSION  OF  UNSPUN  ROUND 

ANGULAR  ACCELERATION 

IS  REDUCED  AGAINST -a LINEAR  AeCELl^TlC>Kl . 
ZERO  LAUNCHING  SPIN. 
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FIG. 19.  FACTOR  R  BY  WHICH  DISPERSION  OF  UNSPUN  ROUND  IS 


REDUCED  AGAINST 


To 
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FIG.20.  MAXIMUM  DISPERSION  OF  DART  DUE  TO  LIFT 
MALALINEMENT  AGAINST  SPIN  (PLOTTED  NON  -  DIMENSIONALLY.) 
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FIG.  22. 
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